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Preface 


The aim of this book is, as its title suggests, fo help someone with little or ne 
knowledge of what thermal analysis can do, to find out briely what the subject 
is all about, to decide whether it will be of use to him or her. and to help in 
getting started on the more common techniques. Some of the less-common 
techniques are mentioned, but more specialized texts should be consulted 
before venturing into these areas. 

This book arose out ofa set of notes prepared for courses on thermal analysis 
given at instrument work shops organized by the $.A. Chemical Institute, It has 
also been useful for similar short courses given various universities and 
technikons. I have made extensive use of the manufacturers’ literature. and Lam 
grateful to them for this information. A wide variety of applications has been 
drawn from the literature to use as examples and these are acknowledged in the 
text. A fuller list of the books, reviews and other literature of thermal analysis is 
given towards the back of this book. The ICTA booklet ‘For Retter Thermal 
Analysis’ is also a valuable source of information. 

Tam particularly grateful to my wife, Cindy, for typing the manuscript. to 
Mrs Heather Wilson for the line drawings, and to Professor David Dollimore of 
the University of Toledo, Ohio, for nuiny helpful suggestions. 


Michael E. Brown 
Grahamstown 1987 


Thermal analysis 


is the measurement of changes in physical 
properties of a substance as a function of 
temperature whilst the substance is subjected to a 
controlled temperature programme 


Chapter 1 


Introduction 


tion and history 


What do bread and chocolate, hair and finger-nail clippings. coal and rubber, 
ointments and Suppositories, explosives, kidney stones and ancient Feyptian 
papyri have in common? Many interesting answers could probably be 
Suggested, but the connection wanted in this context is that they all undergo 
interesting and practically important changes on heating. 

The study of the effect of heat on materials obviously has a long history, from 
man’s earliest attempts at producing pottery, extracting metals ( ~ 8000 ne) and 
making glass (~3400 rnc) through the philosophical discussions of the 
alchemists on the elements of fire, air, earth and water, to carly work on the 
assaying of minerals (~ 1500 Ap), followed by (he development of thermometry 
and calorimetry [1, 2]. Only in the late 19th century did experiments on the 
effect of heat on materials become more controlled and more quantitative. 
Much of this work depended upen the development of the analvtical balance 
which has its own interesting history [3, 4], 

The establishment of the International Confederation of Thermal Analysis 
(ICTA) in Aberdeen in 1965 and the great advances in commercially available 
equipment have resulted in thermal analysis being an extremely active field with 
applications in numerous directions 

There has been some discussion [ 5, 6] on whether thermal analysis can really 
be regarded as an analytical method. It was concluded [5] that thermal analysis 
could be fitted into the general systems theory model of chemical analysis, based 
on three ‘elements’: the sample. the reagent and the signal, by considering heat 
as a reagent (Fig. 1.1(:)). 

It has also been suggested [ ] that the definition of thermal analysis should be 
extended to allow for rapid heating of the sample to some elevated temperature, 
followed by measurement of the Propecty with time under isothermal 
conditions, i.c. at zero heating rate. 
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Figure 1.1 (a) Systems theory model of thermal analysis. (b) Generalized thermal 
analysis instrument, (c) Thermal analysis curve (note: ‘thermogram’ is not an approved 
term), (Based on Levy. P. F, (1971) Int. Lab, Jan/Feb. 61, with the permission of 
International Scientific Communications, Inc.) 


The main techniques 1 
$.2 Thermal analysis instruments 


All thermal analysis instruments have features in common. These are illustrated 
in Fig. 1.1(b). The variety of the techniques to be discussed stems from the 
variety of physical properties that can be measured and the variety of 
transducers that can be used to convert these propertics into electrical signals. 

Measurements are usually continuous and the heating rate is often, but not 
necessarily, linear with time. The result of such measurements is a thermal 
analysis curve (Fig, 1.1(c)) and the features of this curve (peaks, discontinuities. 
changes of slope, etc.) are related fo thermal events in the sample The thermal 
events which may be detected are described in chapter 2. 


1.3 Types of measurement 
Absolute values of the sample property may be recorded, or the difference in the 
property of the sample compared to the same property of a reference material 
may be more convenient to measure, or the rate of change of the sample 
property with temperature may be of interest (derivative measurements) 
1.4 The main techniques 


The main sample properties used and the associated techniques are listed in 
Table 1.1. All of these techniques may be combined with evolved gas analysis 


Table 1.1 The main therma 


alysis techniques 


Property Technique Abbreviation 
Mass Thermogravimetry* TG 
Derivative thermogravimetry DIG 
Temperature Differential thermal analysis DTA 
Enthalpy Differential scanning calorimetry Dsc 
Dimensions Thermodilatometry 
Mechanical properties Thermomechanical analysis TMA 
(Thermomechanometry) 
Dynamic mechanical analysis DMA 
Optical properties Thermoptometry or thermomicroscopy 
Magnetic properties Thermomagnetometry ™ 
Electrical properties Thermoelectrometry 
Acoustic properties Thermosonimetry and 
thermoacoustimetry TS 


Evolution of 
radioactive gas Emanation thermal analysis ETA 
Evolution of particles’ Thermoparticulate analyst TPA 


*The term ‘thermogravirhetric analysis’ and the ahbreviation "TGA' are in common use but are not 
approved by the ICTA Nomenclature Committee 
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(EGA) (see chapter 10), and it is often possible to carry out simultaneous 
measurements of more than one property (see chapter 9). 


ab, =, 
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Chapter 2 


Thermal events 


2.1 Reactions of solids 


When a single pure solid substance, A, is heated in an inert atmosphere, the 
resultant increase in molecular, atomic or ionic motion may lead to changes in 
crystal structure, sintering, melting or sublimation [1,2]. Ifthe intramolecular 
forces are weaker than the intermolecular forces, the substance may decompose 
forming new molecular fragments, some or all of which may be volatile at the 
temperatures reached, Numerous examples of such decompositions exist [3,4] 


e.g. 
BaCl, . 2H,O(s)-»BaCl,(s) 4+ 2H,O0(g) 


CaCO(s)+CaO(s) + CO,(g) 
NH,Cl(s)- NH, (2) + ΟῚ) 
BaN,(s) Ba(s) + 3N,(g) 
More complicated reactions result when the initial solid can react with the 
surrounding atmosphere [5]. e.g. 
2Ag(s) + 30,(2)-Ag,O(s) 
CuO(s) 4 15 (g)-Cu(s) + H,O(g) 
Ni(s)+ 4CO(g)-+NiiCO)4(2) 
C(s) FO(g) + CO,(p) 


When more than one solid substance is present initially, there are 
correspondingly more possibilities for interaction on heating and new phases. 
such as solid solutions and cutectic mixtures, may form as well as new 
compounds formed by addition or double decomposition reactions [6]. ¢.2 


Fe,O,(s)4 MgO(s) *MpFe,O,(s) 
» NaC lis) 4+ AgRr(s)- + AgClis) + NaBr(s) 


The above changes are nearly always accompanied by enthalpy changes, and 
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sometimes also by changes in mass, so they may be studied using one or more of 
the thermal analysis techniques listed in Table 1.1. The main thermal events are 
summarized in Table 2.1, 


Table 2.1 Thermal events 


A(s,) ᾿ Als;) phase transition 
. Atl) melting 
‘ Alp) sublimation 
ὺ R(s) 4 gases thermal 
decomposition 
. gases radiolytic 
A(glass) »A(rubber) glass transition 
Als) ) Bop) + C(s) oxidation 
tarnishing 
Als) | Big) μῆκος f combustion 
ὶ volatilization 
Ais} | (pases), Als) 4 (gases), heterogencous 
catalysis 
A(s) 1 Ris) + ABs) addition 
ΑἸ. | CD(s) +AD(s) + CRIs) double 
decomposition 
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Chapter 3 


Thermogravimetry (TG) 


3.1 Introduction 


Measurements of changes in sample mass with temperature are made using a 
thermobalance. This is a combination of a suitable electronic microbalance 
with a furnace and associated temperature programmer. The balance should be 
ina suitably enclosed system so that the atmosphere can be controlled (Fig. 3.1 
and section 3.4). 


furnace 


sample 


temperature 
programmer 


yy 


Vigure 3.1 A schematic thermobalance. 


3.2 The balance 


Several types of balance mechanism are possible [1-3]. These include beam, 
spring, cantilever and torsion balances. Some operate on measurements of 
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deflection, while others operate in null mode, Null-point weighing mechanisms 
are favoured in TG as they ensure that the sample remains in the same zone of 
the furnace irrespective of changes in mass. 

Various sensors have been used to detect deviations of the balance beam from 
the null-position, ¢.g. the Cahn RG clectrobalance (Fig, 3.2) uses an electro- 
optical device with a shutter attached to the balance beam. The shutter partly 
blocks the light path between a lamp and a photocell. Movement of the beam 
alters the light intensity on the photocell and the amplified output from the 
photocell is used to restore the balance to its null-point and, at the same time, is 
a measure of the mass change. The restoring mechanism is electromagnetic. The 
beam has ribbon suspension and a small coil at the fulcrum, located in the field 
of a permanent magnet. The coil exerts a restoring force on the beam 
proportional to the current from the photocell. 

Provision is also usually made for electrical taring and for scale expansion to 
give an output of mass loss as a percentage of the original sample mass. 

The sensitivity of a thermobalance and the maximum load which it can 
accept, without damage, are related, Typical values are maximum loads of | g 
and sensitivities of the order of | ug. Use can be made of alternative suspension 
points on the balance beam to increase the maximum load. 

The output signal may be differentiated electronically to give a derivative 
thermogravimetric (DTG) curve. 

In addition to conventional TG studies, much use has been made of 
microbalances in studying the kinetics of reactions of solids (see chapter 13) 
under isothermal conditions [4]. Studies of decomposition reactions pose some 
special problems. Many decompositions are reversible if a supply of the gaseous 
products of reaction is maintained, Such studies thus require careful control of 
the surrounding atmosphere and should include runs at reduced pressures. 
Reducing the pressure, though, worsens the heat exchange, causing problems in 
temperature measurement, In an attempt to overcome some of these problems, 
the use of quartz crystal microbalances has been suggested. 

Use of the piezoelectric effect in certain crystals (usually quartz) for 
measuring the mass of material deposited or condensed ona crystal face is well 
documented [5, 6]. There are basically two ways in which such crystals can be 
used in TG studies. The sample may be deposited on a crystal face and the 
sample plus crystal may be heated as required, or the sample may he heated 
separately in one part of a reaction chamber and all or some of its gaseous 
decomposition products may be condensed on the initially clean face ofa crystal 
which is held at a suitably low temperature. Changes in the amount of material 
deposited on the crystal surface show up as changes of frequency of oscillation of 
the crystal, which is usually excited in a conventional series resonance circuit, 
The observed change in frequency depends on the value of the frequency itself 
and the mass and area of the coating on the crystal face. Mass changes of as little 
as 10) '? g can be detected, 

Quartz crystals also show marked changes in resonance frequency with 
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Figure 3.2 Flectronic microbalance with null mechanism, (With the permission of Cahn 
Instruments, Ine.) 


temperature. so most studies have been limited to isothermal conditions. The 
Maximum temperature al which quartz exhibits pyroclectric behaviour is 
57 ©. Lithium niobate crystals can be used to extend this range to beyond 
1000, 

Henderson etal. [7] have used a quartz crystal for non-isothermal TG by 
using a computer to record and correct for the frequency temperature and the 
mass frequency relationships of the crystal to obtain the TG (mass -tempera- 
ture) curve, Descriptions of the TG system, including arrangements for heating 
the crystal plus sample, and of the hardware and software required are given 
[7]. The system was used to record TG curves of several coordination 
compounds, Samples (~ 1 jig) were applied to the quartz crystal face as 
solutions in dichloromethane, 

The less direct approach, that of determining the mass of gaseous 
decomposition products. condensed on a cooled quartz crystal [8,9] avoids the 
frequency temperature relationship problem and also removes some of the 
problems of heating and measuring the temperature of the sample. The sample 
can be heated in a conventional furnace (with contact temperature measure- 
ment) in a system where the pressure can be reduced as desired. Usually the 
difference in frequency between the ‘weighing’ crystal and a matched reference 
crystal (in a separate chamber) is amplified and recorded. The temperature of 
the weighing crystal is also recorded. For quantitative studies the system has to 
be calibrated for different gaseous products and sample geometry. 

A microbalance based on the mass dependence of the resonance frequency of 
a vibrating quartz wire, from which loads are suspended, has been described by 
Mulder [10] 


3.3 Heating the sample 


In the more conventional thermobalances, there are three main variations in the 
position of the sample relative to the furnace as shown in Fig. 3.3. The furnace is 
normally an clectrical resistive heater and may also, as shown, be within the 
balance housing, part of the housing, or external to the housing. It should: (i) be 
non-inductively wound, (ii) be capable of reaching 100 to 200°C above the 
maximunt desired working temperature, (ili) have a uniform hot-zone of 
reasonable length, (iv) reach the required starting temperature as quickly as 
possible (Le havea low heat capacity), and (v) not affect the balance mechanism 
through radiation or convection, Transfer of heat to the balance mechanism 
should be minimized by inclusion of radiation shields and convection baMes 

Furnaces containing electrical resistive heaters [11] rely mainly on heating of 
the sample by conduction, through solid or gas, with inevitable large 
temperature gradients, especially when dealing with samples of low thermal 
conductivity such as polymers and inorganic glasses. Heating by radiation 
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Figure 3.3 Alternative arrangements of furnace. sample and housing. 


becomes significant only at high temperatures in such furnaces, but alternative 
heating systems using either infrared or microwave radiation have been 
considered. 

For infrared heating [12. 13] the light from several halogen lamps is focused 
onto the sample by means of clliptic [12] or parabolic [13] reflectors 
Temperatures of over 1400°C may be achieved at heating rates of up te 
1000°C min! and control within +0.5°C. 

Karmazin ef al. [14,15] have suggested the use Gf microwaves to generate 
heat uniformly within the sample. Such a process would have the immediate 
advantage of allowing for the use of larger and more representative samples 
than are often used in thermal analysis. Temperature measurement and power 
control both present problems in using microwaves. A thermocouple may be 
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positioned in the waveguide if its orientation is such that it collects no energy 
itself from the microwaves, and its size is small compared to the guide so that 
reflected waves are negligible, The thermocouple is then used [14] for power 
control ria a microcomputer system. 

For DTA (chapter 4) the sample and the reference are placed in cylindrical 
Teflon or silica containers, The absorption coefficient for microwaves varies 
from one substance to another and also varies markedly with temperature, thus 
making matching of the sample with a suitable reference difficult. As with 
classical DTA, the sample is often mixed with an excess of inert reference 
material, such as ALO, or TIO, which have high absorption coeflicients, The 
higher the absorption coefficient, the greater the heating rate possible. 

Hydrates show interesting behaviour on heating with microwaves since the 
absorption coeflicient for HO is so high (e.g. ~ 1000 times that of CaSO, in 
gypsum), Conventional DTA shows dehydration as an endothermic process, 
but with microwave heating the process appears to be exothermic. 

Similarly. it was found [15] that ina microwave-heated thermodilatometer 
(chapter 6) the results differed considerably from those obtained using 
conventional heating. 

The possibility has also been suggested [16] of the use of lasers for heating 
and infrared pyrometers for remote temperature measurement, in the in situ 
thermal analysis of bulk materials. 

Steinheil [17] described a vacuum microbalance system, for use at high 
temperatures (1600 to 2400°C) which incorporated induction heating by a high- 
frequency electromagnetic field, but did not allow weighing during actual . 
heating because of the electromagnetic forces acting on the sample. 


3.4 The atmosphere 


Thermohalances are normally housed in glass or metal systems [18] (Figs 3.1 
and 3.4), to allow for operation at Pressures ranging from high vacuum 
(<10 * Pa) to high pressure (> 3000 kPa), of inert, oxidizing, reducing or 
corrosive pases. Microbalances are affected by a variety of disturbances [1,18]. 
One correction which may have to be made is for buoyancy arising from lack of 
symmetry in the weighing system, Ifthe asymmetry can be measured in terms of 
a volume AV, the mass of displaced gas (assuming ideal gas behaviour) is 
m= pMAV/RT (where pis the pressure and M the molar mass). The buoyancy 
thus depends not only on the asymmetry AV, but also on the pressure, 
femperature and nature of the gas. Attempts may be made to reduce Alto zero, 
ora correction may be calculated, or an empirical correction may be applied by 
heating an inert sample under similar conditions to those to be used in the study 
of the sample of interest. 

At low pressures (10? to 270 Pa) [1],a particular problem is thermomolecu- 
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Figure 3.4 Mettler TA3000 ~The TG50 Thermobalance. (With the permission of 
Mettler Instrumente AG.) 


lar flow [1] which results when there is a temperature gradient along the sample 
holder and support. This gradient causes ‘streaming’ of molecules in the 
direction hot-+cold, ic. up the suspension as a rule, giving spurious mass 
changes. Thermomolecular flow may be minimized (i) by working outside the 
Pressure range by adding inert gas, or (ii) by careful furnace design and sample 
placement, including use ofa symmetrical balance design with twin furnaces, or 
(iii) by determination of corrections required, using an inert sample, as is done 
for buoyancy corrections. 

The sample may be heated in a small container with a restricted opening. 
Decomposition then occurs in a self-generated atmosphere [19] of gaseous 
decomposition products, The inhibiting or catalytic effects of these products on 
the decomposition can then be studied 

Garn and co-workers [20, 21] have described the design and use of furnaces 
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for carrying out reactions in controlled atmospheres of ligand. Such systems are 
essential for studying the thermodynamic reversibility of dissociations of 
coordination compounds. 

At atmospheric pressure, the atmosphere can be static or flowing. A flowing 
atmosphere has the advantages that it: (i) reduces condensation of reaction 
products on cooler parts of the weighing mechanism; (ii) flushes out corrosive 
products; (iii) reduces secondary reactions: and (iv) acts as a coolant for the 
balance mechanism. The balance mechanism should, however, not be disturbed 
by the gas flow. It is possible for the balance mechanism to be protected by an 
inert gas atmosphere, while a corrosive gas or vapour is passed over the sample 
(ον. 11,0 vapour for dehydration studies) 

The noise level of Ὁ traces at pressures above ~ 20 kPa usually increases as 
the temperature increases on account of thermal convection [ 1, 22, 231. The use 
of dense carrier gases at high pressures in hot zones with large temperature 
gradients gives the most noise. Variations in the flow rate of the gas do not affect 
the noise level much, but may shift the weighing zero. Noise levels also increase 
as the radius of the hangdown tube increases. Thermal convection, and hence 
noise, can he reduced by altering the gas density gradient by introducing a low 
density gas, such as helium, above the hot region, Alternatively, and more 
practically, baffles can be introduced in the hangdown tube. A series of close- 
fitting convoluted baMes (Fig. 3.5) was found to be most successful [24]. Even 
with the baffle it was found that changes in ambient temperature caused non- 
turbulent gas movements in the hangdown tube and hence apparent changes of 
Mass, 


at 


Figure 3.8 Reduction of convection effects by use of bales in the hangdown tube 


Temperature measurement and calibration 1s 


Sample containers should be of low mass and made of inert material (e.g 
platinum containers may catalyze some reactions). Samples should generally be 
thinly spread to allow for ready removal of evolved gases. Results should be 
checked for sample-holder geometry effects. 

Build up of electrostatic charge on glass housings is another source of 
spurious mass effects. Methods proposed for dealing with this problem include 
use of weak radioactive sources for ionization of the balance atmosphere. 
coating of glassware with a sputtered metal film or other metal shielding. and 
the use of various commercial sprays or solutions. Wiping the outside of the 
glassware with ethanol is reasonably effective 


3.5 The sample 


Although solid samples may be nominally of (he same chemical composition. 
there may be considerable differences in their behaviour on heating, These 
differences arise from structural differences in the solid, such as the defect 
content, the porosity and the surface properties, which are dependent on the 
way in which the sample is prepared and treated after preparation. Forexample. 
very different behaviour will generally be observed for single crystals compared 
to finely ground powders of the same compound [25]. In addition to the 
influence of defects on reactivity [26], the thermal properties of powders differ 
markedly from those of the bulk material. 

As the amount of sample used increases, several problems arise. The 
temperature of the sample becomes non-uniform through slow heat transfer 
and through either self-heating or self-cooling as reaction occurs. Also exchange 
of gas with the surrounding atmosphere is reduced. These factors may lead to 
irreproducibility. Even when the sample material is inhomogeneous and hence 
a larger sample becomes desirable (e.g. coal and mineral samples). the sample 
mass should be kept to a minimum and replicates examined for reproducibility 
if necessary. Small sample masses also protect the apparatus in the event of 
explosion or deflagration. The sample should be powdered where possible and 
spread thinly and uniformly in the container [25]. 


3.6 Temperature measurement and calibration 


The sample temperature, 7,, will usually lag behind the furnace temperature, 7; 
and 7, cannot be measured very readily without interfering with the weighing 
process. The tag, 7;—7,, may be as much as 30°C, depending upon the 
operating conditions. The lag is marked when operating in vacuum or in fast- 
flowing atmospheres and with high heating rates, Temperature measurement is 
usually by thermocouple and it is advisable to have separate thermocouples for 
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measurement of 7) and for furnace regulation. (Platinum resistance ther- 
mometers are used in some controllers.) 

An ingenious method of temperature calibration [27] for small furnaces 
makes use of the Curie points of a range of metals and alloys. On heating a 
ferromagnetic material, it loses its ferromagnetism at a characteristic tempera- 
ture known as the Curie point. Ifa magnet is positioned below the sample of 
ferromagnetic material, as shown in Fig. 3.6(a), the total downward force on the 
sample, at temperatures below the Curie point, is the sum of the sample weight 
and the magnetic force, At the Curie point the magnetic force is reduced to zero 
and an apparent mass loss is observed (Fig. 3.6(b)), By using several’ 
ferromagnetic materials, a multi-point temperature calibration may be 
obtained (Fig. 3.6(c)). It has also been suggested that inert mass-pieces be 
suspended from the balance, using links of fusible wire so that disturbances will 
be caused as the fusion temperature is reached [3133]. 
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Figure 3.6 Curie-point method of temperature calibration [27]. (With the permission of 
Thermochionea Acta, Plsevier, Amsterdam.) 
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Even with careful temperature calibration, 7, may still not be accurately 
known, as slow heat transfer may cause sel-heating or self-cooling from 
strongly exothermic or endothermic processes in relatively large samples. 


3.7 Temperature programmers 


A wide variety of fully-electronic temperature controllers and programmers is 
available commercially, Temperature sensors are cither platinum resistance 
thermometers or thermocouples, Heating rates offered range from fractions ofa 
degree per minute to nearly 1000°C min ᾿ς with facilities for isothermal 
operation. Ina small thermal analysis laboratory, one programmer may serve 
several instruments [28]. 


3.8 Interpretation of TG and DTG curves 


Actual TG curves obtained may be classified into various types [29] as 
illustrated in Fig. 3.7. Possible interpretations are as follows. 

Type (i) curves. The sample undergoes no decomposition with loss of volatile 
products over the temperature range shown. No information is obtained, 
however, on whether solid phase transitions, melting, polymerization or other 
reactions involving no volatile products have occurred, Use of some of the other 
techniques (Table 1.1) is necessary to climinate these possibilities. Assuming 
that these possibilities are eliminated, the sample would then be known to be 
stable over the temperature range considered. This could be good news ifa heat 
resistant material was being sought, or bad news if potential explosives were 
being tested! 

Type (ii) curves. The rapid initial mass loss observed is characteristic of 
desorption or drying. It could also arise, when working at reduced pressures, 
from effects such as thermomolecular flow or convection (section 3.4), To check 
that the mass loss is real, it is advisable to rerun the sample, which should then 
produce a type (i) curve, unless the carrier gas contained moisture or was very 
readily readsorbed on the sample at the lower temperature. 

Type (iii) curves. These represent decomposition of the sa mple in a single stage, 
The curve may be used to define the limits of stability of the reactant, to 
determine the stoichiometry of the reaction, and to investigate the kinetics of 
reaction (see chapter 13). 

Type (iv) curves, These indicate multi-stage decomposition with relatively 
stable intermediates. Again. the temperature limits of stability of the reactant 
and of the intermediates can be determined from the curve, together with the 
more complicated stoichiometry of reaction 

Type (0) curves, These also represcat multtstage decomposition, but in this 
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Figure 3.7 Main types of thermogravimetric (TG) curves. (Based on Duval, C. (1963) 
Inorganic Thermoaravimetric Analysis, Elsevier, Amsterdam, 2nd edn. and Daniels, T. 
(199%) Thermal Analysis. Kogan Page, London, with permission.) 


example stable intermediates are not formed and little information on all but the 
stoichiometry of the overall reaction can be obtained, It is important to check 
the eflect of heating rate on the resolution of such curves. At lower heating rates, 
type (v) curves may tend to resemble type (iv) curves more closely, while at high 
heating rates both type (iv) and type (v) curves may resemble type (iii) curves 
and hence the detail of the complex decomposition is lost, 
Type (vi) curves, Vhese show a gain in mass as a result of reaction of the sample 
with the surrounding atmosphere. A typical example would be the oxidation of 
a metal sample. 
Type (rii) curves, These are not often encountered. The product of an oxidation 
reaction decomposes again at higher temperatures (e.g. 2Ag+50,+Ag,0— 
2Ag+30,) 

Resolution of stages of more complex TG curves can be improved by 
recording DTG curves (Fig. 3.8). Such curves may also be produced from 
digital TG data by numerical differentiation. 
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Figure 3.8 Comparison of TG and DTG curves. 


3.9 Applications of TG 


Application of TG is limited. to some extent, in that not all ofthe thermal events 
listed in chapter 2 are accompanied by changes in mass. For desorption, 
decomposition and oxidation processes, however, much valuable information 
can be obtained from TG alone. Much of the earlier work in TG was on the 
accurate definition of conditions for drying or ignition of analytical precipitates 
[7]. Examples of TG curves for CuSO, . 517,00 and for CaSO, . 2,0 are 
given in Figs 3.9 and 3.10, The mass losses define the stages, and the conditions 
of temperature (and surrounding atmosphere) necessary for preparation of the 
anhydrous compounds, or intermediate hydrates, can be established immedi 
ately. At higher temperatures, the sulphates will decompose further. Knowledge 
of the thermal stability range of materials provides information on problems 
such as the hazards of storing explosives, the shelf-life of drugs and the 
conditions for drying tohacco and other crops. By using an atmosphere of air or 
oxygen, the conditions under which oxidation of metals and degradation of 
polymers become catastrophic can be determined, 

TG curves for more complex materials, such as minerals and polymers, are 
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Figure 3.9 TG curve for CuSO, .5H,O. 
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Figure 3.10 TG curve for CaSO, .2H,0 at different water-vapour pressures. (Wiede- 
mano, ΤΙ, Οὐ, (1977) J Thermal Anal., 12, 147, with the permission of Wiley Heyden 
Ltd.) 


not always immediately interpretable in terms of the exact reactions occurring. 
Such curves can, however, be used for ‘fingerprint’ purposes. The ΤῸ curve 
obtained on a given apparatus, under specified conditions, is compared with a 
bank of reference curves accumulated on the apparatus concerned. Some sets of 
reference curves have been published [30] but, as mentioned earlier, 
comparison is best when curves from the same instrument are available. 

The reactions corresponding to the mass losses can best be determined, or 
confirmed, by simultaneous evolved gas analysis (EGA) (chapter 9). For 
example, in Fig. 3.9, the appearance of traces of SO,,SO, and O, inthe evolved 
gases would indicate the onset of sulphate decomposition. A complementary 
technique, such as hot-stage microscopy (HSM) (chapter 5) may provide 
information on the mechanism of dehydration or decomposition, by showing 
up the formation and growth of decomposition nuclei, or progress of the 
reactant/product interface inwards from crystal surfaces. 

When the process occurring is clearly defined, c.g. the stoichiometric 
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dehydration of a definite hydrate, the kinetics of the reaction can be determined 
from the TG curves (or from a series of isothermal curves of mass loss against 
time, obtained at different temperatures). Details of some of the kinetic analyses 
that have been suggested are given in chapter 13. Values of activation energies. 
obtained in this way, have been used to extrapolate to conditions of very slow 
reaction at low temperatures (in predicting shelf-lives of materials, resistance to 
weathering. and in estimating rates of natural processes, e.g. petroleum genesis 
over geological times) and to very fast reaction at high temperatures (behaviour 
of propellants and explosives), 

The thermobalance may also be used to measure the vapour pressure of a 
sample, e.g. a metal, by determining the rate of mass loss through a calibrated 
orifice in a Knudsen cell. 

It is worth noting, too, that the balance itself may be used, with little or no 
modification, for magnetic susceptibility measurements (chapter 11), accurate 
density measurements, determination of surface areas by adsorption, and 
particle-size analysis by sedimentation. 
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Chapter 4 


Differential thermal analysis (DTA) 
and differential scanning calorimetry 


(DSC). 


4.1 Classical DTA [1, 2] 


DTA is the simplest and most widely used thermal analysis technique. The 
difference in temperature, A7, between the sample and a reference material is 
recorded while both are subjected to the same heating programme. In ‘classical’ 
DTA instruments, represented schematically in Fig. 4.1, a single block with 
symmetrical cavities for the sample and reference is heated in the furnace, The 
block is chosen to act asa suitable heat-sink, and a sample holder of low thermal 
conductivity is included between the block and the sample to ensure an 
adequate differential temperature signal during a thermal event 

Should an endothermic thermal event (AH positive, such as melting) eccur in 
the sample, the temperature of the sample. 7,. willlag behind the temperature of 
the reference, 7), which follows the heating programme. If the output from the 
thermocouples, ΔΊ᾽ -- 7, -- 7,, is recorded against 7; (or the furnace temperature, 
T,~ T,), the result will be similar to Pig. 4.1(c). [fan exothermic process (A/T 
negative, such as oxidation) occurs in the sample, the response will he in the 
opposite direction, Since the definition of ATas T,— 7, is rather arbitrary, each 
DTA curve should be marked with either the endo or exo direction, The 
negative peak, shown in Fig. 4.1(c), is called an endotherm and is characterized 
by its onset temperature. The temperature at which the recorder response is at 
its maximum distance from the baseline, 7),,,. is offen reported but is very 
dependent upon the heating rate. ¢), used in the temperature programme and 
factors such as sample size and thermocouple position. The area under the 
endotherm (or exotherm) is related to the value of the enthalpy change, A//, for 
the thermal event (for a more detailed interpretation, see section 4.7) 

The reference material should have the following characteristics: 

(i) Itshould undergo no thermalevents over the operating temperature range 


24 Differential thermal analysis and scanning calorimetry 


)6.. A T = 
2 
LZ huarmace 
(a) 


AT ᾿ 
͵ 

τῆ | 
! ! 
| a! eee 
ἥρω Vines 

pe 
(c) 


Figure 4.1 Differential thermal analysis (DTA). (a) Classical apparatus (S= sample: 
R = teference); (b) heat-flux; (c) typical DTA curve. (Note the DTA convention that 
endothermic responses are represented as negative, i.e. downward peaks.) 


(ii) It should not react with the sample holder or thermocouple. 
(iii) Both the thermal conductivity and the heat capacity of the reference should 
be similar to those of the sample. 
Alumina, Al,O,, and carborundum, SiC, have been extensively used as 
reference substances for inorganic samples, while for organic compounds. 
especially polymers, use has been made of octyl phthalate and silicone oil, 
Both solid samples and reference materials are usually used in powdered 
form, The particle size and the packing conditions influence results. A common 
technique for matching the thermal properties of the sample to those of the 
reference, 1s to use the reference material as a diluent for the sample. There must 
obviously be no reaction of the sample with the reference material. 
The furnace system is usually purged with an inert gas and the possibilities of 
afmosphere control are similar to these diseussed for TG (chapter 3) 
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4.2 Calorimetric DTA or heat-flux DSC [3] 


In ‘calorimetric’ DTA, also known as Boersma DTA or differential dynamic 
calorimetry, the sample and reference, in similar holders, usually flat pans, are 
placed on individual thermally conducting bases. The thermocouple junctions 
are attached to these bases and are thus not directly in the sample or reference 
material. This configuration has the advantage that the output signal is less 
dependent upon the thermal properties of the sample (see above), but response 
is slower (Fig, 4.1(b)). 

The temperature range of DTA depends on the materials used for the furnace 
windings and for the thermocouples. 


4.3 Differential scanning calorimetry (DSC) [4, 5] 


In power-compensated DSC, the sample and a reference material are 
maintained at the same temperature (AT 7, ~1,=0) throughout the con 
trolled temperature programme. Any energy difference in the independent 
supplies to the sample and the reference is then recorded against the programme 
temperature. The apparatus is shown schematically in Fig. 4.2(a), and an 
example of a resulting DSC curve in Fig. 4.2(b). The Perkin-Elmer DSC-2 is 
shown in Fig. 4.3 and a block diagram in Fig. 4.4. There are many similarities 
between DSC and DTA, including the superficial appearance of the thermal 
analysis curves obtained, but the Principle of power-compensated DSC is 
distinctly different to that of calorimetric DTA, also known as heat-fux DSC 
(section 4.2). 

Thermal events in the sample thus appear as deviations from the DSC 
baseline, in either an endothermic or exothermic direction, depending upon 
whether more or less energy has to be supplied to the sample relative to the 
reference material. Again. these directions should be clearly marked on the 
record, to avoid later confusion. In DSC, endothermic responses are usually 
represented as being positive, ic. above the baseline, corresponding to an 
increased transfer of heat to the sample compared to the reference. Uinfortun- 
ately this is exactly (he opposite convention to that usually used in DTA, where 
endothermic responses are represented as negative temperature differences, 
below the baseline, as the sample temperature lags behind the temperature of 
the reference. 

The operating temperature range of power-compensated DSC instruments is 
generally more restricted than that of DTA instruments. The Perhin-Flme: 
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Figure 4.2 Differential scanning calorimetry (DSC). (a) Apparatus (S=sample; 
R=reference), (b) typical DSC curve. (Note the DSC convention, which is opposite to 
the DTA convention, that an endothermic response is represented by a positive, i.e. 
upward peak.) 
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Figure 4.3 Differential scanning calorimeter, (With the permission of Perkin-Flmer 
Corporation.) 
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Figure 4.4 Model DSC-2 block diagram. (Watson, FE. S.. O'Neill, Μ. J. Justin, 1. and 
Brenner, N. (1964) Anal. Chem 36,1233, with the permission of the American Chemical 


Society.) 


Table 4.1 Recommended standard materials for calibration 


(i) Temperature calibration 


Transition Transition 
Material temperature (C) Material temperature ( (ἢ 
KNO, 127.7 K,SO, 582 
KCIO, 299.5 K,Cro, 665 
Ag,SO, 412 BacO, 810 
SiO, (quartz) §73 srco, 925 


(b) Fathalpy (and temperature) calibration 


Material Mpc) ΜιριΚ) AH pening μι"). AH να κα") 
Indium 156.4 4296 6.80 28.5 

Tin 2319 505.1 14.404-0.01 60.254 004 
Lead 127 4 600.6 5.454001 22.501 00 
Zinc 419.8 692.7 25.9 4-01 1084-04 


Aluminium 664) 4 346 OAD 1 97.} 


Taken from ref, 4 
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4.4 Comparison of the principles of DTA and DSC (uGe ἢ 
A schematic diagram of a differential thermal apparatus is given in Fig. 4,5 


together with definitions of the terms needed for comparison of classical DTA, 
heat-flux DSC and power-compensated DSC. 


Th > temperature of heat source 


Ts ~ actual sample temperature Tr = actualreference temperature 
Tsm > measured sample temperature Trym = measured reference temperature 
Cs ~ heatcapacity of sample + pan Cr ~ heat capacity of reference 4 pan 
Csm “ heatcapacityofmonitoring station Crm - heat capacity of monitoring station 


π΄ RR, and R’, are thermal resistances 
The da/dt terms are heat flows 


Figure 4.5 Differential thermal apparatus. (From Mraw, S.C. (1982) Rev. Sei, Inst., $3, 
228, with the permission of the American Institute of Physics ) 


Foran ideal instrument, the heat capacities and thermal resistances would be 
natched, i.e. C= Cys R= R,= R and Ri = Rl Ξ Κ΄. Note that Καὶ τὶ Κ΄ and 
“7 C,. Wis further assumed that C,>C, and that heat flow is governed by 
Newton's law 


dq/dt=(1/R\AT 


Heat Now to the sample side heats both (i) the sample monitoring station and 
ii) the sample 


Ag. dt = Cygld Ty, /dt) + (470 46) 
Iso 


dq./dt= CAA T,/dt) 
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830 
dq,/dt=C,,,(07,,,/dt) + dqt/dt 
Applying Newton's law 
dq./dt =(1/R)(Ty~ ΤῸ} 
and 
dqudi=(W RYT, T) 


Similar expressions hold for the reference side, 


Classical DTA 
Thermocouples are in the sample and in the reference material so that pe) 
and 7,,,=7,, i.e. R'=0. 
The signal then is 
AT= RIA IMC, - (ἢ 
and depends upon the difference in heat capacities, the heating rate and the 


thermal resistance, R. R is difficult to determine as it depends on both the 
instrument and the properties of the sample and the reference 


Power-compensated DSC 


The power is varied to make ΤῸ το 7,,,. Thus Ty Ty = Tm and R=O. ie. there 
is no thermal resistance. The signal is then Διά 4η -- (47. ) ((ς ©) 


Heat-flux DSC 
The conditions are that 
Ty, # Tem τ ΤΙ and R#R' £0. 
The signal 
ATaT, 


τὴν 


which is of similar form to classical DTA, except that R depends only on the 
instrument and not on the characteristics of the sample. 

Risa function of temperature and the relationship between Rand 7 has to be 
determined by calibration of the instrument (section 4.6) at several different 
temperatures. Some manufacturers of heat-lux DSC instruments have 
incorporated electronic compensation for this temperature dependence and 


Tay = RAT IANN(C, (Ὁ) 
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conversion of the instrument output to units of power, in the instrument control 


module. 


4.5 Sample containers and sampling 


At temperatures below 500°C (773 K) samples are usually contained in 
aluminium sample pans, One type of pan has lids which may be crimped 
(Fig. 4.6) into position, while for volatile samples, pans and press are available 
which enables a cold-welded seal, capable of withstanding 2 3.atm 
(220. 300 kPa) pressure, to be made. It is worth noting that use of aluminium 
pans above 500°C will result in destruction of the DSC sample holder (several 
thousand pounds!). It is essential to ensure that the instrument is fitted with a 
protective circuit, limiting the maximum temperature which may be reached 
without operator intervention, and to inform the operators of the danger. 

For temperatures above 500°C, or for samples which react with the Al pans. 
gold or graphite pans are available (but expensive!). Glass sample holders have 
heen described in the literature but are not available commercially. The 
reference material in most DSC applications is simply an empty sample pan. 
The sample-holder assembly is purged by a gas which may be inert or reactive as 
desired. (The high thermal conductivity of helium makes it unsuitable for 
thermal measurements although it has advantages for evolved gas analysis.) 
Low temperature attachments are available for extending the operating range 
of the DSC to as low as -- 175°C, 

Dense powders or discs, cut out of films with a cork borer, are ideal samples. 
Low density powders, flocks or fibres may be prepacked in a small piece οἵ 
degreased aluminium foil to compress them, A small syringe is used for filling 
pans with liquid samples. When volatile products are to be formed on heating 
the sample. the pan lid should be pierced. This is also necessary when reaction 
with the purge gas is being investigated. 

It is recommended that the total mass of the sample + pan + lid be recorded 
hefore and after a run so that further deductions on the processes occurring can 
he made from any change in mass. 


4.6 Quantitative aspects of DTA and DSC curves 


The features of interest in DTA or DSC curves are the deviations of the signal 
from the baseline. and the baseline is not always easy to establish, Initial 
displacements of the baseline itself from zero result from mismatching of the 
thermal properties of the sample and the reference material, and asymmetry in 
the construction of sample and reference holders, In severe cases this may cause 


Figure 4.6 Crimping and welding presses. (With the permission of Perkin Filmer 
Corporation and Mettler fastraments AG.) 
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a sloping baseline, which may require electrical compensation, After a thermal 
event, the response will not return to the original baseline if the thermal 
Properties of the high-temperature form of the sample are different from those of 
the low-temperature form. Many procedures for bascline construction have 
been supgested. Some of the simpler approximations are illustrated in Vig. 4.7, 


Vigure 4.7 Simpler procedures for baseline extrapolation 


MN should be immediately obvious whether a DSC or DTA feature isa peak or 
a discontinnity, provided that the baseline can be established with certainty 
Whether a peak is in the exothermic or endothermic direction should be 
ascertamed by comparison with a known melting (endothermic) peak, 
Although this seems obvious, reversal of polarity on the recorder switch or after 
femiporiry disconnection of the recorder, can cause a lot of confusion. 

ΠῚ there is any suggestion from initial runs that any of the features involve 
overlapping peaks. then further experiments should be done in which the 
sample mass and the heating rate are varied to see whether the individual peaks 
can be resolved, 

Abrupt changes in slope or position of the baseline usually indicate second- 
order transitions, Fxamples of this lype of transition are the glass transition in 
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polymers (Fig. 4.8) and the Curie point transition in ferromagnetic materials 
(chapter 3), The enthalpy change. A//, for such transitions is zero, but there isa 
change in heat capacity. 

Once a satisfactory bascline has been defined, the area of the endotherm or 
exotherm is determined by: (i) counting of squares: (ii) cutting out and 
weighing; (iii) planimetry: (iv) digitization and numerical integration: (v) 
mechanical or electronic integration. The measured area, 4. is then assumed to 
be proportional to the enthalpy change, AH, for the thermal event represented: 


AH=Ax K/m 


where m is the sample mass and A is the calibration factor, This factor has to be 
determined by relating a known enthalpy change to a measured peak area. 
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Figure 4.8 Glass transition of a printed cireuit hoard (hrominated bisphenol A resin), 
The DSC curve (proportional to the specific heat) is compared with the curve of the 
linear expansion coefficient (proportional to the first derivative of the TMA curve) 
(From Mettler Application No. 3402. with permission ) 
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Usnally the melting ofa pure metal, such as indium, is used for calibration. (The 
melting point ofthe standard is used, at the same time. to check the temperature 
calibration of the instrument.) 

The calibration factor contains contributions from the geometry and thermal 
conductivity of the sample/reference system and is thus specific to a particular 
instrument under one set of operating conditions. For DTA or heat-flux DSC, 
the value of K is markedly dependent upon temperature, so that calibrations 
have to be carried out over the full operating temperature range of the 
instrument, Electronic compensation for the temperature dependence of the 
calibration factor is built into some heat-Nux DSC instruments. A significant 
advantage of power-compensated DSC over DTA is that the proportionality 
constant, K. is virtually independent of temperature, and calibration is thus 
much easier. 

Table 4,1 lists the transition temperatures and enthalpies of materials 
recommended as standards. A range of higher-temperature standards is likely 
to be available soon. through ICTA. Ideally, the thermal properties of the 
standard including its transition temperature, should be as close as possible to 
those of the sample under consideration. 

The procedure for the determination of the enthalpy change for a process 
occurring in the sample, A//,, is then to determine 


R= MF ucioneativraniics X MASS,/A 
and to use this value in 
AH, =(KS. RAMS Rom) 
where 4 is the area and subscripts c and s refer to calibrant and sample. 
respectively; δ᾽ is the recorder chart speed and R the sensitivity range of the 
instrument/recorder, Note that the heating rate used in the temperature 


programme does not enter into the calculations, but. in precise work. should be 
the same for sample and calibrant. 


4.7 Interpretation of DTA and DSC curves 


At the heart of all thermal analysis experiments lies the problem of correlating 
the features recorded with the thermal events taking place in the sample. Some 
aspeets of this correlation for TG results have been discussed in section 3.8, 
DSC or DTA provides different information and if simultaneous measurements 
(chapter 9) are not possible, the results of parallel experiments, using different 
techniques such as TG and DSC or DTA, under conditions (e.g. sample mass. 
heating rate, atmosphere) as closely matched as possible, are most valuable. 
Once the main features ofthe DTA or DSC curve have been established, and 
hascline discontinuities have been examined, attention can he directed to the 
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correlation of the endothermic or exothermic peaks with thermal events in the 
sample. 

A useful procedure is to test whether the event being monitored is readily 
reversible on cooling and reheating, or not. Exothermic processes are not 
usually readily reversible, if at all. in contrast to melting and many solid: solid 
transitions. 

The melting endotherms for pure substances are very sharp (i.e. they occur 
over a narrow temperature interval) and the melting point. 7),,. is usually 
determined, as shown in Fig. 4.2(b), by extrapolation of the steeply rising, 
approximately linear, region of the endotherm back to the baseline. For impure 
substances the endotherms are broader and it is possible to estimate the 
impurity content (up to a maximum of about 3 mol%) from the detailed shape 
of the melting endotherm (chapter 14). 

TG and EGD or EGA information is invaluable at this stage in distinguishing 
between irreversible or slowly-reversible phase transitions and decompositions. 
Quantitative mass losses and detailed FGA may, in addition, lead to 
determination of the stoichiometry of decomposition, Account must naturally 
be taken of the gaseous atmosphere surrounding the heated sample 
Dehydration may be found to be reversible on cooling in a moist atmosphere 
before reheating, and carbonate decompositions are usually reversible in CO, 
atmospheres. Comparison of features observed, under otherwise similar 
conditions, in inert and oxidizing atmospheres is valuable. TG results may show 
increases in mass corresponding to reaction of the sample with the surrounding 
atmosphere. 

It cannot be over emphasized that as many additional techniques available. 
such as hot-stage microscopy (chapter 5), conventional elemental analysis, X- 
ray diffraction (XRD) and the many types of spectroscopy. should he used to 
confirm suggested interpretations of TA features recorded. Table 4.2 summar- 
izes the overall interpretation procedure, but this should not he followed 
blindly. There is no substitute for experience built up by running samples with 
well-documented behaviour on your own instrument, and also examining the 
almost legal-type of ‘proof of events provided by more conscientious authors 
(spurred on by more demanding editors). 


4.8 Measurement of heat capacity 


DSC can be used to measure heat capacities of materials. As heat capacity 
values are fundamental to thermodynamics. this is an important feature, 
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In the absence of a sample, i.c. with empty pans in both holders of Fig. 4.208}. 
the DSC baseline should be a horizontal line. A sloping baseline may be 
observed when the sample and reference holders have different emissivities, 1.e. 
if the amount of energy lost by radiation does not vary in the same way with 7, 
for both sample and reference. This could occur for a black sample and a shiny 
reference, so a sloping baseline is often an indication that the sample-holder 
assembly should he cleaned by heating in air to high temperatures, with 
aluminium pans removed, 

When a sample is introduced, the DSC baseline is displaced in the 
endothermic direction and the displacement, ᾿ς is proportional to the total heat 
capacity, C,. of the sample. (The total heat capacity = mass ¥ specilic heat 
capacity.) This is illustrated in Fig. 4.9(a) 


hee, or h= BAC, 
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pan 
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sample 
pan only 
T——<> 
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dt 
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C, liquid 
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Figure 4.9 Determination of heat capacities using DSC, (a) Displacement of baseline, 
(b) scan of sample compared to standard, (6} shift in baseline after transition, 
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re “is the heating rate and B the calibration factor, The value of B is 
romned using a standard substance, ¢.g. sapphire, scanned under similar 
fitions fo the sample (Pig. 4,9(b)). Although the technique is simple in 
ciple, Suzuki and Wunderlich [8] have shown that a great deal of care is 
ssary Co get aecurate and reproducible results 

cease of the proportionality of the DSC response to the heat capacity of the 
ple. a shift in baseline after a transition is common (Pip, 4.9(c)) and the 
Tine has to he estimated as shown in Fig. 4.7 or by more elaborate 


οτος [9] 


4.9 Measurement of thermal conductivity 


rmal conductivity, 7, is an important physical property, particularly with 
present emphasis on the efficient use of energy. Numerous methods for 
suring thermal conductivity have been published and dedicated commer- 
instruments are available, Several papers have been published (10 13] in 
Ἵν heat-flux DSC instruments have been used to measure thermal 
luctivity. An advantage of using DSC is that specific heat capacities. ¢,,can 
ieasured in the same instrument (section 4.8) and hence, if density values, p, 
wailable or can be measured, thermal diffusivities, D, can be calculated 
ithe relationship, D= 4/pe,. 

hin and Pair | 10} have described an attachment for the DuPont 990 DSC 
(Pig 4.10), which allows for measurement of thermal conductivity of a 
ve (in cylindrical form) without modification of the basic instrument. The 
rerature at the bottom of the sample (7) is measured ria the output of the 
‘while the temperature at the top of the sample (7,) is measured with a 
rate thermocouple in the contact rod. The heat input for the sample is also 
ided hy the DSC The cell is calibrated with a sample of known thermal 
luctivity. The output of the DSC operating in the normal mode is used as 
ecorder zero, With the sample in position, the DSC cell is brought to the 
ed measurement temperature, 7, and when the output to the recorder is 
ly with time, the temperature difference across the sample, ΔΊΣ and the 
‘der displacement, h,, are determined. The sample is then replaced by a 
τ ας οὐκ. a stindard glass, and the corresponding quantities, AT, and h, are 
‘mined, The thermal conductivity of the sample. ἧς, is then calculated from: 


¢ Ris the recorder sensitivity, [the length and d the diameter of the sample 
icript 4} and ealibrant (subscript c) as specified, 
sample length of 10 fo 26 nm is recommended and a precision of better 
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Figure 4.10 Schematic diagram of thermal conductivity cell [10]. (With the permission 
of Thermochimica Acta, Elsevier, Amsterdam.) 


than 3% is claimed. Sircar and Wells [11] have used the above method in 
studying elastomer vulcanizates, and Hillstrom [12], a very similar system. for 
determining the thermal conductivity of explosives and propellants. 
Hakvoort and van Reijen [13] using a DuPont 910 DSC, replaced the 
temperature sensor on the top of the sample by a dise of pure gallium or indium 
metal (Fig. 4.11(a)). The temperature at the top of the sample during the 
melting of the metal is then fixed at the melting point, 7). of the metal. The 
y-output scale calibration factor © (in Wom ἢ. for the instrument is first 
determined by melting a weighed sample of the metal placed directly on the 
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(a) 


170 
160 


(b) 


150 


Figure 4.11 Measurement of thermal conductivity using DSC [13]. (a) sample 
geometry. (b) DSC plot, a: 7) against time—-d7;/dt, b: y against time-»dy/dr; 
(dT /d(dy/dt)=dT/dy= fi. (With the permission of Thermochimica Acta, Flsevier, 
Amsterdam.) 


sensor, without a sample pan, and measuring the DSC signal, y, with time. The 
heat of melting, Q,,, is given by 


On =f QO dt=C f dt. 


so that 


Qenivor = C¥ 


Two identical cylinders of sample (1-3 mm high, 2-4 mm diameter) are then 
‘laced on the sample and reference positions of the DSC (with silicone grease if 
recessary). A disc of the metal to be melted (10-40 mg and diameter similar to 
hat of the sample) is placed on top’ of the cylinder on the sample side and the 
ensor temperature, 7,, and the DSC signal, y are recorded with time as the cell 
s heated at constant rate. During the melting of the metal, the top of the sample 
ylinder remains at constant temperature, 7;,, while the lower side temperature, 
ἐὰν Increases at constant rate, The DSC plot obtained is shown in Fig. 4.1 1(h). 
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Two quantities are determined: the slope of the plot of 7, against time. 
a=d7./dt, and the slope of the DSC curve. bh=dy/dt. Combining these 
quantities, a/h=d7Ti/dy=f. The thermal conduetivity, 2, is then calculated 
from [4] 

2-hC/AR 


where 4 is the cross-sectional area of the sample cylinder and fits height. 

Boddington er al. [14] have carried out a more elaborate modification of a 
Stanton Redcroft DTA model 673 for measurement of the thermal diflusivity of 
pyrotechnic compositions. 


4.10 Determination of phase diagrams 


Tt is not always easy to detect accurately by eye the onset and completion of 
melting of binary or more complicated systems. The records obtainable from 
DSC and DTA runs, under carefully defined sample conditions. and slow 
heating and cooling rates so that equilibrium is approached. provide a more 
accurate way of establishing the phase diagram for the melting and solubility 
behaviour of the system under consideration. 

Ina binary system of immiscih Me solids and their prmneiete miscibly 
with the well-known phase diacram shown in Fig 4.12 πὸ 
acid/naphthalene system [15.Ψ [ΟἹ or the trphenyimethane crans stilbene 
system [17]), the melting curves of the pure components A and B (with melting 
points 7,,, and 7), ) are readily obtained, using DSC or DTA, and should be 
sharply defined. The melting behaviour of mixtures of A and B will depend upon 
the histories of the mixtures and, in this discussion, i will be assumed that the 
mixtures have been prepared in completely molten form initially Qvithout any 
volatilization or oxidation [18 ]). 


met Vanes eee 
B Xp mole fraction A 


Figure 4.12 Phase diagram fora binary system of immiscible solids and their completely 
miscible liquids. 
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A DSC or DTA record of the slaw cooling of such a molten mixture of 


ΠΤ Lobe ALE could 
show no deviation from the baseline until temperature /, when solid B begins to 
orystallize in. an exothermic process. This exotherm is not sharp but tails off as 
trystallization becomes complete (Fig. 4.13). The area under this exotherm will 
lepend upon the amount of B present in the sample of the mixture. As the 
emperature of the mixture falls further, the eutectic temperature, 7, . is reached 
ind solid A crystallizes in a sharp exotherm (Fig. 4.13). If the composition had 
seen chosen to be that corresponding to point F, the DSC or DTA record would 
yave shown only a single sharp exotherm as the eutectic composition solidified. 
At low concentrations of A, and similarly at low concentrations of B, the 
ormation of the cutectic may be difficult to detect and hence to distinguish from 
he formation of solid solutions [18]. 
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“igure 4.13 DSC or DTA record of the slow cooling of a molten mixture of composition 
, intially at temperature 1. in Pig. 4,12. 


Ifthe mixture, whose cooling curve is shown in Fig. 4.13, were to be reheated 
lowly, the DSC or DTA record should ideally be the endothermic mirror image 
that shown. The use of the heating rather than the cooling curve avoids 
sroblems of supercooling. 

Ideally then, the DSC or DTA traces should be readily relatable to the phase 
liapram. as ilustrated schematically in Fig. 4.14. The area under the peak for 
he cutectic melting tsa simple funetion of concentration (Fig. 4.15) and such a 
urve may be used to determine the composition ofan unknown mixture, c.g. an 
Hoy [19]. 

Pope and Judd [18] give an example of a more complicated phase diagram, 
‘ith incongruently and congruently melting compounds and solid solution, 
utectic and liquidus reactions, while Fysel [20] discusses the even more 
omplex Na,SO,/K,SO, system. A variety of systems, including ternary 
ystems and phase studics under high pressure, is discussed by Gutt and 
Majumdar {21}. and the theoretical background is covered by Sestak [22]. The 
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Figure 4.14 Schematic comparison of DSC or DTA traces with the phase diagram 
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Figure 4.15 Relationship between enthalpy of melting and composition 


phase diagrams for the ternary Ge Sb Bi system [23] and the complicated 
Ag, Te-Ag,SSe system, with 17 phase regions [24 |. serve as illustrations of how 
successful such studies can be. 

Combination of DTA with thermomicroscopy (chapter 5 and ref. 25) can 
provide additional information over a more limited temperature range. 


4.11 General applications of DTA and DSC 


The results obtained using DTA and DSC are qualitatively so similar that their 
applications will not be treated separately. It should be noted that DTA can be 
used to higher temperatures than DSC (max 725°C) but that more reliable 
quantitative information is obtained from DSC. 

Once again, DTA and DSC curves may be used solely for ‘fingerprint’ 
comparison with sets of reference curves. It is usually possible, however, to 
extract a great deal more information from the curves, such as the temperatures 
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denthalpy changes for the thermal events occurring. Asan example, the DSC 
rve for CuSO, . SHO ᾿ς shown in Fig. 4,16 for comparison with the TG 
pve given in Pi 3.9. Note the inerease in ΔΗ of dehydration per mole of 11,0 
noved 
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Figure 4.16 DSC curve for CuSO, . 541,0 


5. mentioned earlier, the shape of the melting endotherm can be used to 
nate the purity of the sample. The procedure is discussed in detail in 
Mer 14, 

election of solid: solid phase transitions, and the measurement of AH for 
c transitions, is readily done by DSC or DTA. A low-temperature 
chment permits the range ofsamples which can be examined to he extended, 
rown by the DSC curve for carbon tetrachloride in Vig. 4.17. Satisfactory 
‘ation as low as — 175°C has heen achieved. Condensation of atmospheric 


enco 


Monoclinic) hombohedral cubic liquid 


OSC resceonse 


T(°C)———> 


Figure 4.17 DSC curve for (Ὁ Ἴ: 
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moisture can cause problems and has to be reduced. Redfern [26] has reviewed 
some of the cooling systems available and has discussed some low-temperature 
applications. 

Application of thermal analysis to the study of polymers has been most 
rewarding. A DSC curve for a typical organic polymer is shown in Fig. 4.18. 
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Figure 4.18 DSC curve for a typical organic polymer 


Most solid polymers are formed by rapid cooling to low temperatures 
(quenching) from the melt and are thus initially in the glassy state. The 
transition from a glass to a rubber, the glass trensition, is an example of a 
second-order phase transition. Such transitions are accompanied by a change in 
heat capacity, but no change in enthalpy (A/7=(). The transition thus appears 
on the DSC curve as a discontinuity in the baseline at the glass-transition 
temperature, 7,. As the temperature is slowly increased the polymer may 
recrystallize giving the exotherm shown, before melting occurs. At higher 
temperatures the polymer may decompose (degrade) or oxidize depending 
upon the surrounding atmosphere. 

. The degradation or oxidation of polymers can be studied using DSC in the 
isothermal mode. Figure 4.19 shows the effect of stabilizers on the oxidation of 
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Figure 4.19 Effect of stabilizers on the oxidation of polythene, using isothermal DSC. 
(Cassel, B. (1975) Ind. Res. Aug, 51 56, with permission.) 
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polyethylene at 200°C, The crystallization of polymers can also be studied using 
isothermal DSC, Several approaches are possible, The molten polymer may be 
cooled rapidly to a temperature in the crystallization range, and the 
crystallization from the liquid observed, or the temperature may be raised from 
ambient so that crystallization from the rubber-like state is observed. 

Re-use of plastic waste is obviously most desirable, but is hampered by the 
problems of identification, sorting and collection. DSC may aid in identification 
of the constituents of the scrap, as illustrated in Fig. 4.20. 
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Figure 4.20 DSC of a sample of plastic waste (LDPE τεὶονν density polyethylene. 
HDPE =high density polyethylene; PP=polypropylene; PTFE = poly(tetraluoro- 
ethylene). ‘Teflon’). 


Anexample is given in Fig. 4.21 of the use of DSC in testing for completeness 
of curing of epoxy resins, The second scan shows no residual exotherm but 
indicates the glass-transition temperature for the cured resin. 

The behaviour of liquid crystals is of practical importance because of the 
optical phenomena which occur at the transition points, and their proposed role . 
in biological processes, Instead of going directly from solid to liquid at a sharp 
melting point, liquid crystals form several mesophases between solid and liquid 
states. The compounds which display this behaviour generally have large 
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Ligure 4.21 Curing of an epoxy resin 
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asymmetric molecules with well-separated polar and non-polar regions, In the 
smectic mesophase molecules are aligned and tend to form similarly aligned 
layers; while in the cholesteric mesophase the orientation of the molecular axes 
shifts in a regular way from layer to layer, and in the nematic mesophase 
molecules are aligned, but layers are not formed (Fig. 4.22). Changes from one 
mesophase to another can be detected using DSC, and the enthalpies of 
transition evaluated, as shown for cholesteryl myristate in Fig. 4.23. 
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Figure 4.22 Liquid crystal phases. (Based on Moore. Walter J. (1972) Physical 
Chemistry, 4th edn., with the permission of Prentice-Hall, Inc., Englewood Cliffs, New 
Jersey.) 


DSC can also be used to study the liquid-»vapour and solid vapour 
transitions [27, 28] and to measure the enthalpies of vaporization and of 
sublimation. In open sample pans these changes are spread out over too wide a 
temperature interval for an accurate baseline to be determined. The rate of 
escape of vapour is reduced by partial sealing of the pan. This sharpens up the 
trace as shown in Fig. 4.24, 

Structural changes in fats and waxes also show up clearly, making the 
technique (with low temperature attachments) ideal for food chemistry 
(Fig. 4.25). 
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Figure 4.23 DSC curve of a liquid crystal (chloresteryl myristate) (Rrennan, W. P. and 
Gray, ALP. (1974) Perkin Elmer Thermal Analysis Application Study, No, 13, April. with 
yermission.) 
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igure 4.24 DSC curves for water [27, 28]. (From Perkin-Elmer Thermal Analysis 
ewsletter No. 7, 1967, with permission.) 
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Kigure 4.25 Butter and margarine you can tell the difference! 
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Chapter 5 


Thermoptometry 


5.1 Introduction 


hermoptometry is a group of techniques in which optical properties of a 
imple are measured as a function of temperature [1]. Under this heading, the 
tajor technique is thermomicroscopy, i.e. direct observation of the sample. 
Mher techniques involve measurement of total light reflected or transmitted 
hermophotometry). light of specific wavelength(s) (thermospectrometry), 
‘fractive index (thermorefractometry), or emitted light (thermoluminescence). 


5.2 Thermomicroscopy 


is always interesting that the most obvious property of a material, icc. its 
ward appearance, is so often not examined or monitored during the heating 
ocess and yet can provide so much reliable first-hand evidence of the 
ocesses occurring [2]. Other thermal analysis techniques give results which 
ive to be associated by inference with thermal events and some processes such 
sintering, decrepitation and creeping and foaming of melts are only really 
tectable hy direct observation. Controlled temperature stages for micro- 
ρος are readily available and can cover the temperature range from — 180°C 
"10 2000°C, allowing for study of a wide range of problems, especially phase 
idies (chapter 4), Care has to be taken that suitable heat-resistant objective 
ises of adequate focal length are available. Standard photographic equipment 
video cameras can be used for recording observations. Allowance may have 
be made for heating by the illuminating source, or heat filters may he 
cessary. 

If visual examination can be combined (chapter 9) with a quantitative 
asurement such as DSC or DTA, even more meaningful results can be 
tained, The Mettler system [3,4] (Fig. 5.1) has a heat-flux DSC sensor built 
m thin-film thermopiles deposited on a special glass disc. The sample is 
ced in a sapphire crucible and a similar but empty crucible is used as a 
trence, Sample illumination is by transmitted (normal or polarized) light. 
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Figure 5.1 Instrumental set-up for simultancous hot stage microscopy [3.4]. (With the 
permission of Mettler Instrumente AG.) 


computer 
screen 


The temperature range is from room temperature to 300°C, Wiedemann and 
Rossler [5] have reported on the use of this system for a detailed study of the 
dehydration of gypsum. 


5.3 Thermophotometry 


In thermophotometry provision is made for measurement of the intensity of the 
light reflected or transmitted by the sample, If samples are viewed in transmitted 
light using crossed polars, the only light transmitted is that arising from 
rotation of the plane of polarization of the light. caused by changes in the 
crystalline structure of the sample. Melting or the formation of an isotropic 
structure results in complete extinction of the light. The information provided is 
thus complementary to DTA or DSC measurements. Arneriand Sauer [6] have 
shown how such a technique can be used in the study of polymers, where 
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overlapping transitions, such as melting and degradation or crystallization and 
melting. can lead to ambiguities of interpretation, 

In the Stanton Rederoft system [7]. the sample is contained between two 
Microscope cover glasses which rest on a heating block, within an atmosphere- 
controlled chamber. The block contains a sapphire window for light 
{ransmission, and a platinum resistance thermometer for measurement of the 
temperature and control of the heating or cooling programme. Provision for 
cooling with liquid nitrogen gives an overall femperature range of — 180°C to 
60°C. The intensity of the transmitted light is measured using a silicon 
photodetector. An example of the results obtained on heating KNO, is shown 
in Fig. 5,2, The phase transition from orthorhombic, phase H, to trigonal, 
phase 1 at 128°C and melting at 335°C show up clearly. The system has also 
been used [8] for determining the ignition temperatures of pyrotechnic 
compositions. An earlier version of the hot-stage microscope [9] had a 
femperature range up to 1000°C, but the sample had to be viewed in reflected 
light. 
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ure 5.2 Thermophotometry curves for potassium nitrate [12]. (a) In transmitted 
light. (b) in reflected light. (With the Permission of Thermochimica Acta, Elsevier, 
Amsterdam.) 


Sommer and Jochens [10] have described a micro-thermal analysis system 
based on the use ofa thermocouple as both specimen holder and heating source. 
The system is easily and cheaply constructed and is particularly useful in phase. 
studies of non-metallic samples (11). 

Haines and Skinner [12] have modified a DSC to carry out simultaneous 
measurements of the intensity of light reflected from the sample. Changes in the 
surface of the sample may not be accompanied by measurable enthalpy changes, 
but may show up clearly through changes in reflectance. A binocular 
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microscope was used so that one eyepiece could be used for photographing the 
actual appearance of the sample, while the photodetector was used on the other 
eyepiece. The sample was in either an open DSC pan, or else a mica window was 
crimped into the pan. Changes in the intensity of the reflected light are 
cumulative so that the curves recorded are related to DSC curves as integral to 
derivative (compare TG and DSC). The phase (transition at 128°C for KNO, 
did not show up as clearly as in transmitted light (Fig. 5.2). The transitions in 
liquid crystals were readily detectable. Dehydrations and decompositions were 
also accompanied by significant changes in reflected light intensity. The system 
also proved useful in the study of the browning and ignition of paper samples 


[12]. 


5.4 Thermoluminescence 


The technique involving measurement of light actually being emitted by the 
sample is known as thermoluminescence, TL [13] and has heen used. in 
particular, by Wendland{ in the study of coordination compounds | 14, 15]. ΤΊ 
arises from the annealing of defects in the solid sample as the temperature is 
raised. The defects are present in the sample on account of its previous thermal 
history, or through irradiation of the sample. Light emission from polymers 
[14] may result from the formation of hydroperoxide radicals, while vigorous 
redox reactions between oxidizing groups and reducing ligands in coordination 
compounds may also be accompanied by light emission [14, 15] 

Nuzzio [16] has described the modification of a DuPont 990 Thermal 
Analyzer for TL measurements including arrangements for computer data 
capture. The light-measuring system included fillers to absorb heat and to 
minimize the black-body radiation encountered at temperatures above 
~400°C. A reference light source was used for calibration. TL measurements 
could be made under isothermal or increasing temperature conditions. 
Arrangements for measurements under reduced pressures are also described. 

Manche and Carroll [17] modified a Perkin-Elmer DSC-1B for simultaneous 
TL and DSC, using a matched pair of optical fibres mounted directly above 
each calorimeter cup. The system was tested using mixtures of LiF and KNO,. 
LiF is a well-known and extensively-studied thermoluminescent material and 
KNO, was used as the temperature and calorimetric standard and does not 
normally exhibit thermoluminescence, Ground mixtures were pressed into discs 
and the discs were irradiated with X-rays. The sample discs, with a preheated 
disc of LiF as reference, were heated in the modified DSC. The results obtained 
are shown in Fig. 5.3, Areas under the glow curves and the DSC endotherms 
were measured and were shown to be linear functions of the mass fraction of the 
individual salts. 
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gure 5.3. Thermoluminescence (TL) and DSC curves for mixtures of LiF and KNO, 
7]. (With the permission of the American Chemical Society.) 


5.5 Electron microscopy 


lany of the kinetic studies of mechanisms of decomposition of solids 
hapter 13) invoke processes of formation and growth of product nuclei. These 
ieculations can often be confirmed or refuted using thermomicroscopy. 

Use of scanning electron microscopy, rather than optical microscopy. has the 
{vantages of greater magnification and much increased depth of ficld. There 
δι however, additional problems in the necessity for operating under vacuum 
id the effect of the electron beam on the sample. Many attempts have been 
ade to study reactions directly in the electron microscope. but none of these so 
r has had the temperature control required for thermal analysis. 
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rs Chapter 6 


Ὁ] Thermodilatometry (TD) 


6.1 Basic principles 


Measurement of the expansion of solids and liquids caused by various 
conditions. such as absorption or chemical reaction or other tin 
processes, at fixed temperatures is known as dilatometry. Whe 
emphasis is put on recording such dimensional changes as a fu 
temperature, during a controlled temperature programme, the tec 
labelled thermodilatometry [1] and is often included under the genera 
of thermomechanical analysis (TMA) (chapter 7). For solids, we 
tinguish between linear and volume expansion, and instruments forr 
the length of a solid sample as a function of temperature are the most 

Most solids expand on heating. Exceptions include vitreous silica a 
low temperatures. The change ina linear dimension, L, with Tis gi 

1 


L,=L,(14 ᾿χ 41) 


5 JT 
" * nm” . . 
ise ἢ where [., is the length at ΤΊ. and L, the len that 7, and α is the cor 
on 1 ' 2 2 
wi! linear expansion. Over a small temperature interval, AT, α is appr 
ict constant for many materials. so the dilation 
4 ee 
{ Ξ 
"Ὁ Ly~Ly=AL=L, ΔΤ 
or 
AL/L, =« AT 
ς The value of α is related to the structure and type of bonding in th 
ΧΙ general, strong bonding results in low values of a, so the order of 
fh ταῖς ες mes si 
ail coeflicients is usually; covalent and ionic materials < metals <m 
ar | bonded materials. In crystalline substances, α΄ may vary with 
} ᾿ 1 (anisotropic behaviour). 
4 Ἢ af Measurements of thermal expansion may be classed as either a 
ΠΣ ἶ difference measurements, as illustrated schematically in Fig. 6.1. T 
14 and seale are outside the furnace and the movements of the sample 
White 
Wy ia 


ἀρ | “Foie 
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Figure 6.1 Thermodilatometry. (a) Absolute measurement, (b) difference measuren 


reference are transmitted by identical push-rods. Measurements may be n 
in either the horizontal or vertical position. In horizontal instrum¢ 
arrangement has to be made for suitable spring-loading, so that any contrac 
can be detected. In vertical instruments, the weight of the push-rod 
measuring system has to be counterbalanced (o avoid compressior 
deformation of the sample. 


6.2 Practical details 


Measurement of sample movement obviously has to he more sophisticated 
the pointer and scale used in the schematic diagrams. The transducer (Fig 
is usually a linear variable differential transformer (LVDT), the princiy 
which is illustrated in Fig. 6.2. Movement of the core within the coils prot 
an output with sign dependent upon the direction of the movement, 
magnitude giving the movement of the sample. Amplification of the move 
is claimed to be as much as 20.000 x Calibration of the LVDT is carricc 
initially, using a micrometer screw pAuge incorporated in the measuring Sy 

Safety devices, to prevent damage caused by melting of the sample i 
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Figure 6.2. Principle of a linear variable differential transformer (LVDT) 


yparatns, are often included, These are activated by rapid sample contractions 
id power to the furnace is cut off immediately. 

Expansion of a standard material can be used as a temperature indicator. 
othermal measurements may be made to study creep or recrystallization 
-ocesses in the sample (see applications). Changes in volume of liquid samples, 
sntained in a suitable cylinder, can be converted into equivalent linear 
‘ovement by means of a close-fitting piston attached to a push-rod. A bleed 
alve on the cylinder ensures that no air is trapped in the liquid. 

Linseis [2] has described a laser dilatometer based on the principle of a 
fichelson interferometer. A schematic diagram is given in Fig. 6.3. The laser 
cam is split, after a 90° reflection, into two separate beams. The two beams are 
irected at reflectors attached to the ends of the sample and reference materials. 
fter total reflection, the beams are recombined and the interference patterns 
re examined using photodiode detectors and converted to measurements of 
xpansion or shrinkage. These dilation measurements are accurate to one- 
uarter of the wavelength of the laser and involve no mechanical connections 


furnace 


reference 


Figure 6.3 Laser dilatometer [2]. (With the permission of John Wiley & Sons, 
Chichester.) 


with the sample, so very rapid changes of dimension can he examined. The 
temperature range over which such an optical system can be used is also much 
greater than fora mechanical system and it is suggested [2] that it will prove 
useful in examining the expansion coeflicients of superconducting materials at 
liquid helium temperatures. 

Karmazsin et al. [3] have suggested use of an optoelectronic transducer in 
place ofan LVDT. This has particular advantages for simultaneous thermodila- 
tometric, thermoconductimetric and thermomagnetometric measurements as 
LVDTs are sensitive to electrical perturbations. A stabilized light source 
(Fig. 6.4) iuminates two photosensitive cells through two rectangular holes of 


screen 


photosensitive 
cells 


Figure 6.4 Optoelectronic transducer for thermodilatometry [3]. (With the permission 
of John Wiley & Sons, Chichester.) 


ascreen, which is attached to the push-rods of the dilatometer. Movement of the 
screen alters the relative amount of light received by the two cells. Dimensional 
variations down to 10 ? ym could be measured and speed of response was 
about ten times that ofan LVDT. Examples of simultaneous measurements are 


given [3]. 


6.3 Interpretation of results 


Values of the coefficient of expansion. 7, may be determined from the slope of a 
curve of L against 7, since 


AL/AT=AL/AT=L, % 
Alternatively, the time derivative of the dilation may be recorded, using 
electronic differentiation, and 

(di /dt)(dt/d7T)— Lb, % 
so 


dL/dt=L, mh where ᾧ —d7/dt —heating rate. 
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-engincering coefficient of expansion, α΄. is defined by 
of = (Ly Lyq)/(La(T~ 20”) 


sre 1, and Ly, are the lengths of the sample at 7°C and 20°C, respectively. 
Isnally it is the discontinuities in the curve of Lagainst Τ᾽ rather than actual 
τὸν of y, which are of main interest in thermodilatometry, as these indicate 
urrence of thermal events in the sample. These discontinuities are discussed 
he applications section below. 


6.4 Applications of thermodilatometry 


mentioned earlier, thermodilatometry (TD) and thermomechanical analysis 
MA) (chapter 7) are closely related and they have been applied in similar 
ds (sce also chapter 7). A very practical application of TD is in the selection 
suitable materials for use as brake-linings [4] Obviously a low coefficient of 
yansion, coupled with other qualities such as wear-resistance and a suitable 
‘ficient of friction, is desirable (Fig. 6.5). 
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mre 6.5. Thermal expansion of brake-lining materials [4]. (With the permission of 
ernational Scientific Communications, Inc.) 


Wendlandt [5] has used TD to study the change in structure from octahedral 
tetrahedral in the cobalt(I1) pyridine coordination compound, Co(py),Cl, 
ig. 0.6). 
Another application of TD has been that of Gallagher [6] to the study of the 
netics of sintering of powdered Chromindur alloys (Fe Cr Co), used in 
aking high energy magnets. Arrhenius-type plots (chapter 13) of In(1//o) 
ainst 1/T showed two regions with different activation energies. The values of 
eactivation energies corresponded well with the diffusion energies in the low- 
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Figure 6.6 Thermodilatometry curve for ΕἸ ΡΥ) ΑἹ, [4]. (With the permission of 
Analytica Chimica Acta, Elsevier, Amsterdam.) 


temperature face-centred cubic and the high-temperature (> 1250 ( body 
centred cubic phases. 

TD has also been used [7, 8] to study the sintering behaviour of kaohins and 
kaolinitic clays. On heating. kaolin loses structural water at about 550 © and 
forms the metakaolin structure. Above 960°C this converts to a spinel structure 
and then above 1100°C to mullite. The mullite formed sinters at higher 
temperatures. These reactions are accompanied by the shrinkages shown in 
Fig. 6.7. The TD curves of kaolins and kaolinitic clays may be used for their 
classification [7]. 

Some other recent applications of TD have been to study the effect of 
additives on the shrinkage of cellular concretes [9]: the defects in non- 
stoichiometric oxides [10]; the firing of ceramics and the properties of alloys 
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Figure 6.7 Thermodilatometry (TD) curve on kaolin [7.51 (With the permission of 
Fisher Scientific Co.) 
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Chapter 7 


Thermomechanical analysis (TMA) 


7.1 Basic principles [1] 


Thermodilatometry is carried out by measuring expansion and contraction 
under negligible loads. Further information of interest may be obtained by 
measuring the penetration, 1.6. the expansion or contraction of a sample while 
under compression, as a function of temperature. Alternatively, the extension, 
ie. the expansion or contraction of a sample under tension, may be measured as 
a function of temperature. These techniques, plus flexure and torsional 
measurements, are classified as thermomechanical analysis (TMA) and are 
obviously of great practical importance in materials testing (section 7.3) 


7.2 Practical details 


The apparatus used is based on that for thermodilatometry (chapter 6). The 
principles of penetration, extension, flexure and torsional measurements are 
illustrated in Fig. 7.1. Two main types of experiment are possible: (1) 
measurement of dilation with temperature ata fixed load, or (ii) measurement of 
dilation with load at a fixed temperature. 

The Mettler TMA40 instrument is illustrated in Fig. 7.2. The overall 
measurement range is +S mm about the zere position set by the height 
adjustment ring. The probe force is microprocessor controlled with the help ofa 
linear motor and ranges from 0 to 0.5 N. A dynamic force component can also 
be superimposed on the probe force during measurement (chapter 8) 

The length is calibrated by means of gauge blocks of known size. The 
maximum sensitivity is ~4 nm. The temperature is calibrated using # special 
sandwich sample built up, as shown in Fig. 7.3, from discs of pure metals, 
separated by alumina discs. As the metals melt, the probe shows step-wise 
changes in length. The calibration sample can obviously only be used once. The 
instrument may be used from — 100°C to 1000°C in atmospheres determined by 
the purge gas. The prebe force is calibrated with a known mass. 


Jnermomechanical analysis (LATA) 
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gure 7.1 Thermomechanical analysis (TMA); (a) penetration, (bh) extension, (c) 
xure, and (d) torsional measurements. 


Several different designs are available for the probe e.g. flat, pointed or 
uunded ends, 


7.3 Applications of TMA 


hese techniques have been developed more recently than the others described 
» far. but applications have grown rapidly. TMA and TD (chapter 6) 
leasurements are often carried out on the same sample with the same 
pparatus. Examples of both expansion and penetration measurements on 
coprene rubber [2] are given in Fig. 7.4. The coefficient of linear thermal 
cpansion, αν can be determined from the slope of the expansion curve (see 
vapter 6), Both the glass-transition (7,) and melting (7,) show up clearly. 
nother example [2] 1s the use of penetration measurements for testing 
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Cross section of the TMA40 thermomechanical analyzer. 

(a) measuringsensor (6) linear variable differential transformer 
(Ὁ) sample support () linear motor coil 

(c) purge gas inlet (f) linear motor stator 


Measuring principle of the TMA: 
4. metal furnace assembly with electric heating coil 


2. sample (coated sheet metal) 
3. sample holder {adjustable in height according to 
4. 


sample thickness) 
measuring sensor with ball point (weight is compensated 


through weighing cell built in below, 
position is measured with linear transformer) 


Figure 7.2 Mettler TA 3000 Construction of the TMA4O, (With the permission of 
Mettler Instrumente AG.) 
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gure 7.1 Temperature calibration for TMA, (With the permission of Mettler 
strumente AG) 
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mre 7.4 TMA and TD measurements on neoprene [2]. (With the permission of 
ernational Scientific Communications Inc.) TS 


lycthylene coated paper (Fig. 7.5). The structural transition and the melting 
ow up clearly even though the coating was very thin (<0.3 mm). Similar 
sasurements have been carried out on paint films [3] and can coatings. 
istering of coatings may be revealed [4]. 

The behaviour of polymer films and fibres on heating under load is of great 
actical importance. Figure 7.6 shows the behaviour of a polyimide film [7]. 
ie film starts to shrink near 400°C. A similar measurement [4] on a drawn 
lyester fibre is shown in Fig. 7.7. Daniels [1] has given examples of torsional 
‘asurements, 

Knowledge of the thermomechanical properties, especially the glass- 
insition temperatures, of polymers used as denture bases [5], is important as 
tre are considerable variations in temperature in the mouth, as well as during 
: manufacture of dentures, Huggett et al. [5] have used TMA measurements 
monitor the effect of various modifications of the polymerization of methyl 
thacrylate on the 7, of the resultant material. The factors expected [6] to 
‘rease the 7) of a polymer are: 
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Figure 7.5 TMA measurements on polyethylene-coated paper [2]. (With the permission 
of International Scientific Communications, Inc.) 
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Figure 7.6 TMA measurements on polyimide film under load [7]. (With the permission 
of DuPont Analytical Instruments.) 
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Figure 7.7 TMA and DTMA measurements on a drawn polyester fibre heated at 
10 K min, ἢ in air, under load [4]. (With the permission of Stanton Redcroft Lid.) 
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) The presence of groups in the backbone which increase the energy required 
for rotation, e.g. unsaturation and/or long side chains 

) Secondary bending between chains c.g. hydrogen bonding, 

) Crosslinking between chains. 

) Molecular mass. 

ἡ Copolymerization 

) Plasticization. 
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Chapter 8 


Dynamic mechanical analysis (DMA) 
or dynamic thermomechanometry 


8.1 Basic principles [1, 2] 


One of the techniques of thermomechanical analysis, described earlier in 
chapter 7, involves measurement of the change of a dimension of the sample as 
the load on the sample is increased, while the temperature is held constant. This 
is effectively the determination of a stress strain curve for the sample. The load 
is related to the applied stress (force/area) and the change in dimension is related 
to the strain (elongation/original length). 

Different types of materials give different types of stress-strain curves and. of 
course, the nature of the stress strain curve fora given material will change as 
the temperature at which the measurements are made is changed. 


8.1.1 Stress—strain curves 


The two extremes of behaviour are clastic solids, which obey Hooke’s law, tec.. 
that the applied stress is proportional to resultant strain (but is independent of 
the rate of strain), and liquids, which obey Newton's law. ie. that the applied 
stress is proportional to rate of strain (not the strain itself). Both laws only hold 
for small values of strain or rate of strain 


8.1.2 Tensile stress 


For a sample under tension the tensile stress, σ =force/area = F/VoZq and the 
tensile strain, s=elongation/original length το ἀντ. Hooke's law is ame or 
o= Fr where E= Young's modulus (units: force/area). Foran isotropic body 
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_ properties independent of direction) the changes in length and in width are 
ited by Poisson's ratio 


ἀνὰ (x 
v7 ; 2 (vy, ranges from 0,2 to 0.5) 
p Yo) \dx. - 


Ι.3 Shear stress 


irasample subjected to a shearing stress the shear stress, 7, = F/xy=F/A,and 
Σ shear strain, e,-dx/z=tan 0, Again, Hooke's law is σιαῖπ, Or σ, παι. 
rere Gis the shear modulus. G=o,/e,=F/A tan @ and for small 0), tan @~@so 
= F/A0. 


Fe 


1.4 Compression 


Then a pressure pis applied toa body of original volume, V, causing a volume 
range. AI’, the stress. ¢,=p. and the strain, δα τ -AV/Vo. Hooke’s law is 
"αν OF n= Bey, where B is the bulk modulus, Β΄ τεσ μίδῃτε —PVo/AV. 

The compressibility of the sample is defined as 1/B= —AV/pVo. Of the three 
ioduli, Κι Gand RB, only two are independent. They are related via Poisson's ty 


1tio 
E=38(1—2v,)=2(1 +¥,)G 


1.5. Newton's law for viscous liquids 


he stresses applied to liquids are usually shear stresses and Newton’s law is 


thear stress οἱ rate of shear strain) or 
de, 


στῆ -- 
ὑπ 


there 1 is the viscosity coefficient, 


1.6 Viscoelastic behaviour 


Aany polymers show behaviour intermediate between elastic solids and 
iscous liquids and such materials are classified as viscoelastic, Shear of 
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viscoelastic materials is described by ὃ combination of Hooke's and Newton's 
laws. Strains are additive, so when the total shear stress =4,, 


᾿ς ἀν Ve εἴν tr ἧς vie 
Feet τε σὲ (Hooke’s law) and (tits = 2 (Newton's law) 
cai dt " 
Hence 
ey inst τῷ (V da, 
dt G) dt 
so 


dey ἀπ νι ὄρνις 1\ do, δι 
at dt τ dt ἢ 


When the shear strain is constant, 


; \d 
dey Ὁ Ἐπ᾿ ) do, σ, 
dt σ; 


which rearranges to 


If σ, Ξε σο at t=0 


and o,=¢ at t= 


CAC ae σι 
eae oan 


: (<) G ( G ) 
ie. Inf{--)=——t oF σεέσῃ CxP --t 
σο " " 


This means that at a fixed strain, if the initial stress is a. this will relax 
exponentially with time. The relaxation time is /G. 


8.1.7 Stress-strain measurements on polymers 


Ifa viscoelastic sample is subjected to a tensile force, applied ata uniform rate, 
and the resultant elongation of the sample is measured, a curve of the form 
shown in Fig. 8.1(a) may be obtained. (The detailed shape will depend upon the 
rate of application of the stress.) In the region OL, Hooke's law is obeyed and 
the slope gives Young's modulus, F. Beyond the yield point, Y, viscous flow 
occurs until the maximum elongation is achieved at the break point, B. 
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Tigure 8.1 (a) Typical stress strain curve for a viscoelastic sample [1]. (Ὁ) Effect of 
emperature on stress strain curves [1]. (Based on Cowie. J. M. G. (1973) Polymers 
Chemistry and Physies of Modern Materials, Untertext, New York., Ch. 12. with 
sermission of Blackie and Son Lid., Glasgow.) 


The effects of temperature on the stress-strain curves are illustrated in 
Fig & 10). As the temperature increases (i) the modulus of elasticity, F, (i.e. 
dope of OL in Fig. 8.1(1)) decreases; (ii) the yield strength decreases; and (iii) 
he maximum elongation generally increases. 


RAR Periodic stresses 

In dynamic mechanical analysis (DMA) the sample is subjected to a 
sinusoidally varying stress of angular frequency, ὦ. Fora viscoclastic material, 
the resulting strain will also be sinusoidal, but will be out of phase with the 
applied stress owing fo energy dissipation as heat, or damping. in the sample 
(Fig. $.2(9)). δ is the phase angle between the stress and the strain 
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(c) Measurement of damping [1] 


Figure 8.2 (Based on Cowie. J. M.G. (1973) Polymers: Chemistry and Physics of Modern 
Materials, Intertext. New York.. Ch, 12, with permission of Blackie and Son Ftd. 
Glasgow.) 


8.1.9 The resonance frequency. Ὁ, 


Ifa sampleis vibrated over a range of frequencies and the amplitude of vibration 
is measured, the resonance frequeney is that which produces a maximum ina 
plot of amplitude against frequency. Young's modulus (of elasticity). F. 18 
related to the square of the resenanee frequency. ¥,- 


Exc Lpy?/B? 


where ¢ is a constant. 1. is the sample length between clamps. B is the sample 
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ckness and pis the sample density. Ifthe sample is allowed to oscillate freely. 
will do se at the natural resonance frequency with decreasing amplitude 
ving to damping, The form of the oscillations is shown in Fig. 8.2(b). The 
riod. Γι of the oscillations is constant (= t/1,) and the amplitude (measured 
sma minimum to the preceding maximum) decays exponentially. The rate of 
cay isi measure of how much damping there is. Dampingis usually expressed 
the logarithmic decrement per cycle, A= log) Δι 42) εἰοβιο (A,/Ay)=--- 
¢, Units used are decibels (dB) ie. 10x log (A,/A2). ee. If the amplitude 
ives, log (.4,/A2)— log 2= 0.30. Damping is then 3.0 dB 

Damping may be determined from the width of a curve of relative amplitude 
zainst frequency (Vig. 8.2(c)). h=(vy~¥,)/",. or by measuring the driving 
irce required to maintain a constant amplitude of vibration at the resonance 
equency (forced vibration). 

Damping depends on the physical state of the sample, e.g. at 7» Τὰν polymers 
issipate most of the energy supplied and damping is high, At T< T,. the 
jaterial stores the energy and damping is low. 


1.10 DMA apparatus 


*he DuPont DMA instrumentis illustrated in Fig. 8.3. The sample, which must 
cin film or fibre form, is clamped between the arms of the complex pendulum. 
“he arms-plus-sample system is set in motion by the electromechanical 


arm-locking pins LVDT adjustment screw 


electromagnetic 


driver a 


flexure pivot 


sample arm 


slide lock 
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ὸ control and sample 
Vernier adjustment knob thermocouples 


Figure 8.3 DuPont 98? Dynamic Mechanical Analyser. (With permission of DuPont 
Analytical Instruments.) 
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transducer. A linear variable differential transformer (LVDT, see Fig. 6.2) is 
used to monitor the frequency and the amplitude of vibration and feedback 
from the LVDT is used to control the electromechanical transducer. The preset 
oscillation amplitude can then be maintained and the driving force required to 
do so isa measure of the energy dissipation (damping) of the sample. The DMA 
outputs are thus plots of resonance frequency and of damping as functions of 
temperature. Young's modulus. Ε. for the sample is obtained from the 
resonance frequency, ν᾽, 


(4π2 v2? J, -K) (ἢ 


ὃ "] (5) rl Ἢ 


where J, is the moment of inertia of the arm system, K is the spring constant of 
the pivot, D is the clamping distance, 1” is the sample width, B is the sample 
thickness, and 1, is the sample length. 


Connop et al. [3] have made a careful study of the critical factors in using the 
DuPont 981 DMA, and have shown that to obtain reproducible results it is 
necessary to: (i) keep all sample dimensions within an aspect ratio of from 14 to 
20; (ii) mount the sample consistently with careful centering and tightening of 
the clamps; (iii) maintain a constant gas flow-rate: (iv) maintain a constant 
sample thermocouple position. 

Although samples should ideally be rigid enough to be tested as a sheet or 
rod, softer materials, such as thermoplastics or clastomers can be studied by 
using specially designed supports [4]. With thin films, buckling deformation 
can be a problem. Toth et al. [S| have developed a special sample clamp 
assembly which provides an arched sample cross-section. Miller [6] describes a 
technique for characterizing organic coatings by DMA by coating a multiple 
filament glass strand. The contribution of the substrate alone is subtracted to 
give the contribution of the coating A similar technique using ἃ thin steel strip 
as substrate can be used in examining the mechanical properties of paints. 

The polymer laboratories instrument (Fig. 8.4) normally operates with a bar 
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Figure 8.4 Polymer Laboratories Lid. Dynamic Mechanical Thermal Analyzer. (With 
the permission of Polymer Laboratories Ltd.) 


sample which is clamped at both ends. The central point of the bar is vibrated 
sinusoidally by the drive clamp on the end of a ceramic drive shaft. The stress 
experienced hy the sample is proportional to the current supplicd to the 
vibrator and the strain, which is proportional to the sample displacement, is 
monitored by the transducer, Softer materials are measured in the shear clamp 
assembly shown and liquid samples can be supported on films or absorbed into 


papers or braids, 


8.2 Applications of DMA 


Changes in Young's modulus indicate changes in rigidity and hence st rength of 
the sample. Damping measurements give practical information on glass 
transitions. changes in crystallinity, the occurrence of cross-linking, and also 
show up the features of polymer chains. The information obtained has heen 
used in various very practical areas such as studies of vibration dissipation, 
impact resistance and noise abatement. 

Typical DMA results [7] on two different samples of polyethylene are shown 
in Fig. 8.5(a) and (b), The damping curve for linear polyethylene (Fig. &.5(a)) 
shows peaks at 98°C and 68°C, ‘The lower temperature peak has been 
attributed to long chain ( CH) ), crankshaft relaxations in the amorphous 
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Figure 8.5 DMA curves of polyethylene samples [7]. (With the permission of DuPont 
Analytical Instruments.) 
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phase and the higher temperature peak to similar motion in the crystalline 
phase. The temperatures and relative sizes of the two peaks can be related to the 
degree of crystallinity of the sample. The damping curve for branched 
polyethylene (Fig.8.5(b)) has features at —112°C, --9Ὁ and 45°C. The 
112°C and 45°C peaks are explained as above, while the —9°C peak is 
attributed to CH, relaxations in the amorphous phase. 

In Fig. 8.6, the behaviour of styrene- butadiene-rubber (SBR) has been 
examined [7,8]. Various formulations of SBR.are used in tyre manufacture. 
Different styrene- butadiene ratios may be used, or different butadiene isomers, 
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Figure8.6 DMA curves on rubber blends [7, 8]. (With permission of DuPont Analytical 


Instruments.) 


or different additives, ¢.2. carbon black. A high cis-butadiene content lowers the 
glass transition temperature, T, (to as much as τ 110°C compared τὸ SOC) 
giving greater flexibility at low temperatures. The addition of carbon black 
(Fig, 8.6(c)) increases the modulus of elasticity. The 7, value is also slightly 
increased. The complex damping curve at low temperatures indicates 


length expansion 
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Figure 8.7 DMA curves ofa well and poorly-vulcanized polysulfide product. where the 
plastic flow begins at a lower temperature due to light crosslinking. Heating rate 
10 K min™!, temperature range 90 to 180°C, (Mettler Application Digest. with 


permission.) 
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polymer carbon black interactions and may lead to adverse properties, ¢.2 
heat build-up. 

Huson οἱ al. [4] have used the changes in height of DMA damping peaks as a 
measure of the extent of vulcanization in elastomer blends. Results agreed well 
with torque rheometer measurements. 

Gill [9] has given a general review of the application of DMA to polymer 
composites. 

Toth et al. [5] have used DMA to study the effect ofageing on paper. Samples 
were subjected to accelerated ageing at 100°C in air with 50% relative humidity 
and the DMA curves of original and aged samples were recorded. The increased 
brittleness of the aged samples showed up in the damping curves. 

Schenz et al. [10] have designed a sample mould for preparing frozen 
aqueous solutions for examination by DMA. The DMA technique was found to 
he very sensitive to second-order transitions. Frozen dilute solutions of sucrose 
showed two glass transitions at around — 32°C, suggesting that as the solutions 
are frozen, tiny localized inclusions are formed where the concentration differs 
from the bulk concentration. These inclusions have a 7, higher than that of the 
bulk solution, DSC was not sensitive enough to detect these details. 

TMA instruments (chapter 7) may be capable of imposing a periodic stress 
on the sample. An example of such ἃ measurement made on the Mettler TMA40 
(Fig. 7.2) is shown in Fig. 8.7. The extent of vulcanization of two polysulphide 
samples is compared. The load alternates with cycle time of 12 5. 

Kaiserberger [11] has compared the performance of TMA, DMA and DSC 
for determining the viscoelastic properties of polymers and concludes that 
DMA has superior sensitivity in detecting phase transitions of second and 
higher order. ™> 
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Chapter 9 


Combination of thermal analysis 


techniques 


9.1 Principles 


Results of the thermal analysis of a given sample, under a specified set of 
conditions, on a given instrument, are usually reasonably reproducible 
Agreement with the results on another portion of the same sample, obtained 
using the same technique on an instrument of a different make, may be much 
less satisfactory. It is even more difficult, to. compare results of parallel 
measurements from two or more independent thermal analysis techniques, ¢.8- 
TG and DTA. The advantages of being able to relate results of different 
techniques are obvious. For example TG cannot be used to detect melting. 
while melting and decomposition cannot be distinguished unambiguously using 
DTA. Asubstance which melts with accompanying decomposition must thus be 
studied using both TG and DTA (or DSC). There is usually a synergistic effect 
in that the total amount of information about the sample that is obtained is 
greater than the sum of the information abtained from the individual 
techniques. 

A truly simultaneous technique involves measurements of two or more of the 
properties in Table 1.1 on the same portion of the sample during ἢ single 
temperature programme. Each of the properties may be recorded continuously, 
or they may be sampled ina repetitive sequence 10 allow, for example, the data 
to be captured by computer (chapter 12). Simultancous measurements must 
thus be distinguished from parallel measurements, where different portions of 
the sample are examined using different instr uments, and concurrent measure- 
ments where different portions of the same sample, in different containers, are 
held within a single furnace and are subjected to a common temperature 
programme (Fig. 9.1). This last option is also used in instruments designed for 
carrying out one type of measurement on several samples. 

Simultaneous measurements obviously fake less time than sets of separate 
measurements, but gencrally the sensitivity of the individual techniques ts 
reduced on combinalion, because of compromises in instrumental design. In 
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Figure 9.1 Combinations of thermal analysis techniques. 


considering simultaneous methods, two combinations which probably fit the 
above definition are not considered further. These are the simultaneous 
recording of derivative curves of the basic signal, e.g. TG and DTG, and the 
combination of evolved gas analysis (EGA) with any other thermal analysis 
technique. EGA is discussed in detail in chapter 10 and its effect on t 
sensitivity of the basic technique is usually minimal. 

Leaders in the use of simultaneous methods, based on their instrument called 
the ‘Derivatograph’, have been J. and F. Paulik [1]. The Derivatograph 
(Fig. 9.2) is an example of one of the most successful combinations, i.e. that of 
TG with DTA to give TG-DTA. In the Derivatograph and some other 
instruments, the DTA reference material is not weighed. In the Stanton 
Rederoft TG-DTA system (Fig. 9.3) [2] the whole DTA configuration together 
with its thermocouple leads is incorporated into the balance suspension. A 
different approach is used in the Setaram simultaneous TG-DSC 111 (Fig. 9.4) 
[3]. There is no mechanical contact of the TG sample pans with the detectors of 
the DSC tubes. This is claimed not to affect the quantitative performance of the 
detectors. The limit of detection of the DSC is given as 30 pW in scanning mode. 
The symmetrical configuration of the system and the absence of contact 
problems are both factors in favour ofaccurate weighing (chapter 3). Demensky 
etal. [4] uscd a very similar system based ona Setaram DSC 111 and reported 
that the time constant of the calorimeter was ~ 20 s and that the DSC curve was 
thus distorted, They used a Fast Fourier transform method to obtain true DSC 
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Figure 9.2 1. and F, Paulik’s Derivatograph [1] for simultaneous TG-DTA. 1. Sample 
2. Reference 3. Furnace 4, Thermocouples. (With the permission of Flsevier, Amster 
dam.) 


exploded view of furnace area 


Figure 9.3 The sample holder assembly of the Stanton Redcroft STA-780 simultaneous 
TG-DTA system [2]. (With the permission of Stanton Rederoft Ltd.) 


curves, Gast ef al, [5] have described a magnetically coupled thermobalance 
adapted for simultancous DTA. 

Because of the sophistication of the principle involved, combination of 
power-compensated DSC (chapter 4) with another technique such as TG is not 
likely to be achieved. 


----.-- 
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gas 


igure 9.4 The Setaram simultaneous TG-DSC 111 [Ὁ] (With the permission of 
Setaram Lid) 


DTA may be combined with TMA by inserting thermocouples in the sample 
nd in the reference. DTA has been successfully combined with hot-stage 
ricroscopy (chapter 5) and even small-angle X-ray scattering [6]. Many of the 
‘ss-common properties, whose measurement is described in chapter 11, have 
ecn measured simultancously with more conventional properties. 

An example of a concurrent measurement system is the modification of the 
UPont 951 ΤῸ systems, described by Mikhail [7]. in which a sccond 
hermocouple was added to that already used for the TG, and sample and 
eference pans were placed in contact with the junctions of the two 
hermocouples, DTA measurements were then carried out concurrently with 
Ὁ on a separate portion of the sample, 
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Johnson and Ivansons [8] have described αὶ DuPont multiple-sample DSc 
which can run up to three samples concurrently, [is necessary that all positions 
for the samples and the reference should be thermally equivalent and that there 
should be no thermal communication (crosstalk) between the positions, Such a 
multiple-sample system, 10 addition to allowing for increased throughput of 
samples, also allows for inclusion of reference materials for calibration 
concurrently with the samples. 

While dealing with simultancous measurements, it is worth stressing the 
importance of complementary techniques which can aid in identifying the 
thermal events occurring in the sample, Complementary techniques include all 
forms of spectroscopy, X-ray and electron diffraction, optical and electron 
microscopy and are limited only by the ingenuity of the investigator 
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Chapter 10) 


Evolved gas analysis (EGA) 


10.1 Basic principles 


any samples, on heating. release gases or vapour through desorption or 
composition, This release is accompanied by thermal effects and, obviously, 
tss losses, which, themselves, can be detected by the appropriate thermal 
alysis technique, e.g. DTA or DSC and TG respectively. The thermal analysis 
‘hnique does not, however, identify the gas evolved and, for complex 
compositions, such information is essential. It has thus become fairly routine 
couple the basic techniques already described with a system for cither 
lecting the evolution of gas (or gases) from the sample (evolved gas detection, 
iD) or, more satisfactorily, detecting and identifying the gases evolved 
olved gas analysis, FGA) [1]. The apparatus for EGA will obviously be more 
plex than that required for EGD. 
» 


10.2 Evolved gas detection (EGD) 


iD has the advantage that measurements may be continuous and hence are 
dily related to thermal analysis curves. In most thermal analysis techniques, 
urge gas is used and this then becomes the carrier for sweeping evolved gases 
a detector (Fig. 10.1). The detector should be as close as possible to the 
nple to reduce condensation of vapours, secondary reactions in the gas 
ase, ind time lags between thermal analysis and RGD curves, The most 
nmonly used detectors [1] are those usually found in simpler gas 
omatographs, namely: (i) thermal conductivity detectors (TCD, or kathar- 
eters): (ii) gas-density detectors; and (iii) ionization detectors. In addition. 
has been made of infrared radiometers. 

τοῦ maximum sensitivity the measured property (e.g. thermal conductivity) 
he evolved gas should differ markedly from that of the carrier (Table 10.1). 
and He have very high. and argon very low, thermal conductivities, making 
m suitable carriers, but their influence on the thermal effects being examined, 
st be determined 
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Figure 10.1 Evolved gas detection (EGD). 


Table 10.1 Thermal conductivities of some gases 


Thermal conductivity, 


2(Wm 'K ') Differences, 100 A= 100A yg - Aree 

Gas T=300K T=373K Ref. N, Ref. Ar Ref, He 
Air 0.0262 0.0301 40.01 10.84 -12.47 
Ar 0.0178 0.0213 - 0.83 13.31 
He 0.1509 0.1669 + 12.48 +13.31 

H, 0.1869 0.2090 + 16,08 116.91 3,60 
N, 0.0261 0,0301 +0,83 12 48 
oO, 0.0266 0.0311 10,05 + 0.88 12.43 
co 0.0251 0.0305 ~0.10 40.73 12.88 
co, 0.0166 0.0212 0.95 -Ὁ 12 ~ 13.43 
NH, 0.0246 00,0327 0.15 Ἐ0.68 - 12.63 
N,O 0.0174 0.0212 ~ 0.87 —0.04 - 13.38 
NO 0.0260 0.01 +082 ~12.49 
NO, 0.0372 (328 K) 44.11 41,94 ~—11.37 
so, 0.0153 (est.) 11 ~-0.3 —13.6 
H,O(g) 0.0181 0.0241 ~ 0,80 40.03 ~ 13,2 


Sources: Weast, R. C. (ed.) (1984) Handbook of Chemistry and Physies, CRC Press, Boca Raton, 
64th Edn, p. E-2. Lange. N. A. (ed (1961) Handbook of Chemistry, M w-HHill, New York, ΠΗ 
Edn, p. 1544. Lodding, W. (ed.) (1967) Gas Effluent: Analysis, Fdward Arnold, London, p, 30, 
Thomas, L. ©. (1980) Fundamentals of Heat Transfer, Prentice-Hall, Englewood Cliffs, p. 665 


Flame ionization detectors are particularly sensitive to organic compounds, 
but do not respond to water vapour. 

Some separation of gas mixtures is possible by carrying out several EGD runs 
with suitable cold traps interposed between the sample and the detector. 
Alternatively, two similar detectors may be used, one on either side of the trap. 
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10.3 Evolved gas analysis (EGA) 
10.3.1 Mass spectrometry (MS) 


The most versatile and fastest means of repetitive gas analysis is mass 
spectrometry. so the most obvious solution to the problem of identifying the 
gases evolved from a thermal analysis instrument, is to replace the detector 
(Vig. 10.1) by a mass spectrometer [2, 3]. Such a replacement is, however, 
accompanied by several problems. Mass spectrometers are expensive items of 
equipment and dedicated use on a thermal analysis instrument is unlikely to be 
received enthusiastically by other potential users, The development of 
quadrupole mass spectrometers (QMS) at somewhat lower prices has eased this 
problem and Emmerich and Kaiserberger [4] have suggested that the capital 
already invested in the thermal analysis equipment warrants the further 
investment to get the maximum information per run. A second problem is that 
mass spectrometers require high vacuum for their operation, while most 
thermal analysis experiments are carried out at, or just above, atmospheric 
pressure. Some attempts have been made [ 2, 3] to carry out the thermal analysis 
experiments inside the vacuum system of the MS (Fig. 10.2(a)), Results 
obtained on heating a sample in vacuum are, however, often very difficult to 
relate to the behaviour of the sample under atmospheric pressure. It is thus 
more usual [4] to use some form of molecular leak inlet to introduce samples of 
the carrier gas and the accompanying evolved gases from the thermal analysis 
instrument into the MS (Pig. 10.2(b)). 

Problems that can arise include: (i) condensation of vapours in the sampling 
system attempts to reduce this by heating the system promote secondary 
reactions: (ii) the high concentration of carrier gas may swamp thgqsmaller 
responses of evolved gases - He, with its low mass number, is thus useful as a 
carrier, 

Av alternative to both of the above systems is to use a jet separator in the MS 
inlet (Fig. 10.2(c)). The jet separator allows the evolved gases to be separated 
from helium as carrier before introduction into the MS [5-7]. Dyszel [8] has 
described the interfacing of a TG system to an MS designed for atmospheric 
pressure chemical ionization (APCI). Gases entering the APCT source are 
ionized in a discharge corona. Only selected ions are then admitted to the MS’ 
through an electrostatically controlled inlet system. 

Complete mass spectra may be recorded repetitively, or selected mass 
numbers may be sampled using a suitably programmed system. ora single mass 
number may be monitored continuously. 

In interpreting the results of MS studies, allowance has to be made for the 
fragmentation patterns of the parent product molecules, Excessive fragmen- 
tation of the gas molecules occurs in some cheaper quadrupole mass 
spectrometers which do not have a variable ionization potential. Mixtures of 
CO, and CO, for example, then appear mainly as C' and O° fragments. 
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Figure 10.2 Coupling thermal analysis equipment to a mass spectrometer 


10.3.2 Gas chromatography (GC) 


Another approach to EGA is to precede the detector (Fig. 10.1) by a column of 
a suitable adsorbent, i.e. to pass the evolved gases from the thermal analysis 
through a gas chromatograph [9-11]. Sampling. of course, now becomes 
intermittent as time has to be allowed for the component of the gas mixture with 
the longest retention time to be eluted from the column, before the nex! sample 
is introduced. This is the main disadvantage of GC, since times required for 
adequate separations may be of the order of several minutes, This means that 
only a few samplings may be possible during ἡ rapid thermal event, and such 
events may even be missed completely. The advantage of GC is, however, that 
by suitable choice of column-packing material, most separations can be 
achieved and retention times, once determined, provide a simple means of 
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identification, The vast literature of gas solid chromatography ts available in 
deciding upon a suitable set of analysis conditions for a given sample and its 
expected products. Two useful references on the analysis of mixtures of the more 
commonly-cncountered gases are the books by Hachenburg [12] and 
Thompson [13]. The ultimate, of course, is to separate the evolved gases ona 
column and to confirm their identification through use of an MS as the GC 
detector. The gas chromatograph (which may he assembled from selected 
components) has to include a gas-sampling value (Fig. 10.3) operated by a 
programmable timer. 
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Figure 10.2 Gas-sampling valve (Taylor Servomex). (With the permission of Servomex 
Lad.) 


10.3.3 Infrared spectroscopy 


Morgan [14] has described the use of non-dispersive (i.e. fixed wavelength) 
infrared analyzers, coupled in series to a DTA apparatus, to measure the 
amounts of CO, and H,0 evolved by minerals and mineral mixtures. SO, is 
another common gaseous decomposition product and this was determined with 
a special electrolytic cell (see below) also in series. The infrared analyzers are 
shown schematically in Fig. 10.4. Two intermittent beams of infrared radiation 
of equal energy are passed simultaneously through the sample and reference 
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Figure 10.4 Non-dispersive infrared analyzer [14]. (With the permission of Wiley 
Heyden 1.14.) 


cells and into a detector. The detector consists of two symmetrical chambers, 
divided by a pressure-sensitive diaphragm and filled with a sample of the type of 
gas to be measured. The reference cell is filled with a non-absorbing gas. When 
the sample gas is passed through either of the analysis cells (depending upon the 
concentration range) it absorbs cnergy and this imbalance causes deflection of 
the detector diaphragm. The movement of the diaphragm is converted to an 
electrical signal. The evolution profiles are integrated and the apparatus is 
calibrated with standard gas mixtures. A relatively fast flow rate of cart ier gas is 
used and the delay in the EGA signal for CO, and SO, relative to the DTA 
signal was found to be negligible (~ § s) with a longer delay (10 5) for 11,0 
caused by temporary condensation between furnace and analysis cell 

Both inert and oxidizing atmospheres can be used so that solid gas 
interactions and combustion reactions can be studied, Filsinger and Bourrie 
[15] have used a very similar system consisting of a Netzsch thermal analyzer 
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with a Lira Model 202 IR analyzer for CO, to study the gases evolved in the 
decomposition of calcium oxalate. It was thought that the decomposition of 
CaCO, to CaCO, would be useful as.a calibration standard for CO, but the 
dispropertionation of the CO turned out to be significant at the temperatures 
used, even though thermodynamic equilibrium was not reached. 

The use of Fourier transform infrared spectrometry (FT-IR) for evolved gas 
analysis has been reviewed by Lephardt [16]. The potential use of FT-IR for the 
identification of products of thermal degradation had been suggested by 
Liebman et al. [17], who remarked on the specificity and the short 
measurement times required, as well as warning of the possibilities for 
interference from the rotational-vibrational fine structure of spectra of small 
molecules such as H,O and HCI. 

In considering the apparatus required [16] for FT-IR-EGA one is faced with 
some of the general problems discussed in chapter 9 on simultancous methods. 
Compromises have to be made between the ideal requirements for operation of 
individual instruments. c.g. high heating rates and low carrier-gas flow rates will 
give greater concentrations of degradation products in the analysis cell, but 
these conditions are not always suitable for accurate thermal analysis. 
Secondary reactions between gaseous products are also enhanced at low flow 
rates, Lephardt [16] deals with both the coupling of commercially available TA 
equipment to FT-IR instrumentation and the development of specialized 
thermal units for FT-IR-EGA. 

Lichman et al. [17] used a flow-through, heated ligkggpipe gas cell (30 cm 
long and 44 mm cross section). Collimated IR radiation was passed into an 
interferometer and the exit beam from the interferometer was focused on to the 
window of the cell. The beam emerging from the light-pipe was focused on to the 
IR detector. A predetermined number of interferograms (usually 50), were 
signal averaged and each set was stored for later recall and transformation, The 
resulting spectra were ratioed against a single-heam spectrum of the 
background, Spectra were measured at 4cm~! resolution. The time interval 
between spectra. for 50 scan sets, was 2 min, but this could be reduced by 
reducing the number of scans. This increased the noise, especially above 
2500 em ὅν, The spectra obtained [17] showed that FT-IR allows all volatile 
products fo be monitored simultaneously. 

Lephardt and Fenner [18] ensured that interferograms were accurately 
associated with sample temperatures. They also arranged for transfer of large 
arrays of spectra from the spectrometer to a secondary computing system, so 
that the dedicated computer could continue the data capture, while the larger 
computer did the data processing. Output was available as absorbance against 
temperature plots. The system was applicd to the pyrolysis of tobacco [18] 
Later developments | 19] included standard routines for smoothing (Savitsky 
Golay) and integration (trapezoidal rule) (see chapter 12) as well as more 
advanced processing of the EGA interferograms, using comparison with 
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interferograms of compounds of interest, to deconvolute overlapping profiles 
into their components. In addition to compound-specific profiles, reaction- 
specific profiles could be produced. 

Davidson and Mathys [20] describe the use of a ratio recording dispersive 
infrared spectrometer for FGA. Gas cells of from [0 cm to multipass cells of Fm 
path length and low volume (45 ml) were used and digitized spectra at a 
resolution of 5 cm~! were recorded in a scan time of | to 2 min and smoothed. 
Heating rates of 5K min τα were used to allow adequate definition of gas 
evolution profiles. Most effluents yield relatively simple spectra, and evolution 
profiles were constructed from measurements on selected absorption bands 
‘The temperature corresponding to each measured absorption was determined 
from the scan time and the heating rate, with allowance for times required for 
writing data to disk. The gas transfer lines and sample cell were at the 
temperature of the spectrometer (32°C) and water condensation problems were 
not experienced at the low levels used. The pyrolysate from some polymers, ¢.2. 
polyethylene, could cause contamination problems. 


10.3.4 Special-purpose detectors 


There is no reason, in principle, why any means of analyzing gases should not be 
coupled to a thermal analysis instrument, and there are reports of the use of 
absorbents and volumetric methods for the determination of total amounts of 
evolved gases. Detailed information on the evolution of each component with 
time is, however, most desirable. 

Water vapour is a very common product of thermal decompositions, as well 
as being a product of the reduction of metal oxides with hydrogen, and special 
detectors have heen developed to monitor evolution of water, in addition to the 
use of non-dispersive infrared analyzers described above [14]. Many of the 
common methods for determining the water content of gases are not entirely 
suitable when the water content is varying fairly rapidly, on account of slow 
and/or non-linear detector responses. For EGA the response of the detector 
should also be selective for water as other products may be evolved 
simultaneously. 

Warrington and Barnes [21] have shown that an electrolytic hygrometer is 
suitable for continuous water analyses. The hygrometer has two fine platinum 
wires wound closely. but not in contact, on ἃ PTFE former, The wires are 
coated with phosphoric acid and the whole element is enclosed in a flow- 
through glass tube which may be coupled to the outlet ofa TA apparatus. Atthe 
start, the acid is electrolyzed to dryness using a potential of 100 V hetween the 
two wires. and, in this state, there is then a negligible current between the two 
wires. The acid coating absorbs any water from the gas stream passing through 
the tube and this water is then clectrolyzed. The electrolysis current is 
proportional to the water concentration in the gas. Precautions necessary to 
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avoid spurious effects are described [21], and the system was tested on the 
dehydration of several hydrates and the decompositions of several carboxylic 
acids, 

Gallagher et al. [22, 23] have developed an EGA system based on a 
Panametrics Model 700 Moisture Analyser in which the moisture content is 
determined from the dew point of the gas stream and the flow rate. The 
relationship is non-linear and the computations are described [22], Allowance 
has to he made for background moisture in the gas stream and for degassing of 
the TA apparatus as the temperature rises, as well as loss of moisture by 
adsorption on cooler surfaces of the system. 

Morgan [14] has described an electrolytic detector for SO, analysis 
(Fig. 10.5) based on a fuel cell. SO, in the carrier gas diffuses through a 
semi-permeable membrane and is adsorbed ona sensing electrode, producing a 
current in the circuit which is proportional to the partial pressure of SO, in the 
carrier gas. Calibration using a standard gas mixture is required. 
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Figure 10.5 Electrolytic cell for SO, analysis [14]. (With the permission of Wiley~ 
Heyden Ltd.) 


Gauthier et al, [24] have discussed the use of solid electrolytes of K SO, for 
analysis of SO, and K,CO, for CO,. The latter has been used ina study of the 
decomposition of CaCO,. 

A review in Analytical Chemistry [25] points to some of the revolutionary 
new sensors under development, some of which will undoubtedly find 
application in EGA systems of the future. 

Chemical conversion agents have been suggested to simplify the final 
analysis. c.g. use of 1,0, to convert CO to CO, and of CaC, to convert H,0 to 
CH, for easier GC analysis. The kinetic aspects of these processes introduce 
undesirable uncertainties in EGA, 


Applications of EGI7 ana ΤΑῚ 
10.4 Applications of EGD and EGA 


The main use of EGD is to distinguish between phase transitions and 
endothermic decompositions. This is of particular use in studies on coordina- 
tion compounds, where loss of ligand is usually endothermic and many 
structural changes, including melting, are possible. 

In Fig. 10.6(a), EGD curves using a thermal conductivity detector with and 
without a preceding liquid air (90 K) trap, and recorded simultancously with 
the DSC trace, are shown for the decomposition of nickel oxalate dihydrate. 
These responses were interpreted as dehydration followed by ἃ complex 
endothermic decomposition with evolution of mainly CO, and traces of CO 
(not condensed at liquid air temperature). 

A closely-related study on the isothermal decomposition of nickel formate 
dihydrate in the DSC at 480 K with simultancous EGA by GC is illustrated in 
Fig. 10.6(b). There is concurrent evolution of CO,, CO and HO, The relative 
proportions would have to be determined by calibration e.g. by dehydrating a 
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Figure 10.6 FGA coupled with DSC used in the study of the decomposition of metal 
carboxylates. (a) EGD for nickel oxalate dihydrate, (b) EGA by GC for anhydrous 
nickel formate (isothermal DSC; 480 K) 
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known mass of hydrate or decomposing compounds known to yield only one 
product e.g, CO, only or CO only. 

Schmid and Felsche [26] used TG and DTA coupled with MS to study the 
thermal decompositions of hydrated copper(II) and cobalt(II) tartrates. Their 
results for the cobalt salt are illustrated in Fig. 10.7, The m/e traces of the 
strongest signals are shown and the initial endothermic process is clearly 
identified as dehydration only, while the higher temperature endothermic 
decomposition of the tartrate group gives rise to several fragments as shown. It 
is not possible to determine definitely whether these fragments are primary 
decomposition products or the result of secondary processes in the ion source of 
the MS, Detailed decomposition mechanisms were proposed [26]. 
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Figure 10.7 1G, DTA and MS curves for the thermal decomposition of hydrated cobalt 
tartrate, Co(C,H,O,) 2.5H,0, heated at 5K min”! in an argon atmosphere [26] 
(With the permission of Thermochimica Acta. Elsevier, Amsterdam.) 
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Figure 10.8 DTA FGA of CuSO, . 5H,0 using non-dispersive infrared analyzers [14] 
for analysis of H,O and SO,. (With the permission of Wiley Heyden Ltd) 


The DTA and FGA curves for CuSO, . SHO using non-dispersive infrared 
analyzers [14] for H,0 and SO, are shown in Pig, 10.8. Once again the 
increased information obtainable by adding FGA toa TA system is illustrated 

Most of the references to the instrumental systems described above also give 
some illustrations of their application. The possibilities are virtually unlimited 
and in many situations (e.g. the degradation of polymers) it is essential, for 
practical applications. to have information on the gaseous products, especially 
hazardous products such as HCLand HCN. 
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Chapter 11 


Less-common techniques 


Under this heading, referring to Table 1.1, come emanation thermal analysis 
(ETA), thermomagnetometry, thermoelectrometry and thermosonimetry 
These techniques are all gaining in popularity but usually require more 
specialized equipment which is often not commercially available as complete 
systems. 


11.1 Emanation thermal analysis (ETA) 
11.1.1 Introduction [1, 2] 


Emanation thermal analysis (ETA) involves the measurement of the release of 
inert (and usually radioactive) gas from an initially solid sample, on heating. 
The rate of release of gas is used as an indication of the changes taking place in 
the sample. Most of the solids to he studied do not naturally contain inert gas 
There are thus several techniques for incorporating inert gas. 


11.1.2 Sample preparation 


Techniques for incorporating inert gas ina solid sample may be divided broadly 
into two groups: (A) techniques for introducing the parent nuclide of the inert 
gas, αι». 
225 ΤΉ - ᾿ς 3 * ,220Rn, 

(B) techniques for introducing the inert gas itself, Labelling with parent nuclide 
gives a sample which is stable for longer as far as a source of inert gas ts 
concerned. Gas is being formed continuously and will not all be lost in a single 
run on the sample. Under these headings, the following methods have been 


used. 
Al Coprecipitation of parent isotopes. The disintegration process ensures a 
random distribution of the products. 
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2. Impreqnation of the sample with a solution of parent isotope is used when 
coprecipitation is not possible, Parent atoms are distributed only on the 
surface, Disintegration results in penctration of the daughter atoms into 
crystallites, but unless the crystallites are small the distribution is non- 
uniform, The distribution may be altered by annealing the sample. 

(2 Nuclear reactions may be used to produce the parent nuclide or the inert 
pas, c.g. 

AAG Ey IK rs 44Ca ee Ar 


1277 mn » 125] ν "Χο 


21 Jon bombardment is ἃ common method of implanting inert gas atoms in 
solids. The ions may be generated in vacuum using an electrical discharge or 
amicrowave plasma. The quantity of gas absorbed depends on the gas used, 
the energy of bombardment and the nature of the solid. 

32. Diffusion at high temperature and/or pressure is also sometimes possible. 
Ἀπ Καὶ is usually used. (Solids labelled with **Kr are called kryptonates. 
1,910.76 years for -emission.) 

31. Crystallization or sublimation of the sample in an inert gas atmosphere may 
be more effective in incorporating gas. 


[1.1.3 Measurement of gas release 


Details of the gas release are very dependent upon the way in which the sample 
was originally labelled, Release may involve either bulk or defect diffusion, 
There may also be recoil ejection, The emanating power, Εὶ is defined as 


E= Nye! Nroem 


where Nee is the release rate and Nie is the formation rate. Neorm May be 
measured by dissolution of the sample in acid or other solvent and measurement 
of the gas released. 

ETA is usually carried out in conjunction with other TA techniques, e.g. 
Netzch markets EVA/DTA/EGA apparatus. Carrier gas, at an accurately 
sontrolled flow rate is used to carry released gas to suitable counting chambers, 
Rn, an o-emitter, requires a scintillation counter, or an ionization chamber or 
zemi-conductor detector, while Geiger counters are used for Kr, Xe and Ar (f- 
emitters). For experiments lasting a long time compared with the half-life of the 
inert gas used, the decay of the measured gas should be taken into account. 


11.1.4 Applications of ETA 


In general, there is ood agreement between FTA results and those of DTA and 
DTG. One of the major uses of ETA is in the characterization of powders. 
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Emanating power is related to surface area and hence changes in grain size and 
the occurrence of sintering during heating may be detected. 

Figure 11.1 shows a series of isothermal PTA curves for NiO. samples. 
prepared from the carbonate. These curves, obtained ina nitrogen atmosphere, 
were used [3] to determine the kinetic parameters for sintering of the samples. 
The apparent activation energy of 276 kJ mol | was reduced to 155 ΚΙ mol 
when a similar series of runs were done in oxygen. Phase changes also show up 
as changes in emanating powcr, 6.8. the orthorhombic to rhombohedral phase 
transition of KNO, at 128°C (Table 4.2) is seen in Fig. 11:2. 

ETA curves for decompositions (including dehydration) are similar to those 
obtained using DTA or DTG. Solid gas reactions have also been studied, e.g. 
the oxidation of labelled metals or reduction of labelled oxides. FTA and EGA 
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Figure 11.1 Isothermal ETA curves for NiO powders in nitrogen [2. 3]. With the 
permission of Flsevier Scientific Publishers. Amsterdam.) 
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Figure 11.2 ETA traces for the heating (1) and cooling (2) of KNO, through the 
orthorhombic-rhombohedral phase transition at 128°C [2.4], (With the permission of 
Elsevier Scientific Publishers, Amsterdam.) 
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curves are in good agreement. Information on solid--liquid reactions, e.g. the 
hydration of cement, has been obtained. In both solid-gas and solid-liquid 
reactions the formation of a layer of solid product may hinder emanation. 
For solid solid reactions such as spinel formation [2, 5], 6.8... ZnO(s) + 
Fe,O,(s)>ZnFe,0,(s), with ZnO labelled with 22*Th, the emanation of the 
individual components is checked first. ETA traces for the reaction are shown in 
Fig. 11.3, The oxides interact in a series of stages. Product begins to form by , 
surface diffusion between 250 and 400°C as shown by the increased emanation 
(curve 3), The DTA and TD curves (curves 2 and 1 in Fig. 11.3) show no 
changes at these temperatures. Sharp changes in all of the curves in Fig. I1.3are 
abserved at 670 to 700°C, The exotherm on the DTA curve is small, but the 
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Figure 11.3 Spinel formation by the reaction of ZnO-Fe,O, mixtures. Curve | is the 
TD. curve 2 the DTA and curve 3 the ETA trace, all at 10 K min” '. The ETA trace 
during reheating of the reacted mixture is shown in curve 4 [2, 5]. (With the permission 
of Fisevier Scientific Publishers, Amsterdam.) 


increased emanation (curve 3) shows up clearly. This corresponds to interac- 
tion by volume diffusion. The reaction is complete by about 800°C. The dilation 
shown in curve 1 is suggested [2] to be caused by the formation of a very finely 
powdered product which sinters at higher temperatures. The ETA curve during 
the second heating of the reaction mixture (curve 4) shows that reaction was 
complete during the first heating, The reactivity of different preparations of 
Fe,O, has been compared using this technique [2]. 

The main drawback of ETA is that preparation and handling of samples 
requires all the usual radiochemical facilities and precautions. The amounts 
used in samples are so small that evolved gas, after dilution with carrier, does 
not form a hazard, Inert gases are not incorporated biologically and the decay 
products are stable. so hazards are reduced, 
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11.2 Thermomagnetometry (TM) 
11.2.1 Introduction 


When a sample is placed in a magnetic field, it may experience either attractive 
or repulsive forces. Attractive forces arise from three types of property of the 
sample: antiferromagnetism, paramagnetism and ferromagnetism (with ferro- 
magnetism causing the strongest interaction). The diamagnetic properties of a 
sample, arising from the orbital motion of electrons, and hence present in all 
samples, but sometimes swamped by the other magnetic properties, give rise to 
weak repulsive forces. Antiferromagnetism occurs in only a few transition-metal 
compounds. 

The volume magnetic susceptibility « of a sample is defined as «= M/F. 
where M is the magnetization and /1 is the magnetic field strength. The mass 
magnetic susceptibility 7 is then —x/p, where p is the density of the sample. 
Materials in which χ is positive are called paramagnetic, and those for which χ is 
negative are called diamagnetic. (An electronic microbalance is ideally suited 
for measurement of magnetic susceptibility [6].) 

Different types of magnetic behaviour can be distinguished by their 
temperature dependence. For paramagnetic samples, the magnetic susceptibi- 
lity decreases with absolute temperature according to the Curie Weiss law: 
y= C/(T—9), where Cis the Curie constant and @ is the Weiss constant. The 
susceptibility of diamagnetic materials does not change much with temperature, 
while ferromagnetic materials have high susceptibilities at temperatures below 
what is known as the Curie point. Above their Curie points, ferromagnetic 
materials hecome paramagnetic and these magnetic transitions are used for 
temperature calibration of TG systems (chapter 3). It should also be noted that 
Curie points can be determined by impedance measurements [7]. 


11.2.2 Apparatus 


Both the apparatus used, and the techniques and precautions required for 
thermomagnetometry are basically those required for thermogravimetry (TG) 
(chapter 3), with the addition of means of producing a strong magnetic field 
(> 10* gauss) around the sample. The apparent mass of the sample is then the 
sum of the actual sample mass and the magnetic force. Allowance obviously has 
to be made for the effect of the field on the sample container, which should thus 
have a low magnetic susceptibility ¢.g. quartz. The magnetic field may he 
applied periodically so that measurements of actual mass (TG) and apparent 
mass (TM) can be compared. Interference of the magnetic field with the balance 
mechanism is avoided by use of long suspension wires and suitable magnetic 
screening. 
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2.3 Applications of TM [12] 


ie major applications of TM have been to ferromagnetic samples where the 
ignetic effects are strongest. For example, 1M has been used [8] to examine 
> spinel phases formed during the thermal decomposition of siderite (FeCO). 
ie TM and TG curves shown in Fig. 11.4 for siderite samples heated in 
rogen, show the onset of decomposition at around 400°C, As the wustite 
ce) originally formed is oxidized by evolved CO,, magnetite (Fe,O,) is 
med, The apparent mass gainin the TM curve occurs as the magnetite nuclei 
ow and crystallize, When the temperature rises beyond the Curie point of 
ngnetite, the TM trace coincides with the TG trace in the absence of a 
agnetic field, In oxygen, the wustite is oxidized so rapidly to haematite that the 
rongly magnetic spinel phase never has a chance to form. 
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(a) TG of siderite (b) TM curves for siderite 


igure 11.4 TM and TG curves of siderite (FeCO ,) heated in N, and in O, [8] (With 
ie permission of Thermochimica Acta, Elsevier, Amsterdam.) 


The conditions for optimization of the formation of the ferrite NiFe,O, from 
ecomposition of the NiFe citrate have also been determined through use of 
M [9]. 

TM provides a means of studying reactions in sealed systems [5]. δῷ: 
orrosion, provided that there is at least one ferromagnetic reactant or product. 
‘igure 11S shows the apparatus used and the results obtained for the corrosion 
f carbon steel in an ammoniated EDTA solution. 

Another application is in the determination of pyrite in coal [11], usually in 
ombination with EGA, The residual ash from normal proximate analysis is 
eated in a reducing atmosphere (H,/N,) in the presence of a magnetic field 
Fig. 11.6), Prom the apparent changes in mass, as the Fe,O, is reduced to Fe 
nd the product is cooled, the pyrites content can be calculated [11]. 
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Figure 11.5 TM apparatus used, and results obtained 
of John Wiley & Sons. 
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11.3 Thermoelectrometry 


11.3.1 Introduction 


The main electrical properties of the sample which may be measured as a 
function of temperature are resistance (or conductance) and capacitance. A 
further technique involving measurement of the EMF generated when two 
dissimilar metal electrodes are in contact with the sample during a heating 
programme, is known as thermovoltaic detection [13]. Most thermoelectro- 
metry studies are carried out simultancously with other techniques, especially 
DTA. TG studies on the decomposition of solids have also been carried out in 
the presence of applied electrical fields. 


11.3.2 Apparatus 


Both DTA and DSC cells have been modified for simultaneous measurements 
of electrical resistivity [14]. The system used [14] with the DuPont DSC is 
shown in Fig. 11.7. The DSC cell base is insulated from the rest of the cell with a 
thin glass slide. A thin piece of asbestos was used in the reference position to 


glass capillary Pt electrode 
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asbestos pad nT oo 
 ———— 


Pt foil 
glass slide’ —— 
ref, platform sample platform 
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Figure 11.7 DuPont DSC cell modified for clectrical resistivity measurements [14] 
(With the permission of Thermochimica Acta, Flsevier, Amsterdam.) 


compensate for the heat capacity of the sample. Platinum foil was used to form 
an electrical connection from the bottom of the sample to the top of the 
reference. Glass-insulated platinum wire electrodes contact the tops of the 
sample and the reference. The current through the sample was recorded as a 
function of temperature. The dimensions of the sample were then used to 
calculate the resistivity. Sample surfaces were coated with colloidal graphite to 
overcome contact resistance, The sensitivity of the DSC was reduced 
considerably by the presence of the glass slide. 

Measurements of dielectric constant [ 15, 16] are based upon the measure- 
ment of the phase shift and attenuation of a sinusoidal AC signal after passage 
through the sample. The sample cell for concurrent dielectric measurement and 
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PA [15] is illustrated in Fig. 11.8. A coaxial-type two-terminal electrode 
afiguration is used. The inner electrode is a silver rod positioned 
nmetrically with respect to the outer electrode, which is a thin silver foil 
sssed against the walls of the cavity in the nickel block. The sample may he a 
uid, a powder or a machined solid. Approximately 500 mg of sample is 
juired. A separate smaller (50 mg) sample is used for concurrent DTA. The 
nple is connected to the inverting input of an operational amplifier, 
nfigured as a current-to-voltage converter with capacitive feedback. The 
‘asured phase-shift and attenuation of this network can be related [15] to the 
Jectric properties of the sample. 

The apparatus can be operated in two modes: either using a continuous range 
frequencies (50 Hz to | MEIz)at selected temperatures, or over a continuous 
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ure FE8 Cell for measurement of dielectric constant [15]. (With the permission of the 
terican Chemical Society) 
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range of temperatures al selected frequencies. Changes in the dielectric 
properties of the sample. brought about by phase transitions or chemical 
reactions are better resolved using the second mode, The instrument was 
calibrated with standard organic liquids. 

In the earlier apparatus of Chiu [16] based ona DuPont 941 TMA module, 
an automatic digital capacitance bridge with three fixed frequencies (120, 400 
and 1000 Hz), was used. Operation of the apparatus was checked using the 
transitions of KNO,. 


11.3.3 Applications of thermoclectrometry 


Electrical conductivity measurements are useful in detecting the appearance of 
liquid phases in the reactions between initially solid reactants. It is also possible 
to monitor the loss of chemisorbed OH groups in quantities too small to be 
detected by TG. Changes in electrical conductivity during heating have also 
been attributed to changes in concentrations of crystal defects. 

The dehydration of CuSO, SHO has been studied by measuring the 
electrical resistance of a pressed sample gimultancously with thermodilate- 
metry. The evolution of the five molecules of water in three steps is observed, as 
was found using TG (Fig, 3.9) and DSC (Fig. 4.16). 

Measurements of resistivity against temperature for polymers [14] show 
sharp drops in resistivity at about the glass-transition temperature. The 
formation of conductive carbon chains by cross-linking and cyclization within 
the polymer can be detected by the gradual decrease of resistivity with 
increasing temperature. Such processes are not shown up by TG or DSC. A 
curve for PVC is shown in Fig. 1 1.9, with the glass transition and the region of 
dehydrochlorination marked, Thermoelectrometry is also useful in studying the 
effect of carbon black additives on the properties of polymers [14]. 

The thermal decomposition of oil shales has been shown [ 17], using 
thermoelectrometry, to be a two-step process involving the breakdown ofan 
outer shell polar bridge (180 350°C) and cleavage of an inner naphthenic 
structure (350 500°C). 

AC electrical conductivity measurements (up to 350°C) on the Group I metal 
and ammonium perchlorates [17] showed that identical charge conduction 
mechanisms attributed fo movement of interstitial cations were present mn all 
these salts. 

Changes in the dielectric constant. rv, ofa material with temperature can arise 
from changes in molecular orientation brought about by phase transitions oF 
chemical reactions. 

Acurve of dielectric constant against temperature for CuSO, . SHO] 18] has 
the same form as the TG curve (Fig. 3.9). The curve of NaNO, (120-200 [9] 
shows ἃ peak αἱ ~ 165°C. of a form similar to that observed in the DTA curve, 
for the order disorder crystal transition. This peak is attributed (ὁ Debye 
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et 11.9 Thermoclectrometry curve for PVC [14]. (With the permission of 
sochimica Acta, Elsevier, Amsterdam.) 


ation behaviour as N atoms diftusMlong the b-axis in the crystal. 
paring (he two curves, it was proposed that removal of H,0 from the 
er sulphate crystals is a very fast diffusion process, so no dielectric 
ation is observed. Thus significant additional information is provided by 
lielectric measurements. 

tension of the above studies [18] to include another known order-disorder 
electric crystal, K4Fe(CN), . 3H,0, gave an endothermic DTA peak at 
5°C with accompanying mass loss showing removal of 3H,O. Curves of 
ctric constant against temperature show that prior to the dehydration, the 
ν molecules are involved in diffusion, similar to the order—disorder 
tition in NaNO,, but escape of the HO leads to an ordered anhydrous 


pound 


11.4 Thermosonimetry (TS) and thermoacoustimetry 
.1 Introduction 


lermosonimetry (TS), sound waves emitted by a sample are measured as a 
tion of temperature during a controlled heating programme, The sounds 
ted arise from the release of thermal stresses in the sample and range in 
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frequency from audio to several MHz. Stresses may be relieved by any of the 
thermal events described in Table 2.1, as well as by motion and creation of 
crystal defects. These last processes are not generally detectable by the more 
conventional thermal analysis techniques, on ὁ ccount of their low energy. and 
thus TS may be used, amongst other applications (section 11.4.4), to assess 
radiation damage, defect content and the degree of annealing of samples 

In thermoacoustimetry, the characteristics of imposed acoustic waves after 
passing through a sample are measured as ἢ function of temperature while the 
sample is subjected to ἃ controlled temperature programme, Further details of 
this technique are given briefly in section 11.4.3 


11.4.2 Apparatus for TS 


The sounds are emitted as mechanical vibrations prior to and during thermal 
events in the sample. This sonic activity in the sample is picked up and 
transmitted by means ofa specially adapted stethoscope. The mechanical waves 
are converted to electrical signals by conventional piezoelectric transducers 
The output consists of a rapid cascade of decaying signals, which may be 
recorded in three ways: (i) as the mean rate (counts per second); (ii) as the mean 
amplitude level (energy); or (iii) by direct frequency monitoring. 

The stethoscope is made of quartz (up 10 1000°C) or ceramics or noble metals 
for higher temperatures. The sample is held in the sample head which acts as an 
acoustic transformer and is connected ria a transmitting rod to a piezoelectric 
cell, fixed on a heavy recoil foundation and a scismic mount to prevent 
interference from external noise. A schematic diagram of the apparatus [19,20] 
is shown in Fig. 11.10. The propertics of the transducer may vary with 
temperature, so the waveguiding system is used to transmit the acoustic 
emissions from the heated sample to the transducer at ambient temperature. 
Contact surfaces must be well polished and thin films of silicone oil improve 
signal transfer. Care has to be taken that the switching of the temperature 
programmer does not register as noise on the transducer. Direct insertion of a 
thermocouple in the sample can cause severe mechanical damping, so the 
thermocouple is usually placed as close as possible to the sample without 
actually touching it. Alternatively a DTA head may be placed in the furnace for 
simultaneous TS-DTA measurements. Such a system has been used [20] to 
examine the transitions in the ICTA standards (Table 4.1). 

Frequency distributions of the TS signals are obtained [20] by timing the 
intervals between the amplitude components of the decaying signal. The time 
intervals are converted into pulse heights which are fed to a multi-channel 
analyzer to give ἃ display of the frequency distribution. Relationships are then 
sought between the frequency distributions and the processes occurring in the 
sample. 
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ire 11,10 Apparatus for thermosonimetry [19]. (With the permission of John Wiley 


Sons. Chichester.) 


4.3 Apparatus for thermoacoustimetry 


the apparatus used by Mraz et al. [21] a pair of lithium niobate transducers 
in contact (under a constant pressure of ~ 3 atm) with opposite faces of the 
aple (Fig. HEL), One transducer induces the incident acoustic signal and the 
ier detects the transmitted signal. Thermal expansion of the sample during 

heating programme is monitored continuously with a linear variable 
erential transformer {ΠΝ ΤῊ), so that changes in sample dimensions can be 
ved for in the calculations. Arrangements are made for atmosphere control 


wind the sample. 


tei 


Thermosonimetry (TS) and thermoacoustimetry 113 


(ἢ 


top support 
bar. 


insulator 
block 


top stack -- τς 
heater | aluminum 
contact_ 
driver 


aluminum 
spacer 


receiver 
. crystal 
ceramic}——4— - aluminum 


spacer τ --|pattom stack 
= heater 


base plate 


Figure 11.11 Apparatus for thermoacoustimetry [21] (With the permission of 
Thermochimica Acta, Elsevier. Amsterdam.) 


The incident signal is generated by a pulse generator. The transmitted signal 
received by the second transducer is inverted and amplified and an attenuated 
version of the driving pulse is added to this output signal. This summation 
procedure enables detection of the first-arrival times of both the compressional 
(P) and shear (S) waves. The velocities of the P and S waves and hence the 
various elastic moduli are then computed at set temperature intervals. and the 
final result is a plot of velocity or modulus against temperature, The instrument 
was calibrated with an aluminium standard for which the P and S wave 
velocities were accurately known, 


11.4.4 Applications of TS and thermoacoustimetry 


TS curves are usually used in combination with other thermoanalytical results 
to characterize a sample. TS curves on kaolins [19] show two regions of 
increased activity (Fig. 11.12), Comparing results obtained by TS with those 
obtained by TG and DTA, these regions have heen identified as dehydroxyla- 
tion (S00 600°C) followed by reerystallization (980 985°C) to the metakaolin 
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we 11.12 TS curve for a kaolin sample [19]. (With the permission of John Wiley & 
s, Chichester.) 


1200 


icture, TS examination of the thermal decomposition of a powdered sample 
srucite (Mp(OU),) [19]. suggested that dehydroxylation occurs ina stepwise 
rion with successive bursts of reaction as the sample breaks up. The barrier 
οἱ which has to be overcome was not identified. 

Che frequency distributions obtained [400] on heating samples of K,SO, (a) 
1 KCIO, (b) to their transition femperatures of 582°C and 299°C, 
pectively, are shown in Fig. 11.13. The pre-transition activity of KCIO, 
οἶνος release of included fluid, while that in K,SO, corresponds to 
eration of micro-cracks. The frequency distribution for K ,SO, shows fewer 
‘frequency components. Frequency distributions for the dehydration steps 
"USO, . 5H,O are qualitatively similar to those of KCIO, again through the 
d loss processes. 

Chermoacoustimetry has been used [21] to distinguish between grades of oil 
les. Both the P and the § wave velocities decrease with increasing 
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Figure 11.14 Thermoacoustimetry curves of oil shales [21]. (With the permission of 
Thermochimica Acta, Elsevier, Amsterdam) 
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Figure 11.15 Thermoacoustimetry curves of synthetic fibres [22]. (With the permission 
of John Wiley & Sons, Chichester.) 
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nperature and with increasing organic content (Fig. 11.14), Results are very 
sroducible, Discontinuities and peaks in the plots are related to loss of water 
d decomposition of some hydrocarbon fractions. 

Thermoacoustimetry has also been used [22]. in combination with DTA, to 
amine the characteristics of synthetic fibres (Fig. 11.15). The increases in 
imal occur firstly at the glass-transition temperatures and then prior to 
ching, The glass libre shows no changes in this temperature range. 
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Chapter 12 


The use of microcomputers in thermal 


is 


analy 


12.1 Introduction 


As with other fields, thermal analysis has been greatly influenced by the 
microcomputer revolution [1]. Wunderlich [2] has dealt with some of the more 
historical aspects of this impact. It is ofinterest (hat only just over ten years ago 
Wendlandt [3] reviewed the applications of digital and analogue computers in 
thermal analysis. Since that time analogue computers have faded from the scene 
and changes have heen so rapid that equipment described in papers published 
even a few years ago can easily be obsolete. In dealing with obsolescence. one is. 
of course, referring to the on-line use of computers for direct data capture and 
interactive control of the instrument, as distinet fram the more established off- 
line uses in data processing (Fig. 12.1). 

Many of the manufacturers of thermal analysis equipment now offer both the 
hardware and the software necessary to carry oul most thermal analyses with a 
high degree of automation, and to calculate the usual parameters from the 
captured data. As might be expected. these complete systems are not cheap, as 
prices include a component for software development. which is notoriously 
expensive. These systems are also specific to the particular manufacturer's 
equipment and, in some cases, the eoftware may be difficult, or even impossible 
to modify for one’s own requirements. Some calculations. such as the derivation 
of kinetic parameters under non-isothermal conditions [4], are areas of 
controversy and continuing development and it is necessary to know the 
approach and the algorithms being used in such software before any 
significance can be attached τὸ the output, 

For routine applications and/or in the absence of programming expertise, the 
computer systems provided by the manufacturers of the thermal analysis 
equipment are probably worth the additional expense. Some of these systems 
have been described in the literature, ¢.2- the DuPont 1090 system [5]. the 
Mettler TA3000 system [6], and the Perkin-Flner TADS system [7]. but the 
latest information is best obtained from the manufacturers. 
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Figure 12.1 Use of computers [13]: (a) off-line. (b) on-line. (With the permission of the 
Journal of Chemical Education.) 


For research purposes and/or when funds are not readily available but 
someone with an interest in computers is, an initial reduction in sophistication 
can probably be traded against increased flexibility, by coupling one of the 
many general-purpose microcomputers to whatever thermal analysis equip- 
ment is available. The prices of microcomputers have fallen to levels 
comparable with those for good chart recorders and their dedicated use should 
be casily justified in terms of the increased efficiency. 


12.2 Hardware requirements [8, 9] 


The hardware requirements for coupling of a computer to an instrument include 
the following aspects 
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(1) Suitable analogue signals should be available from the instrument for 
analogue-to-digital conversion. Modern instruments may have special 
outputs in the usual ranges of Oto SV, 0to 10 V or ~ Sto +5 V for most 
analogue-to-digital converters (ADCs) (Fig. 12.2(a)). The output signals 
from older instruments were designed for potentiometric recorders in the 
mV range and will need high-gain, low-noise amplification. Interfaces 
containing variable gain amplifiers (which may be software controlled) in 
addition to ADCs are commercially available [8]. Certain instruments. 
usually with digital readouts, may have some outputs, ¢.g. temperature, in 
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(a) Four-bit analogue-to-digital converter (ADC) 
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(b) Basic digital-to-analogue converter (DAC) 
Figure 12.2 Computer interfacing [13]. (From Perone, S. P. and Jones. D, O. (1973) 
Digital Computers in Scientific Instrumentation, McGraw-Hill Book Co., New York, 
with permission.) 
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hinary-coded decimal (BCD) form [9] to permit interfacing with other 
devices, including computers, Alternatively, the analogue signal may be 
converted to BCD through use ofa suitable digital voltmeter (DVM). The 
interfacing of a BCD output to a computer is described by Blanck [10]. 

When only linear heating programmes are used in the thermal analysis, ἃ 
direct reading of the temperature is not essential, provided that the time 
scale can be accurately related to the temperature (see (4) below). 

(2) The resolution of ADCs [8] depends on the number of bits in the digital 
output e.g. an 8-bit ADC may convert the maximum input, 5 V say, into a 
number 255 (binary | 1 1111 1). Maximum resolution is then 1 in 
256 or 19,5 mV. Variations of signal smaller than this value will thus not be 
detected. A 12-bit ADC gives ἃ more acceptable maximum resolution of tin 
4096, which is sufficient for most purposes, although 16-bit converters are 
available at much higher cost. For sequential capturing of several signals, 
an ADC with an adequate number of channels is required. 

(3) In most thermal analyses there is no demand for rapid (<1 5 interval) 
sampling of the signals, so a very fast (and hence more expensive) 
conversion time in the js range is not necessary. 

(4) Timing of sampling requires that the computer should have a real-time 
clock. In simpler systems, sampling may then be carried out by r@peated 
reading of the clock until the set time is reached, at which time that ADC is 
read, after which the program returns to reading the clock until the next set 
time is reached. In more sophisticated systems an interrupt facility is 
available which signals the computer at preset intervals that sampling 
should take place. For the remainder of the interval the computer is free to 
carry out other operations and thus may do real-time processing of the data 
already collected, or may service the requirements of other instruments. 

(5) For the control of instruments, for automation of analyses, the computer is 
required to do a calculation based on values of signals captured and to 
makea decision on the basis of the result (e.g. in ἃ ΤῸ run, no further loss of 
mass may be occutring). The decision, in the form of a number, then has to 
he translated into some kind of control action through use of a digital-to- 
analogue converter (IMAC) (Fig. 12.2(h)). The voltage output. correspond- 
ing to the number, may be used to operate a switch or valve, or to drive a 
stepping motor etc. For the DAC 8-bit resolution is usually ample, 

Interactive control of the instrument during an analysis may be arranged 
hy using one or more additional channels of the ADC to input signals, upon 
receipt of which action is taken. A channel of the DAC may be used to 
actuate an alarm, if necessary. 


It is seldom necessary to do any of the design or building of interfacing 
equipment, Anything required can usually be quickly assembled from standard 
commercially-available modules [8]. 
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12.3 Software requirements 


When it comes to considering software requirements, it is very important, as 
James [1] puts it, ‘not to reinvent wheels in computing’. Much software is 
available, admittedly of variable quality, but at prices that are low, through 
economics of scale, compared to the cost of development from scratch. The 
availability of a wide range of high-quality, general-purpose software should be 
one of the criteria in the original selection of the hardware. 

Because of the generally slow rate of sampling required, the program for data 
acquisition can generally he written ina high level language such as BASIC. For 
faster sampling. ASSEMBLER or MACHINE CODE programs are usually 
necessary. 

Again, because of the slow rate of sampling. some sort of averaging procedure 
can often be introduced at this stage to smooth the captured signal. Several 
rapid readings of the ADC output may be taken and the average stored 


12.4 Data storage 


Data will generally be initially stored as an array within the random-access 
memory (RAM) of the computer (the size of which will limit the number of 
points that can be collected) and must then be stored more permanently, either 
before and/or after further processing. Storage may be on magnetic tape or. 
more efficiently, on disk. An alternative which has been suggested [1 1]. is to link 
the dedicated microcomputer, via a communications interface, to a mainframe 
computer and to transfer the data to the larger computer to make use of its 
usually superior storage facilities, calculating ability and powerful software, as 
well as its peripherals, such as printer and plotter. This approach avoids the 
expense of having the dedicated microcomputer connected to ils own storage 
device (although this is needed for so many purposes), printer and plotter. 
During the storage process. all the relevant run parameters should be 
included for later reference. Back-up copies of important data are essential, 
especially when storage is on floppy disk. The ease and convenience of handling 
of data files may be a consideration in deciding upon the computer system to be 
purchased. (This is one of the weaknesses of Apple microcomputers.) 


12.5 Data processing 


The data processing required for thermal analysis in general, and especially for 
DSC and DTA, bears many resemblances to that developed for gas 
chromatography, c.g. the establishment of and correction of the signal for the 
baseline; detection of onset of peaks, neak maxima and the return to the 
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scline; resolution of overlapping peaks and determination of areas under 

aks by numerical integration [12, 13]. The routines required for TG usually 

dolve smoothing and numerical differentiation [14, 15]. buoyancy correc: 
yns, and possibly polynomial regression. All of the above procedures are well 

ycumented (c.g. ref. 16). 

Some programs in Applesoft BASIC are listed and described in Appendix B. 

gure 12.3(a) shows a noisy DSC trace (obtained using proffim DSC) after 
irrection for the rescanned baseline, using program OPERFILE. The trace 
as smoothed using one pass of program SMOOTHEFILE (5-point Savitsky 

olay filter) (Fig. 12.3(b)) and five passes (Fig. 12.3(c)). Program BINARY 

ILTER was also used for smoothing of the original trace. The results for one, 
ve and ten passes are shown in Figs 12.3(d), (e) and (f), respectively. The areas 
nder the peaks were determined using program TRAPSUM., and are given on 
ig. 12.3, to show the effect of smoothing on the measured area, Spikes on a 
‘ace are best removed manually or using program ROBUST FILTER. For 
tore details see Appendix B. Other more specialized procedures, such as 
inetic analysis under non-isothermal conditions, determination of purity by 
etailed analysis of the shape ofmeltingendotherms, the determination of glass- 
ransition temperatures and of heat capacities, may be developed to suit the 
ser's needs, General procedures such as graphics display and plotting routines 
re obviously dependent upon the specific equipment in use. 

The choice of programming language [17] will usually be some form of 
SASIC, because of the amount of software available, but a knowledge of 
FORTRAN is useful for implementation of published data processing routines 
16], and more and more programs are being written in the more structured 
anguages such as PASCAL, C, APL, ADA εἰς. [17]. Compilers are available 
or use with BASIC which give a considerable increase in speed. 


12.6 Additional benefits 


Once a microcomputer is available, even though much of its time may be 
dedicated to data capture, there will always be opportunities for use in 
additional (bonus) areas such as word processing, reference storage and 
possibly even searching of literature and data banks, as well as laboratory 
stock-control and computer-assisted instruction (CAT) [2]. 


12.7 Automation of thermal analysis equipment 


Thermal analysis equipment does not generally lend itself to full automation in 
the sense of analysis of sample after sample without operator intervention [18]. 
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sat advances are however being made in the naam of complex 
nipulations and developments can be expected in this direction. If there is 
ne means of changing the sample automatically, there is usually no problem 
arranging for automation of the heating/cooling programme and the 
itching of carrier gas as required. Most modern instruments allow for the 
rage, recall and implementation of such procedures. A multi-specimen 
rousel’ thermobalance has been designed by Ferguson ef al. [19] (Fig. 12.4), 
ich allows twenty specimens, mounted on arms radiating from a central 
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gure 12.4 Schematic diagram of multi-specimen carousel thermobalance [19], (With 
ὁ permission of John Wiley & Sons, Chichester.) 
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shaft, to be weighed sequentially. ‘There is obviously a limit to the sorts of 
treatment which the samples can reccive between weighings, but they can all be 
exposed to a similar environment. 

The Perkin-Elmer Corporation has recently announced a robotic system 
which enables 48 samples to be run consecutively on their DSC without 
operator intervention. The system, which is shown in Fig. 12.5, hasa removable 
48-position sample carousel and a pneumatically controlled sampling arm. The 
arm automatically selects any desired sample from the carousel, places it in the 
DSC sample holder and closes the cover. Control is ria the data station's 
microcomputer, so that once the sample is loaded, it is subjected to 
programmed treatment as desired. On completion of this treatment. the sample 
is removed from the DSC and replaced in the carousel, The sequence in which 
the samples are handled and their individual treatments can be programmed 
and the data obtained are automatically stored on disk. The system is obviously 
at its most efficient when used for quality control with one procedure applied to 
many similar samples. 


Figure 12.5 Perkin-Flmer DSC 4 Robotic System, (With the permission of Perkin- 
Elmer Corporation.) 
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Chapter 13 
Reaction kinetics from thermal 


analysis 


13.1 Introduction 


The rate of a general homogeneous reaction of the form 
A+B+C 


is conventionally measured by following the decrease in concentration of 
reactant A or the increase in concentration of cither product B or ( at constant 
temperature. A rate equation of the form 


Rate=k ficoncentrations of reactants and products) (T constant) 


is then determined from experiment. The rate coefficient, k, is a function of 


temperature: 
kane ΚΤ (δε k= ATC IM) 


and by carrying out a series of experiments over a range of different but constant 
temperatures, the Arrhenius parameters, F the activation energy and A the pre- 
exponential or frequency factor, can be determined. 

In thermal analysis the reactions studied are almost invariably heterogeneous 
reactions and the reaction temperature is usually being continuously increased 
or decreased according to some set (usually linear) programme. In spite of these 
differences, much effort has been directed at obtaining kinetic information from 
the results of thermal analysis experiments. Numerous papers have appeared 
and are still appearing on this topic and itis a field of considerable controversy. 
What follows is only an introductory and fairly selective account of some of the 
problems and the approaches suggested for their solution, More detailed 
accounts can be found in the references quoted. 


13.2 Heterogeneous reactions 


When a solid sample is heated, one of the many possible changes which it may 
undergo is decomposition (Table 2.1). Information on the kinetics and 
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inisms of solid decompositions is of both practical and theoretical 
tance [1]. For the heterogeneous reaction 


Als) B(s) + Clg) 
eg. CaCO,(s) +CaO(s)+CO,(g) 


yneept of concentration no longer has the same significance and the 
ess of reaction has to be measured in some other way. Usually the 
anal reaction, ας is defined in terms of the change in mass of the sample 
No <= m)/(mg <M) Where Mo is the initial mass and m, the mass of the 
le when reaction is complete) or equivalent definitions in terms of amounts 
sevolved or heat absorbed or evolved. The heterogeneity of such react ions 
“restricted (o the participation of the different phases represented above. 
+ are also the variations in properties with direction (anisotropy) in the 
al structure of a single pure solid, as well as the presence of numerous 
rities and structural defects, including surfaces. edges, dislocations and 
1 defects, in any real solid sample. Such factors have no analogy in 
ogencous reactions, but have marked effects on the physical properties of 
s, especially their thermal stability [2]. 

umerous observations confirm that decomposition of solid reactants 
rally is initiated at defective regions of the crystal such as the surface or, 
> specifically. points of emergence of dislocations at the surface. Nuclei of 
| product, B, are thus formed, the gaseous product escapes (sometimes with 
sulty) and the resulting disruption causes strain in the neighbouring regions 
jreacted A, resulting in growth of the nuclei (Fig. 13.1). The shape of these 
ci will he governed by the crystal structure, in that decomposition in some 
stions will occur more readily than in others. The detailed geometry of the 
esses of nucleation and growth [3, 4] leads to specific predictions of the rate 
hich the product gas is evolved, In addition to considering the geometry of 
eactant/product interface, account has to be taken of the chemical reactions 
ng place at or near the interface [4] and physical processes such as diffusion 
heat transfer [5]. 


13.3 Formulation of the problem 


inctic study thus involves measurement of acither asa function of time, f, at 
stant temperature, or asa function of temperature, 7, which is increased 
ording to some heating programme (usually linear), Φ -- 4ΤΊ αι, The 
hermal method, α rs. ἦν corresponds to the conventional concentration rs. ἢ 
ve of homogencous kinetics, while the dynamic method, i.e. measurement of 
+. Τ᾽. ἃς the basis of thermal analysis (Figs 13.2(a) and (b)). 

εἴποις analysis of both tsethermal and dynamic results involves attempting 
elate the experimentally observed a, For, Tvalues with values predicted for 


Formulation of the problem 129 


Figure 13.1 Formation and growth of nuclei of product in the decomposition of solids: 
(a) nucleation sites, (b) first nuclei formed, (c) growth and further nucleation, (d) overlap 
of nuclei, (¢) ingestion of a nucleation site, (f) continued growth. 


a limited set of models [1] based on processes of nucleation and growth, 
diffusion, or some simpler geometrical forms of progress of the reactant/product 
interface. The expressions derived from these models can all be written in their 
integral forms (at constant T) 


f(a)=k(t to) 
or differential forms 
da/dt—k g(a) 
(see Table 13.1 and Fig. 13.3). 
The effect of temperature is introduced through use of the Arrhenius equation, 
k= Ae~*/®T (or sometimes k= A Te TR! 16 R]) so that 


Fee 


ἀν αι τ Ae g(a) 
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Figure 13.2 (a) A series of isothermal a-time curves at different temperatures [56]. (b) 
The constant time contours derived from (a). (With the permission of the Journal of 
Chemical Education.) 


The validity of applying the Arrhenius equation to heterogeneous reactions has 
been questioned [5,9], but the parameters Eand A do have practical value even 
if their theoretical interpretation is difficult. 

For reversible reactions, the rate of reaction will depend on the partial 
pressure of the gaseous product. The rate equation should thus include 
allowance for this, in terms of some function h(p) 


da/dt = Ae~*'®T g(a) h(p) 
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Table 13.1. Broad classification of solid-state rate expressions 


1 da 
{(a)—kt MO= Fa 
1. Acceleratory « time curves 
P1 power law faa n(x)” 
El exponential law Ino ” 
2. Sigmoid a-time curves 
Α2 Avrami: Erofe'ev [Indl -α}}""}ἢ AL—ay--In(t ay? 
ΛΔ Avrami- Erofe'ev [ind |!" lay Inge)? 
A4 Avrami-Erofe'ev {Intl |" 41 -- αὐ Ing α}}" 
BI Prout- Τοπιρκίης Infa/(1 -0}} a(l-—o) 
3, Deceleratory a-time curves 
3.1 based on geometrical models 
Ε2 contracting area 1— (1a)! 2(1 —a)"!? 
R3_ contracting volume 1—( ~ay? 3(1-- αὐ 
3.2 based on diffusion mechanisms 
D1 one-dimensional σ΄ 1,2 
D2. two-dimensional (1 -αὐπί! - αὐτὰ —In(t-a)) ' 
D3 three-dimensional (1 (tay? (Vay 5} 


D4 Ginstling Brounshtein (1 27/3) (1 ay” {πῶ 5 ..} 
3.3 based on ‘order of reaction’ 


ΕῚ first order In(1-) 1--α 
Ε2 second order 1ΛῚ - a) ({--α 
ΕΔ ιμἰγά order [ΠῚ α]}} (1-- α)Ν 


This complication is usually ignored and this may he justifiable when working 
in vacuum or with a strong flow of inert gas through the sample. 
For dynamic measurements, the usual approach is to write 


do/dT =(do/dt) (dt/d7) = (1h) (da/dt) 


where @=d7/dt is the heating rate. The heating rate is usually maintained 
constant, although other programmes, c.g. hyperbolic with 1/7=u - rt where u 
and v are constants, have been considered [ 10. 11] and Rouquerol οἵ al. [12] 
have proposed heating the sample in such a way that reaction takes place al 
constant rate. 

The relationship between da/d7 and da/dr is one of the areas of controversy 
[13-17] but ifthe above procedure is accepted. and most treatments proceed on 
that basis, then 


dor /AT = (1p) (d/dt) = (Alp) ς΄ F/T g(a) (13.1) 
Separating variables 
ἀν (υ)- (4.φ)ς "TdT 


min) 
t(min) 


Figure 13.3 o-time plots for the model expressions (Table 13.2): (a) sigmoid group, (b) 


deceleratory group (based on «=0.98 at = 100). 
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The unknowns are thus A, ΠΕ ἀπά the form of either fla) or g(a). (Note th 
is no standardization on this terminology and in some papers {5} and 
used in exactly the opposite sense to that given here.) Often, f 
homogencous kinetics, g(a) is taken to be (1 - 0)". so that n. the apparer 
of reaction’, becomes the third unknown, A more general function, sugg 
Sestak and Berggren [18] is 


ayo" 0)" In =)" 


which increases the number of unknowns, but allows for all the usual πὶ 
solid state reactions. The unknowns then have to be determine 
experimental measurements which can be converted to values of α and/ 
at temperatures 7. obtained at a set heating rate ᾧ. The methods use 
analysis of the non-isothermal kinetic data are then usually divid 
differential methods, based on use of equation (13.1), or integral method 
on use of equation (13.2). 

A brief history of the development of non-isothermal kinetic analysis | 
given [19] and contributions to the field continue to appear in the 
literature, The International Confederation for Thermal Analysis (IC 
appointed a committee to investigate and report on the field 
introductory account only the main approaches will be outlined. Mo 
can be found in various reviews [1, 17, 19 21]. Some attention will be 
those approaches which have become ‘named’ methods and are often 
and sometimes incorporated in software packages (chapter 12). Such1 
will be amongst those to he fully scrutinized by the above-mentione 
committee [58]. Their mention here implies no more than that they hi 
widely used. Objections that have been raised to the use of these 
methods will be mentioned where possible. In addition, the more 
objections to the non-isothermal approach should be kept in mind 
these objections have heen summarized by Garn [5] and concern the p 
of measurement of sample temperature when allowance has to be made 
transfer from the furnace to the outer regions of the sample and then 
sample; the self-cooling or sell-heating of the sample during reaction. 
removal of evolved product gases from the vicinity of the sample 
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degree of reversibility. 

Attempts have been made to determine all three unknowns, A, Zand the form 
of {{ν} or a(x), from a single dynamic (i.e. «, 7) run. This is only possible if 
assumptions are made about the form of f(x) or g(a). It has often been assumed, 
without much justification, that reaction is analogous to a first-order process 
(q(a)= 1 -- α and fla) = —In(l —«)). 

If one is prepared to relax the constraint that data from a single run only are 
to he used, considerable simplification is possible. The three possibilities for two 
runs are that they should be: (i) both isothermal - which corresponds to 
conventional kinetic studies and, although it is probably the most reliable 
method, will not be discussed further here; (ii) one dynamic and one isothermal 
run; and (iii) two dynamic runs at different heating rates. Once more than one 
run is permitted, the approach is usually to eliminate the unknown form of the 
model f(x) or g(a) by comparing measurements made at a common value of « 
under the two sets of different conditions. 

In summary then, the approaches may be classified as either differential or 
integral methods and, within each group, as either single-run or multiple-run 
methods 


13.4 Differential methods 


These are methods based on equation (13.1) 
da/dT'=(A/p) e~F/®™ g(x) (13.1) 


A direct approach [22] has been to measure a series οἵα, 7 and (da/d7) values at 
a given heating rate ᾧ (with the practical problems of integration or 
differentiation and their associated errors, depending upon the initial form of 
the experimental data). Ifa form of g(a) is assumed, then from equation (13.1) 


y= In| (da/d7)/g(a)] = In(A/h)— (F/R) (1/7) 


so that a plot of y against (1/7) leads to values of F and A. 

Piloyan et al. [23] suggested the simplification that for values of «<0.5, q(x) 
could be taken as approximately constant and Καὶ and A could be determined 
from a plot of In(da/d7) against 1/7. 

An extension of this approach, when more than one set of data is available, is 
the method of Carroll and Manche [24] in which da/d7 is measured at a fixed 
value of α from a series of thermal analyses at different heating rates. ᾧ. Then 
In(¢ da/d7) is plotted against 1/7 and £ is calculated from the slope. 

The often-quoted Freeman. Carroll method [25] assumes g(x) = (1 —«)" and 
considers incremental differences in (da/d7), (1 ~) and (1/7) and leads to the 
expression 


"ον “αν ὧν 


————— 


te ee ih ee i nel Me ee SD 


This expression can be used to determine the value of E by plotting 
[A In(da/d7)/A In(t oJ against, [A (ΤΊ ΛΔ Int a) or [A In(ds 
A(1/T)] against [A In(1--%)/A(1/7)]. Sestak οἱ al. [19] have discussed 
improvements, including extension to other models. It is worth noting tha 
constant heating rate ᾧ. A In(do/d7)= A In(do/dt). 

Criado et al. [26] have demonstrated that the Freeman and Carroll trea 
does not allow an nth order kinctic equation to be distinguished succe: 
from other kinetic models. This was confirmed by Jerez [27] who point 
the large errors involved in the regression procedure, and sugges 
modification in the calculation of / and n involving the use of the point at 
the rate is a maximum and the centre of gravity of the cluster of experir 
points. 

Van Dooren and Miiller [28] have examined the effects of sample 
particle size and heating rate on the determination of n, Zand A from DS! 
using the Freeman and Carroll procedure. Plots using points from the σοὶ 
DSC peaks were curved. so different sections of the DSC peak were cons 
separately. Different values of the kinetic parameters were obtained { 
different section of the peak and these values also varied with sample 
particle size and heating rate. 

Several approaches involve use of the second derivative of equation (1: 
the version with g(a)—(1 —a)", with respect to temperature [19, 20] ¢ 
respect to α [20], in spite of the problems of obtaining accurate values of: 
derivatives. Using 

da/dT=(A/p) cP"? (1 - αν" 


τῇ E (9) nm 
dT RT? 4171 - τὶ 


gives 


and since this derivative must be zero at the inlexion point ofa TG curve 
maximum of a DSC peak 
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from which E may be calculated ifn is known and Tyg (12/47 inne anda 
measured [29]. Combining equations (13.3) and (13.4) gives 
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and since (1—0%m,,) is a constant for a given value of n, the Kissing: 
method of obtaining a value for Fis to plot Ππ Φ ΤᾺ,.} Against V/Thyae FOF: 
of runs at different heating rates, ᾧ. The slope of such a plot is -- F/R. 

Augis and Bennett [31] have modified the Kissinger treatment for τι 


the Avrami--Frofe’ev model, which applies in so many solid-state rea 
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They plot In(G/(Tyrny -10}} AGAINST T/Fiyay WHETE 10 18 The INTUAT Temperature at 
the start of the heating programme, instead of In(/72,,,) against Τ 7, κι: 

Hlder [32] has generalized the Kissinger treatment to make it applicable to 
the full range of kinetic models, The generalized equation is 


In(h/TM2)=In(AR/E) + In (L~ E/RT yy) 


max 


where mis the temperature exponent of the pre-exponential term in the modified 
Arrhenius equation and is often taken as zero, and 


L= —Q(Omay)/(1 4+ κι ΚΤΊ κε  ΕῚ 


This correction term was found to be relatively small, but helps in distinguishing 
between similar models, The values of F obtained were not very sensitive (o 
incorrect choice of model. 

Although originally derived as an integral method, the Ozawa treatment [33] 
is also applicable to derivative curves and is similar to the Kissinger method. In 
the Ozawa method In ¢ is plotted against 1/7;,,, and the slope of this plot is 
again -- F/R. 

Van Dooren and Miiller [34] studied the effects of sample mass and particle 
size on the determination of kinetic parameters from DSC runs using the 
methods of Kissinger and of Ozawa. It was found that both sample mass and 
particle size could influence the values of the kinetic parameters, but the extent 
of these effects varied from one substance examined to another. The two 
methods gave similar values for E with slightly lower precision for the Kissinger 
method. It was suggested that temperatures at =0.5 (half conversion) should 
be used in place of T),9.: 

The Borchardt and Daniels method [35. 36], originally developed for DTA 
studies of homogencous liquid-phase reactions, involves calculating the rate 
coefficient, k, from the expression [37] 


k={(KSV/mp)""! (C dAT/dt-+ K ATYV/[(K(S—s)—C AT)]" 


where V is the volume of the sample and mp the initial amount of sample in 

moles, Sis the total DTA peak area, sis the partial area up to the times at which 

the peak height 7 and the slope dA7/di are measured, K is the heat transfer . 
coefficient and C is the overall heat capacity of the sample and the holder. 

Because of the problems of obtaining values of K and C, the assumption that 

n= land the approximation [36] that C(dA7/dt) and CAT are small compared 

to the terms to which they are added. or from which they are subtracted, are 

made, so that 


k=AT/(S~—s) 
For DSC the equivalent form would be 
k=(dH/dt)/(S~s) 
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Shishkin [38] has discussed the similaritics between the Borchardt and Dar 


and the Kissinger approaches. 
Another method based ona second derivative is that of Flynn and Wall [ 


They differentiated equation (13.3), rearranged as 
Tda/AT)=(AT?/p) € FIRT (1 .- α)" 


with respect to a, to give 


Gh (er eo : n da 
a (7 aa) <crims2r( τ (aun) 


At the beginning of decomposition, i.c. low «, the last term on the right-| 
side is negligible and even 27 <E/R. ‘They then took finite differences inste 


the derivative, so that 
A(T? da/dT)/Aa=(E/R)+ 27 
where T was the average temperature over the interval. The value of E was 


calculated from the slope of a plot of (7227/7) against Ac for the beginni 
the reaction only, or, ifone starts from σ᾽ = Owhere da/dT=0, T?(do/d7T)m 


plotted against «. 
Some of the more common differential methods are summariz¢ 


Table 13.2. 


13.5 Integral methods 


Use of equation (13.2) 
f(n)= i, (Πα) «ταν 
0 


obviously involves evaluation of the temperature integral. Usually the lit 
integration are changed to fe "/* 4Τ ἵΓα is zero up to To. The proh 
simplified by introducing the variable x= £/RT so that 


Τ ” 
| oF ar=(E/R) | (e*/x?) dx =(E/R)p(x). 


0 


ε B/RT dT 


τ 
r, 


Then equation (13.2) becomes 
f(=(A E/Re)p(x) 


Tables of values of the integral p(x) have been provided [40, 41]. A 
attention has been directed at finding suitable approximations for the 
temperature integrals [1. 19, 20, 42, 43]. Gorbachev [44] has suggeste 
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the nature of the original ας 7 data. 

Zsak6 [20] has suggested sub-classification of integral methods on tt 
of the means of evaluation of the temperature integral. The thre 
approaches are; (i) use of numerical values of p(x); (ii) use ol 
approximations for p(x); (iii) use of approximations to obtain an exp 
which can be integrated. The series most commonly used for approxi 
p(x) are [20] the asymptotic expansion : 


p(x) =(e7 =x?) [1 —(28/x) + (1 γι) = (4x8) + Ἐ{- ΠΝ. Ὁ Nx) 
or, with better results [45], the Schlémlich expansion 
p(x) =(e>*/x(x +1) [1 = (1/(x + 2)) + (2/(x + 2x +9) 
—(4/(x +2) (44-3) Ὲ4)....7 


Taking only the first two terms of each series 


p(x)X(x—2)e */x? 


or 
plx)re > “/Oxte + 2)) 


Other approximations are the empirical formula [46] 
p(x)xe*/((x~d)(x-~ 2)) 
with d=16/(x? —4x +84), or [47] 
p(x)ve */(x(x?-+4x)"?) 
Doyle's approximation [48] (for x > 20) is 
log op(X)® ~ 2.315 —-0.4567% 


Coats and Redfern [49] suggested an asymptotic expansion of the temp 
integral, giving 

fla) =(ART?/PE)[1—(2RT/E)] ee PRT 
which, assuming 2R7/F «1. leads to 


In fla) =In( A R/PE) +2 nT (E/RT) 
or 
In( f(x)/T?)=—In( AR/HE) -- (ΕΠ ΚΤῚ 


One approach for arriving at an integrable expression [20. 43] is ta 
expansion for 1/Tin equation (13.2), with the aim of obtaining a simple Γ 
of T, (7) in equation (13.2) rewritten as [43] 


In fia)=In Cy + C, OCT) 
where C, is a constant containing 4 and C, a constant containing FE. 
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Van Krevelen ef al. [DU] usea 
(1/T)=(1/T,) [1 -- (7 -- T/T, + °°] 
while Horowitz and Metzger [51] used 
(1/T)=(1/T,) [1 —In(7/T,)+ °° '] 


where Τὶ is an arbitrarily defined reference temperature somewhere within the 
extremes of temperature covered by the experiment [43]. Van Krevelen’s 
treatment leads to C,=E/RT,+1 and @(7)=In 7, while Horowitz and 
Metzger’s treatment gives C, = E/R(T.)? and 0(7)= Τ᾿ The problem with these 
approaches is the arbitrariness of the reference temperature. Use has been made 
of the temperature at which the reaction rate is a maximum, Kassman [43, 52] 
has shown that equation (13.2) can be made integrable by using 0(7) = 1/7. This 
leads to the integrated form 


Slo) =[Ae?(T,)?/(@((E/R) + 27, )ye 277 


hence 
In f(a) = —((E/R) + 27,)/T + constant. 


Kassman also showed that the reference temperature can be optimized in the 
above expression, or those of Van Krevelen or Horowitz. and Metzger above, by 
calculating the geometric mean of the temperature interval 


Tn = (ΤΊ ονενι Rigen) 


and then using 
T= To, 01 = (Tyen/ C1) + 2(Tyme/C1)7] 


where C, is the observed slope of the fitted line when In f(a) is plotted against 
1/T. Kassman [43] claims that this method is superior to the others in 
determining values of & and A. He also examined the problems of 
discriminating amongst the models, f(x), and showed that it was not possible to 
distinguish models within the various groups in Table 13.1 from a single non- 
isothermal data set. 

Zsak6 (20, 53] suggested the following practical procedure starting from 
equation (13.5) in the form 


In( f(x) —In(p(x))=In(AE/RG)= B 


where Bisa constant fora given reaction and constant heating rate, @. Various 
functions of a were assumed and combined with selected values of E. For each 
experimental « value, measured at temperature 7, one thus obtains a value of 
f(a) and hence In(f,(«)). Combination of the corresponding experimental Τ᾽ 
value with an assumed E value= E,, provides a value of x, = £,/R7 and hence 
p(x,) and In(p(x,)). Tf the choice of both the model, f;(a) and the activation 
energy, £,, was correct 


In(f\(a)) ~ In(p(x,)) should be a constant, =In(AE,/Ro), 


of 


ar, "94 


using computers, such a search is time consuming. 

Satava and Skvara [54] developed a graphical method to simplify the sea 
Tables of log(/(a)) for the various models are available [40], or can be rea 
generated, together with tables of log(p(x)) against 7 for fixed values of E 
40]. The experimental α values are thus converted to log( f(a) values for var 
functions of «, and are then plotted against the experimental Τ᾽ values, or 
same scale as plots of —log(p(x)) against T for various values of E (Fig. 1 
One of these sets of plots should be on transparent paper. With the 7s: 
matched, the plots are moved relative to cach other in the y-direction ἀπε] 
of the log( {1α}} against 7curves matches one of the ~-log(p(x)) against Tcu 
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Figure 13.4 The plots required for the Satava method [3]. One plot is usually d 
transparent paper so that they can be superimposed for matching [56]. (W 
permission of the Journal of Chemical Education.) 


The displacement between the y-scales is then equal to B, from 
log(/(a)) —log(p(x))=B 

Hence the value of A can be calculated. The drawbacks are: (i) the difficulties in 
distinguishing between models [4], and (ii) the use of discrete E values. 

Relaxation of the constraints to allow data from one dynamic and one 
isothermal run, introduces considerable simplification. Under isothermal 
conditions (T= Tp), «=o, at t=t, (Fig. 13.5) and f(x)= 4 ε΄ ΚΑ τος. 

For the dynamic run at heating rate @, x-0, at T=7, and f(a%)= 
(4 ΕΙΚΦῚ p(x)) where x, = F/RT,. Hence 

Ae PRTot = (AE/Rp) p(x;) 
or 
t= (E/RO)e*Fp(x,) 


and 
log t, = log(E/Rp) + (F/2.303 RT)) + log(p(x,)) 


Evaluation of p(x,) requires an estimate of E, but using the Doyle [40] 
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Figure 13.5 An isothermal curve for the contracting volume expression, R3 [56]. (With 
the permission of the Journal of Chemical Education.) 
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approximation that log(p(x))= 2.519 —U.490/% (UW x > ay δὴν tie wre 
ture interval scanned in the dynamic runs is restricted to «100 K) 


log t= log(E/R¢) + (F/2.303 RT) ~ 2.315 — (0.4567 E/R\1/T)=m(1/T) 


A plot of log t, against (1/7)) should thus be linear (Fig. 13.6) with 
m= —0.4567E/R, from which a value of Ecan be calculated. Thus at the ς 
an additional run it is possible to eliminate assumptions about the form ¢ 
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Figure 13.6 Doyle's method [5] of plotting log ἡ, against 1/T, (see text) [56]. (Ὁ 
permission of the Journal of Chemical Education ) 


Ifinstead of using data from one dynamic and one isothermal run, wet 
dynamic runs at different heating rates p, and #,, then, following Gyu 
Greenhow [55], equation (13.2) is written as 

Υ 1 
s10)=(4/6) | extre ara | e FRT AT =(A/p)h 
To 0 
where 
τ τὺ 
b= [ ec FRT AT and | ε΄ 7.0 
0 0 


Jo 
The method is again based on comparison of runs at common values 0 
At heating rate @, ,¢—=a,at T= Τὶ but at heating rate ,,a=0,atT=T, 


Sf (,)=(A/b,)hy =(A/P 2) 2. A whole series of h,/h ratios is possible 


'y 


2. ene ae 


OO00O000 


{) 


ΕΥ̓ αν 


a> 


= 


- 


cy 
Sede μον οἰ τξν. ἃ τὸς 


value of αὐ has not been specified, i.e. (b,/h2)=(b,/b2). Tables of —log ΒΊΟΥ 
various values of Tand £ are given [55]. Pairs of corresponding temperatures 
(7, and 7), at common values of a, are determined from the two dynamic runs 
(Fig. 13.7). For some value of E=E, say, then combination with (7, and 71); 
results in values (h, and b,),,. Log(b2/h, ), is then plotted against δὶ (Fig. 13.8) 
and the unknown value of E is determined from the point at which 
log(h,/h,)—logthy/h, )y- 
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Figure 13.7 Non-isothermal curves for the contracting volume expression, R3 [56]. 
Heating rates: curve A 1,000 K min™!, curve Β 0.257 K min~ '. (With the permission of 
the Journal of Chemical Education.) 


13.6 The influences of various parameters on the shapes of 
theoretical thermal analysis curves 


With there being so many parameters in equations (13.1) or (13.2), it is useful to 
sce what effect varying onc of these parameters at a time has on the shape of a 
TG or DSC curve | 19, 20, 53]. The most fundamental variable is the form of the 
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Figure 13.8 Gyulai and Greenhow’s method [6]. Curves vary slightly (a- d) with val 
a chosen [56]. (With the permission of the Journal of Chemical Education.) 


mechanistic model fl) or σία). The shapes of the theoretical isothermal α 
curves for various models [1, 56] are shown in Fig. 13.3. Itis not always an 
task to distinguish amongst the models, even under isothermal conditions 
The shapes of these curves are, of course, considerably altered under 
isothermal conditions and theoretical TG curves [20, 32] for various mode 
given in Figs 13.9 and 13.10, The models based on apparent order of reactiv 
(even if fractional) i.c. Fl, F2, R2and R3, Fig. 13.9, are difficult to distingu 
low values οἵα. Distinguishability improves for higher orders at higher vali 
a. The diffusion models, D2, D3 and D4, give generally lower 
temperatures and flatter curves (Fig. 13.10) than the nth-order group, whi 
Avrami-Erofe’ev models have higher onset temperatures and steeper ct 
Figure 13.11 shows the differential curves corresponding to the integral c 
given in Fig. 13.10. 

The next stage is to choose one model and use this to examine the influe! 
the other variables, the heating rate @, the activation energy, &, and thi 
exponential factor, 4. Zsak6 [20] has done this for the first-order (Fl) πὶ 
Although this is nota very realistic model it is often assumed to apply as 
approximation. Figure 13.12 shows the regular effect on the theoretical 
order TG curve, of doubling the heating rate in the range ᾧ = I to ΙοΚπ 
Decreasing the pre-exponential factor by factors of 10 in the range A= 1 
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Figure 13.9 Influence of reaction order, n, on the shape of theoretical TG curves [20]. 
Ε-- 142 ΚΙ mol~',A=10'*s>',6=4K min '. Order m:a=0; b=1/3; c= 1/2; d= 2/3; 
ἐπι i f=2. (From Zivkovic, Z. Ὁ. (Ed.) (1984), Thermal Analysis, University of Beograd, 
Yugoslavia, with permission.) 
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Figure 13.10 Integral curves for various models [32]. ¢=10K min™', E= 


270 ΚΙ mol ', 4= 10'S min’ '. 1 to 9: models ΕἸ, A2, A3, R2, R3, D2, D3 and D4, 
respectively. (With the permission of Wiley Heyden Ltd.) 
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Figure 13.11 Differential curves for various models [32]. @=10K min 
220 ΚΙ mol”, A=10'!S min‘. | to 9: models ΕἸ. A2, ΑΔ, ΕΖ, R3, D2, ΠῚ anc 
respectively. (With the permission of Wiley- Heyden Ltd) 
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Figure 13.12 Influence of heating rate, ᾧ, on the shape of theoretical TG curves 
n=1, Ε-- 142 ΚΙ mol ". 4=10'*s '. Φ(Κ min '): a=l; b=2; c= 4: d=: ὁ 
(From Zivkovic, 7. D. (Bd.) (1984) Thermal Analysis, University of Beo 
Yugoslavia, with permission.) 
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Figure 13,13 Influence of pre-exponential factor, 4, on the shape of theoretical TG 
curves [20]. n=l, F=142kJ mol', φ -- 4 Κὶ min ἈΠ ἢ; a=10'*; b= 10"; 
e—10'4; d=10'3; e=10'?. (From Zivkovic, 7. D. (Fd.) (1984) Thermal Analysis, 
University of Beograd, Yugoslavia, with permission.) 
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Figure 13.14 Influence of activation energy, £. on the shape of theoretical TG curves 
[20]. n=1, A=10'48°', p=4K min |. Ε (kJ mol” "): a=126; b= 134; c= 142; 
{τι 1541: ¢=159. (From Zivkovic, Z. Ὁ. (Ed.) (1984) Thermal Analysis, University of 
Beograd, Yugoslavia, with permission.) 
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of the curves (Fig. 13.13) with the remaining segments being almost par 
Very similar behaviour is observed for curves with the activation er 
increasing in steps of 8.3 kJ mol! (2 kcal mol ') as shown in Fig. 13.1¢ 
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mechanistic model applying, while the position of this curve on the temper: 
axis is determined by the values of F, 4 and, toa lesser extent, the heating ra 


Figure 13.15 
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Schematic representation of steps involved in kinetic analysis of soli 


reactions [3]. (With the permission of John Wiley & Sons, Chichester.) 
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13.7 Conclusion 


Galwey [3] has given a flow diagram (Fig. 13.15) summarizing the steps 
involved in the kinetic analysis of solid-state reactions. Only non-isothermal 
methods have been discussed here, and the methods have been selected to 
illustrate the approaches that may be used. Numerous variations and 
modifications have been and are being published. More information will be 
found in refs 1 and 7-9. Garn [5] has provided a comprehensive and critical 
review of the many factors which must be taken into account when assessing the 
reliability of kinetic parameters derived from dynamic measurements, Two 
reasons given [5] for possible failure of a method of kinetic analysis are: (i) the 
approximations used may not be valid, and (ii) the model used in deriving the 
mathematics may not take into account physically real factors of the experiment 
such as the heat transfer and atmosphere control problems. 


13.8 References 


1. Brown, M. E., Dollimore, D. and Galwey, A. K. (1980) Reactions in the Solid State, 
Comprehensive Chemical Kinetics, Vol. 22, Elsevier, Amsterdam. 
2. Boldyrev, V. V. (1986) Thermochim. Acta, 100, 315. 
3. Galwey, A, K. (1982) Proc. 7th Int. Conf. Thermal Analysis, Wiley-Heyden, 
Chichester, Vol. t, p. 38. 
4. Galwey, A. K. (1985) Thermochim. Acta, 96, 259. 
5. Garn, P. Ὁ. (1972) Crit. Rev. Anal. Chem., 3, 65. 
6. Dollimore, D., Gamlen, G. A. and Taylor, T. J. (1982) Thermochim. Acta, 54, 181. 
7. Petty, H.R., Arakawa, E. Τὶ and Baird, 1. K. (1977) J. Thermal Anal., 11, 417. 
ἃ, House, J. E. Jr., (1981) Thermochim. Acta, 48, 165, 
9, Arnold, M,, Veress, G. E., Paulik, J. and Paulik, F. (1982) Thermochim, Acta, 52, 67, 
(1981) Anal, Chim, Acta, 124, 341, 
10. Varhegyi, G. (1982) Thermochim. Acta, 57, 247. 
11. Popescu, C. and Segal, E. (1982) J. Thermal Anal., 24, 309. 
12. Reading, M., Dollimore, D., Rouquerol, J. and Rouquerol, F. (1984) J. Thermal 
Anal., 29, 775. 
13. Gyulai. G. and Greenhow, E. J. (1974) Talanta, 21, 131. 
14. MacCullum, J. R. (1982) Thermochim, Acta, 53, 375. 
15. Garn, P, Ὁ. (1974) J. Thermal Anal., 6, 237. 
16. Sestak, J. and Kratochvil, J. (1973) J. Thermal Anal., 5, 193. 
17. Sestak, J. (1984) Thermophysical Properties of Solids, Comprehensive Analytical 
Chemistry, Vol. XUD, Elsevier, Amsterdam. 
18. Sestak, J. and Berggren. G. (1971) Thermochim. Acta, St. 
19, Sestak, J.. Satava, V. and Wendlandt, W. Ν᾿. (1973) Thermochim, Acta, 7, 447. 
20. Zsak6, J. (1984) Thermal Analysis, (ed. 7. D. Zivkovic) University of Beograd, Bor. 
Yugoslavia, p. 167. 
21. Blazejowski, J. (1981) Thermochim, Acta, 48, 109; (1984) 76, 359. 
22. Achar, B. N., Brindley, G. W. and Sharp, J. H. (1966) Proc. Int. Clay Conf. 
Jerusalem, 1, 67; Sharp, 1. H. and Wentworth, I. A. (1969) Anal. Chem., 41, 2060. 


πὸ IE ee ο΄ Te me 


1229. 


. Carroll, B. and Manche, E. P. (1972) Thermochim. Acta, 3, 449, 
. Freeman, E. S. and Carroll, B, (1958) J. Phys. Chem., 62, 394; (1969) 73, 751. 


Criado, J. M., Dollimore, D. and Heal, G. R. (1982) Thermochim. Acta, 54, 1S! 


. Jerez, A. (1983) J. Thermal Anal., 26, 315, 

Van Dooren, A. and Miiller, B. W. (1983) Thermochim. Acta, 65, 269. 

_ Fuoss, R. M., Sayler, O. and Wilson, H. 5. (1964), J. Polym. Sei., 2, 1147. 

. Kissinger, H. E. (1956) J. Res. Nat. Bur. Stand., 57, 217; (1957) Anal. Chem., 


1702. 


. Augis, J. A. and Bennett, J. FE. (1978) J. Thermal Anal., 13, 283. 
. Elder, J. P. (1985) J. Thermal Anal., 30, 657; (1984) Analytical Calorimetry, Vol 


(eds. P. 5. Gill and J. F. Johnson) Plenum, New York, p. 269. 


. Ozawa, Τὶ (1965) Bull. Chem. Soc Jpn., 38, 1881: (1970) J. Thermal Anal., 2,3 
"Van Dooren, A. A. and Miiller, B. W, (1983) Thermochim. Acta, 65, 257 

. Borchardt, H. J. and Daniels, F. (1957) J. Am, Chem. Soc., 79, 41. 

, Reed, R. L., Weber, L. and Gottfried, B. 5. (1965) Ind. Eng. Chem. Fundamental. 


38. 


. Sharp, J. H. (1972) Differential Thermal Analysis, Vol. 2, (ed. R. Ὁ. Macken 


Academic Press, London, p. 47. 


. Shishkin, Y. L. (1985) J Thermal Anal., 30, 557 

. Flynn, J, Η, and Wall, L. A, (1966) J. Res. Nat. Bur, Stand., 0A, 487. 
. Doyle, C. Ὁ. (1961) J. Appl. Polym. δῖ. 5, 285. 

_ Ζεακό, J. (1968) J. Phys. Chem., 72, 2406 

. Blazejowski, J. (1981) Thermochim. Acta, 48, 125. 

. Kassman, A. J. (1985) Thermochim, Acta, 84, 89. 

. Gorbachev, V. M. (1982) J. Thermal ἐπα], 25, 60% 

. Doyle, C. Ὁ. (1965) Nature (London), 207, 290 

. Zsako, J. (1975) J. Thermal Anal., 8, 593. 

. Balarin, M. (1977) J. Thermal Anal., 12, 169. 

. Doyle, C. Ὁ. (1962) J. Appl. Polym, Sci., 6, 639 

. Coats, A. W. and Redfern, J. P. (1964) Nature (London), 201, 68; (1965) Po 


Lett., 3, 917. 


. Van Krevelen, D, W., Van Heerden, C. and Huntjens, F. J. (1951) Fuel, 30, 7 
᾿ Horowitz, H. H. and Metzger, G. (1963) Anal. Chem., 35, 1464; (1963) Fuel, 42, 
. Kassman, A. J. (1980) J. Thermal Anal., 18, 199 

~ Zsak6, J. (1970) J. Thermal Anal., 2. 145; (1968) J. Phys. Chem., 72, 2406. 

. Satava, V. and Skvara, F. (1969) J. Am, Ceram. Soc,, 52, 591; Satava, V. (197 


Thermal Anal., 5, 217. 


. Gyulai, G. and Greenhow, E. J. (1974) J. Thermal Anal., 6, 279; (1973) Therma 


Acta, 6, 239. 


. Brown, M. E. and Phillpotts, C. A. R, (1978) J. Chem. Educ., $5, 556. 
. Brown, M. E. and Galwey, A. Κ. (1979) Thermochim, Acta, 29, 129. 
. Flynn, J. H., Brown, M. E. and Sestak, J. (1987) Thermochim. Acta, 110, 101 


Chapter 14 


Purity determination using DSC 


14.1 Introduction 


Measurements of the depression of the melting point of a sample are often used 
to determine its purity. Calculations are based on the assumptions that solid 
solutions are not formed and that the melt is an ideal solution. The practical aim 
of purity determinations is usually to decide whether or not the sample meets 
certain specifications, determined by the intended further uses of the sample. 

The melting endotherm for a pure substance recorded ona DSC is illustrated 
in Fig. 14.1. Tp is the melting point of the sample and the area ABC is 
proportional to the enthalpy of fusion, AH/,, of the sample. The presence of an 
impurity in the sample (the solvent) gencrally lowers the melting point of the 
solvent and also broadens the melting range, giving a broader DSC endotherm 
as illustrated in Fig, 14.2 (inset), From endotherms such as illustrated in 
Figs 14.1 and 14.2, melting points and enthalpies of fusion may readily be 
determined, With more effort an estimate of the purity of a compound can be 
obtained, from analysis of the detailed shape of its melting endotherm, e.g. 
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Figure 14.1 Idealized DSC record of melting of pure indium, Slope AB (=1/R,) is used 


to correct for thermal lag. Τὴν is the melting point. Area ABC represents the enthalpy of 
fusion, AI. [6]. (With the permission of the Journal of Chemical Education.) 
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Figure 14.2 Idealized DSC record of melting of an impure sample. The slope 
(=slope of AB in Fig. 14.1) is used to correct the programmed temperature, 1. 
sample temperature. 75. Το is the melting point. The area ABC represents the enth: 
fusion, ΑΙ Γι . The fraction melted, /,, at temperature 7, = area ADF area ARC 
[6]. (With the permission of the Journal of Chemical Education.) 


Fig. 14.2, without reference to compounds containing known amou 
impurities. 


14.2 Phase equilibria 


The simplest system to consider is that in which the impurity does not! 
solid solution with the solvent, but forms an ideal solution in the melt 
eutectic system, If the impurity is labelled, in (he customary way, as comp 
2 and the solvent as component 1, then for equilibrium (at constant pre 
between pure 1 in the solid and | in the solution (or ‘melt’) at activily a, 
must be equality of the chemical potentials 6} of | in the two phases’ 
wp (s)— 1) 

i.€,, pe(ss= rr (+ RT In ay 
(where the superscript “ refers to standard conditions, i.e. unit acti 
Differentiating equation (14.1) with respect to temperature, G 
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G45 (8) HN) Rina, + RT ay) 

47 dT dT 


and, since djt/d7= — 5 (where the bar represents a molar quantity) 
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Rin αι τίμ (5) <1 (DVT 
so equation (14,2) on rearranging, becomes 
d(Inay)_ = Car (s)= 10 (1 SPs) $7 0) 
dT RT? RT 
Us TSO =Car) + TST 
7 RT? 
or 
din ay) _ AP) Η 70) AH? 
47 RT RT? 
ace ΞΟ TSANG =). 5 Integrating equation (14.3) between the limits 
αν, τοὶ at Τ τε Το (since solid solutions are not formed) and a,=a, at T= Ἧς 


assuming that AH, is independent of temperature over the range, 


ΔΗ 1 1 
νον αὐτῶ (7-7) 
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For an ideal solution, a, =%, (the mole fraction of 1). Hence 


AHe, | ! 
In x, =In(l -- χγ) = ΓΝ (5 = i) 


(14.3) 


For a dilute solution, i.e. small values of x>, 


In(1—x,)> τ Χ2 
ΔῊΝ [Ὁ 
ie (᾽ x) (14.4) 


Equation (14.4) forms the basis of melting-point depression calculations, as 
follows. At T=T7,,. the melting point of the impure sample 


AH? ( - τὶ _ Ae, AT 


— 14.5 
R Vola {π|9] 
If AT; is small, Τὸ -- Το and ToTm% Τὰ, Also x2="3/(ny +n,)~m M,/1000, 
where m is the molality of solute and M, the molar mass of the solvent. Hence 


2M 
ΕΣ a jae k m 
AH?, .1000 


where Αἱ is termed the cryoscopic constant. 
Only when the sample is completely melted, i.e. at ΤῸ Tn: is the mole fraction 
of impurity in the liquid, x), the same as that in the original sample, x}. From 


value x3 is the minimuUM vatuc wren x7 sess π AE τοι» 
fraction of sample that has melted, F, is less than unity, the composition οἱ 
melt is closer to that of the eutectic, i... %2> x. When melting commences 
first liquid has the eutectic composition 
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Figure 14.3 Low concentration region of a simple cutect ic phase diagram (inset). 
lever rule: BC/AB=Nyonia/Miguia: © is the composition of the melt in equilibriun 
pure solid at 7. Here x, >}. F=n,/(n, ἐπ τ AB/(AB + BC)=AB/BC = v3/x, 
the permission of the Journal of Chemical Education.) 


If Fis the fraction melted at temperature 7, then, assuming ἃ linear 
segment of the liquidus curve (Fig. 14.3), and using equation (14.4) 
pee eee RTS _ 
% lon tT ΔΗ (- ἢ 
Rearrangement yields 
παρ S88] 
Ne le 
If Fcan be determined at various temperatures, T.a plot of Tagainst ΠῚ 
yield a straight line, provided that AH, is independent of temperature 
values of Tp and ΔΙῚ τ, are known, x can he determined from the me: 
slope of the line. The DSC curve is capable of providing values ὁ 
temperatures 7 for use in such a plot. 


14.3 The DSC melting curve 


The DSC measures the thermal energy per unit time, dH/dt, transferre: 


from the sample as the temperature of the sample holder, 7, is chang: 


SPS 


τὸ τ 


== 


constant rated V/dt = >. Thus the output from the DSC is directly proportional 
to the heat capacity of the system, 471,6 7° 
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ae (ar) ar, τὴν 
For an absolutely pure compound with zero melting range, d/7/dt would 
hecome infinite at the melting point, Το. For an impure compound, d///dT is 
finite and is a function of T. 

When the fraction melted, F, is zero, the heat capacity of the sample is that of 

the solid mixture, and when /=1 the heat capacity of the sample is that of the 
ideal solution, Intermediate behaviour is obtained as follows 


dH _ dH) (ἃ 
dT \dF/\4AT, 
47 41 is obtained from equation (14.6) as 
dF XRT 
47 AH? (To TY 
It is also assumed that, because of the restriction to consideration of ideal 
eutectic systems and to the formation of ideal solutions on melting, that 


Ἡ,Ξ- ΔΉ FE 
and therefore + 
ΔΠ ΡΞ ΛΠ Γι 
Combining these results 
dH! ERTS 
ape ι τ τὴς 14.9 
dT {0- μὰ 
Equation (14.9) then gives the variation of the heat capacity of the sample 
during melting as a function of 7. The upper limit of the melting process is 
T=T., (when F=1). Therefore, equation (14.7) becomes 
T= TR 
AH?, 


The lower limit of the melting process is 7.470. when dH/d7~xt R=C, ic. the 
heat capacity of the sample is approximately constant. In the idealized DSC 
curves given in Figs 14.1 and 14.2, it has been assumed that the heat capacity of 
the liquid just above the melting temperature is the same as that of the solid at 
lower temperatures (i.c., both equal to ὍΝ; 

Equation (14.9) can be written as 


aH χἱ R(To/T? 
dP {{π|1-ἰ|ὴ 
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during melting depends only on the mole fraction of impuri 
temperature relative to the melting point of the pure substance, 7/7). 

Since d/1/d1 is proportional to d/1/d7, plots of dH/dT against 7 represen 
initial part of an idealized DSC melting curve. Such curves for phenacetin 
for benzamide, with values of x from 0.0050 to 0.2000, have been give 
Marti et al. [1, 2]. 

The real DSC melting curve, because of factors such as thermal lag, whi 
discussed in more detail below, will look more like the curve in Fig. 14.2.1 
The total area under the curve, i.e. area ABC, is proportional to the enthal 
fusion, AH’, . The actual value of A//;", can be obtained by calibration ¢ 
instrument with a standard of known AM” (Fig. 14.1). The feature se 
for the present discussion, the fraction of the sample melted, Γ᾿ at temper: 
T, is obtained directly from the fractional area under the curve 
F=area ADE/area ABC. The range of F values used in practice is us 
restricted to 0.1 << 0.4, Even with this restricted range, the linearity of ph 
T against 1/F (equation (14.7)) is often poor. Corrections have to be mac 
thermal lag and for undetected premelting as discussed below, 
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14.4 Corrections 
14.4.1 Thermal lag 


Flow of thermal energy from the holder at the programmed temperature, 
the sample at a slightly lower temperature. T.. is governed by Newton's 

eee OH 

pte ge Re a 
where Ry is the thermal resistance, The value of R, for the instrument m 
obtained [3] from the melting curve of a high purity standard (Fig. 14.1) 
will melt over a very narrow temperature range. so that as the progran 
temperature continues to increase with time, (he sample temperature. 7, 
remain constant, i.e. d7,/dt=0. From equation (14.10) 
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Hence 


and 


R, may thus be determined from the slope, AB, of (he DSC. metung curve rr τὴν 
high purity standard (Fig. 14.1). The value of Ry is then used graphically or 
analytically to correct the programmed temperature, T,, to the true sample 
temperature, Τὶ. 7, then, rather than 77, is plotted against 1/F. 


14.4.2 Undetermined premelting 


Even with correction for thermal lag, the linearity of the plots of 7, against 1/Fis 
often not good, and corrections have to be made to the measured areas for 
melting which has occurred at lower temperatures and which is difficult or 
impossible to measure. This is evident in the small but cumulatively significant 
deviation from the ideal baseline illustrated in Fig, 14.2, 

Ifthe undetermined area under the curve is δ, the measured partial areas up to 


temperatures Τὶ. 7γ... ὁ ΤΊ ATE Ay, 424+ - + Ans respectively, and the measured 
total area is A, the true value of F, is 

F,=(a,+8)/(A +8) 
or 


1/F,=(A+e)/(a,+#)~A/(a, +), 

since A>e, so the effect of including ¢ is to reduce the value of 1/F, (Fig. 14.4). 

In practice eis treated as a parameter whose value is adjusted so that a plot of 
T, against the corrected 1/F is linear. The restraints are that the final value of 
(4: 6) should correspond directly to the correct value of AH, (if known) 
and that the value of Τὸν determined from the intercept on the 7, axis, 
should be correct. Once these conditions have been met, the value of 
xt (=slopex ΔΗ͂ ΓΙ /RT2) and hence the purity of the sample can be 
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Figure 14.4 Correction for undetermined premelting. The correction = 100% 2/A and 
F=(a+e)/(A +e). (With the permission of the Journal of Chemical Education.) 
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may be far below the range of temperatures being examined, the correctir 
may sometimes be quite large and values of as much as 30% of the total are 
not uncommon [9]. Obviously the approximation, Α Ἐπ ον 4, cannot th 
used, Sondack [5] has suggested an alternative procedure, 

The whole procedure for purity determination is summarized in the low 
(Fig. 14.5). Computer programs for these calculations have been describe 
2, 4]. More detail on the procedure is given in τε Ὁ; 


IDSC curve for impure) 
(x}) sample 


measure areas, an, of 
endotherm up to 
temperatures, (Tp)n 


calculate total area 
A= an 
correct (Tp)n to (Ts)n 


calculate Fy = an/A 


plot(Ts)n against 1/Fy 


correct for 
premelting. 
Add € to an 
andto A 


ῷ 


yes 


no 


does the value 
of A correspond to the 
correct value of ΔΗ ἡ 
(if known)? fe 


yes 


calculate x! from 
slope = xSAT?/AHF , 


DSC curve for pure 
reference 


determine thermal 
resistance,Ro 


measure total area of 
endotherm 


calculate 


calibration factor 
from known ΛΗ 


Figure 14.5 Flowchart representing the procedure used in purity determination 
[6]. (With the permission of the Journal of Chemical Education.) 
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Staub and Perron [7] have shown that a stepped heating technique can extend 
the working region for purity determination. The sample is heated through the 
melting region in steps of a few tenths of a degree, thus allowing a closer 
approach to true thermodynamic equilibrium. The results of such a procedure 
[8] on an impure phenacetin sample, using a modified DSC-1B are shown in 
Fig. 14.6. There are 6 steps/5 K and the areas of cach peak. obtained by 
integration. are given. The first 1! peaks and the last 3 arise from the difference 
in the heat capacity of the sample and reference. These peaks have 
approximately constant area and the melting peaks can be corrected for this 
difference. Because thermal equilibrium is established after ench peak, no 
correction is required for the thermal resistance, Ry, of the system. The 
corrected areas of all the melting steps are then summed and converted to 
fractional areas, F, and 1/F is plotted against 7, as before, except that with this 
procedure [8] all the data are used including the latter part of the melting 
process. The linearization process for the undetermined premelting still has to 
he carried out, 

Gray and Fyans [9] have suggested an alternative procedure to that given 
ahove, in which the mole fraction of impurity (x4) is calculated from the areas of 
two consecutive stepped peaks (x, and ας. ,). the magnitude of the stepping 
interval (AT). and the molar mass (M) and the melting point (74) of the pure 
major component. This method depends upon the applicability of the van't Hof 
equation and the relationship derived [9] is 


The melting point (7)) of the pure solvent may be determined [9] from the areas 
σα, and o, ,. the step interval AT and the final temperature of the step 7), 


It is recommended [9] that as large step increments as possible should he used. 


14.6 Conclusions 


Raskin [10]. ina critical review of methods of purity determination using DSC, 
concludes that the accuracy of the method is generally overestimated. A realistic 
estimate of the accuracy in the range 0.005 < τῇ <0.02 is given [10] as 30.50% 
when the sample mass is less than 3 mg and the scanning rate is less than 
2K min ! 
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Garn et al, [11 13] have discussed the problems that arise when there is 
appreciable solid solubility. They used NMR to detect the solidus and 
compared DSC and NMR results. They showed that lack of thermal 
equilibrium is not a principal source oferror in the method and also found that 
the measured impurity content is sometimes dependent upon the nature of the 
impurity as well as its concentration. 
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Chapter 15 


The literature and nomenclature ὁ 


thermal analysis 


A selection of references is given below to serve as an introduction tot 
literature of thermal analysis. A more comprehensive list is given in For 
Thermal Analysis, published by ICTA, and soon 10 come out ina new edi 
Dr J. O. Hill. 
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15.2 Reviews 
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Many other more specialized reviews have appeared in the Journals 
mentioned in section 15.4. 
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and the European Symposium on Thermal Analysis and Calorimetry (FS 
are valuable surveys of the ‘state-of-the-art’ (Tables 15.1 and 15 ray 

The annual conferences of the North American Thermal Analysis § 
(NATAS) are published in short form and are very useful. 

A series on Analytical Calorimetry based on the ACS Symposia, editec 
S. Porter and J, F. Johnson. is published by Plenum Press: Vol. 1. 1968; 
1970: Vol. 3. 1974, Vol. 4, 1978 and Vol. 5, 1984. 

The series Progress in Vacuum Microbalance Techniques, publist 
Heyden, and the earlier volumes Vacuum Microbalance Techniques. pul 
by Plenum Press, are sources of information on TG. 


Table 15.1 International Conferences on Thermal Analysis (ICTA) 


ICTA Date Place Editor Publisher 

8th 1985 Crechoslovakia A. Blavck Thermochim 
Vols 92, 93, 

7th 1982 Canada ἢ. Miller Wiley 

6th 1980 Germany ΗΠ G. Wiedemann Birkhauser 

Sth 1977 Japan HW. Chihara Heyden 

4th 1974 Hungary 1, Buzas Heyden 

3rd 1971 Switzerland H. ὦ. Wiedemann Birkhauser | 

2nd 1968 United States R.F. Schwenker Jr. Academic P 

and PD. Garn 
Ist 1965 United Kingdom J. P. Redfern Macmillan 


Table 15.2 European Symposia on Thermal Analysis and Calorimetry (ES 


ESTAC Date Place Editor Publisher 

3rd 1984 Switzerland B. Ε΄. Rordorf Thermochin 
Acta Val & 
1-527 

2nd 1981 United Kingdom Ὁ. Pollimore Heyden 

Ist 1976 United Kingdom D. Dotlimore Heyden 


15.4 Journals 


The Journal of Thermal Analysis (Wiley) and Thermochimica Acta (Flse 
the main specialist journals, although results of thermal analyses a 
many other journals, especially Analytical Chemistry, Talanta, ¢ 


Chimica Acta, International Laboratory, Laboratory Practice, Analyst, and the 
many polymer journals. The Journal of Thermal Analysis has ἃς useful 
bibliography section and Wiley also publishes Thermal Analysis Abstracts 
(Fig. 15.1), 

Computer-based searching services, ¢.g.. CA Selects: Thermal Analysis, 
published by the American Chemical Society, are a convenient way of keeping 
abreast of the literature (Fig. 15.2). 


15.5 Manufacturer’s literature 


The major manufacturers (listed in chapter 16) provide extensive information 
on thermal analysis. Special mention must be made of the Thermal Analysis 
Newsletters. Thermal Analysis Application Studies, Instrument News and 
bibliographies produced by Perkin-Elmer; the Applications Briefs and Du Pont 
Thermogram; The Cahnsultant; The Technical Information Sheets from Stanton 
Redcroft. and the Thermal Techniques Series from Mettler. The manufacturers 
often provide reprints of articles, published in the literature, but based on their 
equipment, 


15.6 The nomenclature of thermal analysis 


The International Confederation for Thermal Analysis (ICTA) has published 
several recommendations for the standardizing and reporting of results of 
thermal analysis. Many of these are contained in a booklet For Better Thermal 
Analysis published by ICTA, Other references are: Mackenzie, R. C. (1969) 
Talanta, 16, 1227; (1974) Pure Appl. Chem., 37, 439; (1972) Talanta, 19, 1079; 
(1975) J. Therm, Anal., 8; (1979) Thermochim, Acta., 28, 1 and (1981) 46, 333. 
See also The Metrication of Thermal Analysis or Conversion to 4] Units, Blaine, 
R. L. (1978) Thermochim. Acta, 26, 217-228. 


LILCTEMIGL Διιαιμσιῦ ot 
An official organ of the International Confederation for Thermal Analysis 


Sample Abstract 


Abstract Number ΜΙΝ) 
Title TEMPERATURE DEPENDENCE OF THE HEAT CAPACITY 
CHANGE FOR THE DISSOCIATION OF ACETIC ACID ANID 
PROTIONIG ACH IN WATER 


Author(s) Olotssan, G 
Address Thermochemistry Laboratory, Ghemieal Contre Ὁ Row? 
S.22007 Lund, 5! n 
Reference | Chem. Thermodyn . 1611). 3044. 1984 (Fre) 
(Language) 
Abatract The malin enthalpy of praton dissociation, determined t 
Various temperatures hatween 274 and 378 bh. varted 
Hneary with femperature for both acide Phe two acids 
yave virtually (he same heat capacity. the average vat 
OE aCp am bein 142 1K mol 


Keyword It) 109, 151, 1ΜῈ 220, 227 την, 995 
Numbers Methodology Trennique, Materials ΠΣ 
Application Subjects Compaund 


Abstractor 0S 


KEY 10 -- the method used was calorimetry 
1oo = determination of heat of raaction, ety 


131 = determination of specific heat ar heat canary 
136 — determination of thermodynamic constants 
27) — the materials investigated were liquid 


227 — the materinis investigated were organic 

+19 — the matoriate investigated Involved organic acid 
or salts 

505 — the materlals Investigated Involved water 


Figure 15.1 Sample entry from Thermal Analysis Abstracts, (With the permis: 
John Wiley & Sons, Chichester.) 


ITEM NO 013127 CA VOL 105 NU O02 * PRUFILE 06 


AUTHORS: WENOLANODT. We welled 
(4) DEP. CHEM«s UNIV. HOUSTON. HOUSTONe TX,» 77004% 


USA 


TITLE: THE DEVELOPMENT OF THERMAL ANALYSIS ([NSTRUMENTATION 
1955-1985 


REFERENCE: THERMOCHEIMs ACTA (THACAD vou τοῦ NO Lb PP 1-22 
(1986) ISSN 00406031 


INOEX TERMS: 
SECT LUNS: CA069000 THERMODYNAMICS» THERMOCHEMIST AY 


CA020; CAOTI; CAOTIS caoao 


DESCRIPTORS: 
THERMAL ANALYSIS¢ APPecee (OEVELOPMENT OF) 


IDENTIFIERS: 
REVIEM THERMAL ANALYSIS INSTRUMENT» OTA LNSTRUMENT πενι 


THERMOELEC ΤΗΌΜΕΤΗΥ INSTRUMENT REVIEWs LIGHT EMISSION 
INSTRUMENT THERMAL ANALYSES REVIEW. THERMUBALANCE 
INSTRUMENT THERMAL ANALYSIS REVIEW 


SEARCH TERMS PRESENT: THERMAL ANAL 3 CAOT9 
NO ABSTRACT GIVEN 


Figure 15.2 Sample output from a computer-based literature searching service ct 
Chemical Abstracts. (National Institute for Informatics, CSIR, Pretoria ) 


Saat tees a Aaa αν 


Chapter 16 


Thermal analysis equipment 


16.1 Choosing thermal analysis equipment 


A bewildering array of equipment is available to anyone starting out in thermal 
analysis and selecting a suitable system is not an easy decision. Tt is advisable 
first to collect the pamphlets and specifications of most systems from the 
suppliers (section 16.2), and then attempt to make a preliminary selection by 
considering the following factors. 

(1) What types of sample are going to be examined both immediately and as far 
as can be predicted, in the future? " 

(2) What sort of information on the sample is required, e.g. thermal stability, 
glass-transitions, percentage crystallinity, mechanical properties, details of 
gases evolved etc.? 

(3) Over what temperature ranges are the changes that are of interest likely to 
occur? 

Answers to the above questions can be used to eliminate definitely unsuitable : 

systems. Unless there are very special requirements such as extreme ranges of 

temperature, or use of very corrosive atmospheres or strongly exothermic or 
even explosive samples, there will usually still be a wide choice of instruments. 

Prices related to budget available will obviously be a further restriction, and 

modular systems, which can be added to, are attractive. Most systems 

incorporate their own computers and it is worth looking carefully at this 

(chapter 12) to see whether a system can be interfaced to existing laboratory 

computers. This should only be contemplated if you have a source of software. 

Some systems cannot operate without the manufacturer's computer. 

A most important question to be answered before choosing a system, is the 
training and service available from the suppliers (section 16.2). This involves 
seeking out other users and checking on their experiences. At the same time it is 
worth checking on the prices of spares and accessories (which can be 
ridiculously high), Even the best equipment (which need not be the most 
expensive) will be difficult to operate and maintain if the local agents cannot 
provide informed and rapid service. 


Manufacturer 


Cahn-Ventron Corporation 


DuPont Company 
Instrument Systems 


Harrop 


Linseis Messgerate GmbH 


Maple Instruments 


Mettler Instrumente AG 


Netzch-Geratebau GmbH 


Perkin-Elmer Corporation 
Analytical Instruments 


Polymer Laboratories 
Rigaku-Denki Co. 


Setaram 


Shimadzu Corporation 


Stanton Redcroft 


pana ea eee 


Address 
16207 South Carmenita Road, Cerritos, CA " 
USA. 


Concord Plaza, Quillen Building, Wilmingto 
19898, USA. 


Harrop Laboratories, 3470 East Fifth Avenu 
Columbus, Ohio 43219, USA, 


Vielitverstrasse 43, 8672 Selb, FRG 


Pissummerweg 1. 6114 AH Susteren, The 
Netherlands. 


CH.-8606 Greifensee, Switzerland, 
P.O. Box 71, Hightstown, NJ 08520, USA, 
P.O. Box 1460, D-8672 Selb, FRG 


Main Avenue (MS-12), Norwalk, Conneeticr 
06856, USA. 


Fssex Road, Church Stretton, SY6 6AX, Eny 


9-8, 2-Chome, Sotokanda, Chiyoda-Ku, Tok 
101, Japan. 


7, rue de l'Oratoire, 69300 Caluire, France. 


Shinjuku Mitsui Bldg, 1-1, Nishi-Shinjuku 
2-Chome, Shinjuku-ku, Tokyo 160, Japan. 


Copper Mill Lane. London SW17 OBN, Eng 


Chapter 17 


Conclusion 


In the preceding sections, only a few selected examples have been given of the 
vast number of applications of thermal analysis. To get a better idea of the range 
of applications over all branches of inorganic, organic, physical and industrial 
chemistry. metallurgy. polymer science, glass and ceramic science, food science, 
biology, geology, pharmacy, medicine, agriculture and engineering, it is worth 
scanning the contents pages of some of the proceedings of the International 
Conferences on Thermal Analysis (ICTA), as well as the biennial reviews by W. 
W. Wendlandt in Analytical Chemistry (detailed references are given in 
chapter 15), Although applications have been grouped under the main 
techniques used, many of the studies reported have been, or eould have been, 
extended by the use of simultaneous and/or complementary measurements 
(chapter 9). 

As suggested in the opening paragraphs of this book, there are very few 
materials which will not show interesting changes on heating. The materials 
studied haye ranged from kidney stones to synthetic diamonds, and from 
ancient papyri to the latest polymers and ceramics for space research. The 
information obtained has been used to solve old problems and to develop new 
processes for the future. 


Appendix A 


Introductory experiments in therme 


analysis 


ven below. What can be do 


A selection of introductory experiments is gi 
available, and constant ref 


obviously be determined by the apparatus 
should be made to the manuals supplied with the instruments. Μ 
computerized instruments will usually operate interactively, providing s 
the directions given in the procedures below. 


A.1 Differential scanning calorimetry (DSC) 


A.1.1 Calibration 

Calibrate the DSC with respect to temperature and heat flow. Check 
reproducibility. Compare the results obtained on different instrum 
possible. 


A.1.2 Determination of peak areas 


Use a planimeter to determine areas recorded on chart paper. Use cc 
programs to carry out numerical integration. 


A.1.3 Dehydration 


Determine the temperatures anden 
CuSO, .5H,0. 
Carry out similar measurements on some other hydrates. e.g. Bat 1 


Ni(HCOO), .2H,0 or Ni(COO),.2H,0 and. see how the cnt! 


dehydration per mole of H,0 varies from salt to salt. 


thalpy changes for the dehydration 


A.1.4 Decompositions 


ἢ Study the decompositions οἵ οὁπις metal carboxylates, ¢.g.. 
: Ni(H1COO), .2H,O or Ni(COO), .2H,0, after first having dehydrated them 
(sce section A.1.3). ‘Try to carry out an isothermal DSC run at a suitable 
temperature, Carry out kinetic measurements 


A.1.5 Phase transitions 


=e 
EES SS 


Determine the temperatures and enthalpy changes of phase transitions in salts 
such as NH,NO,, KNO4, KCIO, or Ag,SO,,. Study the reversibility of these 
changes on cooling and comment on the tse of these transition temperatures as 
temperature standards for instrument calibration. Use a hot-stage microscope 
(FISM) for visual detection of the phase changes. 


= a) 
“+ 
a 


4 
ryt ν a 
ie ἢ A.1.6 Specific heat capacity 
LA { 
Ha ἡ Determine the specific heat capacity of an inert substance, such as Al,O,, 
pe ᾿ relative to that οἵ aluminium metal (0.900 J K ~'g~! at 25°C). Comment on the 
ri at variation of the specific heat capacity with temperature. 
eo 


A.1.7 Purity determination 


Compare the melting endotherm of pure indium with that for benzoic acid and 
estimate the purity of the benzoic acid. 


A.1.8 Polymer stability 


- 


ast A 


Do runs on some polymer samples in nitrogen and in oxygen. Examine the 
results to determine the glass-transition temperatures and the temperatures at 
Ll which melting. degradation and oxidation occur. 


A.1.9 Polymer crystallinity 


Polyethylene (PF) is a semicrystalline thermoplastic. From a DSC curve, 
determine the temperature range over which melting occurs as well as the 
enthalpy of melting. The percentage crystallinity is calculated by comparing the 
measured value with that for 100% crystalline material (290 J g Lith 


A.1.10 Curing of an epoxy resin 


Do a DSC run on epoxy resin. Determine the enthalpy change and the kinetics 
of the exothermic curing process. Rescan the product and get the glass- 
transition temperature of the polymer. 


A.2 Thermogravimetry (ly) 


A.2.1 Temperature calibration 


Use magnetic standards of known Curie point, to calibrate the fu 
temperature. 


A.2.2 Dehydration 
Determine the temperatures and mass losses accompanying the dehydi 


stages of CuSO, .5H,0. Compare your results with those from the 
experiment A.1.3 and draw up a full description of the dehydration pro 


A.2.3 Decompositions 


The studies carried out in the DSC experiment A.1.4 can be complemen 
determining the mass losses. Kinetics may he determined from a 
isothermal run or from a series of isothermal runs. 


Α.2.4 Percentage filler in a polymer 


Decompose a sample of an epoxy putty in nitrogen. The residue is the 


A.2.5 Analysis of coal or of a rubber 


(a) Heat a sample in nitrogen until no further mass loss occurs. This gi 
proportion of volatile material 

(0) Change the purge gas to oxygen while holding the sample at the high 
the temperature range, The mass loss corresponds to the proper 


carbon residue. 
(c) The mass of the residue in oxygen corresponds to the inorganic as 


Appendix B 


Thermal analysis software 


Listings of several programs, written in Applesoft BASIC by the author, for use 
on APPLE ΤΙ PLUS microcomputers with at least 48K, are given below. 
Although every care has been taken, no responsibility is assumed for the correct 
operation of these programs. Copies of the programs may be obtained by 
sending a blank diskette to the author. The programs should in most cases be 
fairly readily adaptable to other microcomputers. δ 

The sequence of operations in the capture and processing of data from 
thermal analysis experiments is usually: 
Β.1 Data capture and temporary storage, 
B.2 Permanent storage; 
B.3 Display of data for preliminary examination; : 
B.4 Modification of the data; (a) baseline correction; (b) smoothing, (c) scaling; , 
B.S Processing of the data; (a) numerical differentiation; (b) peak integration; 
B.6 Other calculations, e.g., (a) kinetic analysis; (b) purity determination, etc. 


B.1 and 2 Data capture and storage 


Program DSC is given here as an example. This program has been designed for 
operation with a Perkin-Elmer DSC2 but is readily adaptable to other 
instruments and to other TA techniques. The program given here is based on 
the use ofa Datel Systems Inc., ADC-ET12 12-bit analogue-to-digital converter 
(ADC) (sce program lines 315 to 330) and a Mountain Hardware clockhoard 
(sec lines 175 (0 314). Full details of the hardware are available from the author. 


Input: ADC channel number, filename for data storage, starting and 
stopping temperatures, heating rate and number of points to be 
collected. 

Output: Temperature and DSC response stored in a sequential text file. Data 
is plotted on the screen and may be dumped to a printer (Epson). 


REM CHECK: THAT FROGRAM 16 LOADEY Rruve 


IF PEEK (104) = 96 THEN 5 
POKE 103,12 FOKE 104,94: FORE 4096 * 6,9 
FRINT CHR (4)"Rl JN DSC" 


HOME 
PRINT (aE ρυ  μἕἔἐἘἔμἔμ ρον I 
PRINT "* FROGRAM DSC 29/05/87 MER *" 


PRINT RENEE IESE IG EOE EER a EEER ISITE 
PRINT + PRINT 

AD = S:SLOT = 2: REM AD =A/D SLOT NO., AND SLOT 
Lock SsloT NO. 

PRINT "PROGRAM TO CAPTURE DATA FROM THE PERE TN 
MER DSC2 ": PRINT 

FRINT "THIS VERSION USES THE MEDIUM SPEED 12 BI 
Q/D ROARD SUPPLIED HY RHODES ELECTRONICS AND FL 
IN δι τεῆς" SET IN LINE 15, AND THE MOUNTAIN 
Lock IN SLOTH":SLaT:;" ALSO SET IN LINE 15 “= FF 


PRINT 
REM oceania ea | δὶ 
INVERSE 

INPUT "ENTER A/D CHANNEL NO. ὦ TO S Pe 
we< OORC DS THEN τὸ 

PRINT +: INVERSE : INFUT "DO YOU WANT TO STORE 
DATA IN A DISK FILE 7 Y OR N ΕΣ “sks 

4 


IF κα > "Ἢ" THEN 39 

PRINT τ INFUT "ENTER FILENAME "3G: PRINT 

HOME 

PRINT : INPUT "ENTER STARTING TEMP/KELVIN "sKO 
IF KO < 300 THEN FRINT "TOO Low ! "τε GOTO 40 
PRINT +: INPUT “ENTER STOPPING TEMP/KELVIN "skKF 
IF KF > 750 THEN PRINT "TOO HIGH ! "κα GOTO 42 
PRINT : INPUT "ENTER HEATING RATE/K FER MIN, E 


R O FOR ISOTHERMAL "sHR : 
IF HR < © OR HR > 320 THEN PRINT "OUTSIDE SET 
GS AVAILABLE "+ GoTo 44 
ΤΕ HR < > 9 THEN 50 
PRINT : [NFUT "ENTER TIME INTERVAL BETWEEN 1st 
RMAL READINGS, IN SECS "3H 
PRINT : INPUT "ENTER NO. OF POINTS TO ΒΕ COLLEC 
N.B. MAXIMUM IS 500 "Wz: PRINT 
IF W > 500 THEN So 
ἐκ < 3 ye THEN 64 
INPUT "IS THE RIGHT DISK IN DRIVE 1 FOR DATA 
AGE 2 Y OR N τ weve 
IF Y® < >» "Y" THEN 62 
PRINT : TF HR = © THEN 68 
H = 60 κ (KF - Ko) ἡ CHR * W) 
REM H = TIME INTERVAL = 60 X T RANGE / ( HEA 
RATE X NO, OF FOINTS ) 


6A IF H < 0.5 THEN PRINT "TIME INTERVAL BETWEEN σοι. 
LECTION OF FOINTS 15 TOO SHORT. REENTER ": GOTO 4 
7 
70 DIM R(W + S),VCW + 5) 
71 DIM TiW + 5) 
90 Dt = CHR (4): REM "CONTROL ἢ" 
100 PRINT D#;"NOMONI,O,C": REM -KEEP DISK COMMANDS F 
ROM PRINTING 
110 HOME : REM -CLEAR SCREEN 
170 1 = Of REM - FLAG FOR START AND INITIALISE COUNT 
FR 
29 NORMAL 
πὸ PRINT "POINT TIME/S", "TEMF/K", "DSC OUTF" 


= 


᾿σὸ ΡΓΠΙΝΥ eee aeration on ΞΟ ΌΡΝΜΕ ἢν 
IF X # 1 THEN 510 
REM -= X IS USED AS THE FLAG FOR FRINTOUT. IT IS 


SET IN LINE S00 
159 FRINT : FLASH 


ae. 
a 


Wye, 160 INPUT “READY TO START 7 Y OR N 7 “:85 

hee i 141 NORMAL 

CW 170 IF δὲ < > "Y¥" THEN 160 

ἡ ἃ 175 REM ς 

EC, 180 PRINT D#;"IN#"sSLOT: REM -SETS INPUT TO CLOCK 
Ya 190 PRINT D#:"PR#";SLOT: REM -SETS OUTPUT TO THE CL. 
len ofr 
Vent | 200 INPUT " ";T#: REM -OBTAIN THE TIME 

᾿ \ 202 Ht = MIDS (T#,7,2) 

ear 203 REM H# CONTAINS THE HOUR PART OF Τῇ 


210 Mf = MID* (T#,10,2) 

REM Μὲ CONTAINS THE "MINUTES" FART OF Tt 

Sf = MID® (T#,13,2) 

REM 584 CONTAINS THE INTEGER "SECONDS" FART OF T 
£ 

230 F# = RIGHT# (T#,3) 

O31 REM πᾷ CONTAINS THE FRACTIONAL PART OF "SECONDS 
* FROM TS 

240 St — St το" 4+ FS 

341 REM INTEGER AND FRACTIONAL FARTS OF "SFCONDS" A 
RE COMBINED WITH A DECIMAL POINT 

250 T = 7400 * VAL. (HE) + 60 * VAL (Mf) + VAL (S#) 


0.06 


ye 
it 
teats 
NNN 
Ν᾿ 
- oe 


| 
“st 
ae 


‘£3 
f 


προτάξας τ᾽ 
πόδι ας 
{γι ) 
ae στο τες ἈΞ, ines Sa ae ee ᾿ξ ναῦς 
N » 
Ὁ a 
x “ἘΦ 


REM CONVERTS HOURS AND MINS ~-*SECONDS AND AD 

i DS SECONDS 
it IF 1 <¢ > © THEN 280 
AS 261 REM -- CHECKS FLAG FOR START OR MIDDLE OF RUN 

yon 270 10 = Tr = TF 

‘Ox | 271 REM - STORES STARTING TIME 

NG 272 ΤΙ) = Os REM INTTTALISE RUN TIME 

Ua 200 Rk = T - TO 
We 2861 REM - MEASURES ELAPSED TIME 


290 %F J = 0 THEN 310 


᾿ς 


πες τὸ 
= 


“Ort 


4 


ec 


ἘΞ ΤΕ πα Ὁ aE ee 
REM — CHECK WHETHER TIME INTERVAL HAS El APSE 
N.B.H=TIME INTERVAL IN SECONDS CALCULATED IN Ut. 
NE 65 


THE), ὦ ΤΟΙΣ 4 
R(I) = (R κ HR / 60) τ KO 
R(I) = INT (1000 * R(1)) 7 1000: REM 3 DEC PL 
CES 
PRINT D#;"IN#O" 
PRINT D#;"FRHO" 


POKE - 16256 + AD * 16,C » 721 POKE ~ 1224 1 
AD * 16,12 

FOR N = 1 TO 20: NEXT N 
vit) = ( PEEK (© - 16256 + AD # 16) # 256 ὁ ΡΈΕΙ 


( = 16255 + AD * 16)) - (Cc κα 8192): REM RFAD ci 
ννει. C 
TOD) = INT (TCE) ἢ 1000) / 1000:R(T) = INT (R 
» * 1000) ἡ yooo:V¢r) = INT (VET) κα 1000) 7 190 
PRINT T:" “ΤΟΙ ROL) VCD)? PRINT 
ye Pe A IREM, - COUNTER 
IF I = W + 1 THEN Ao0: REM N.B. WeNO. OF POTN 
TO RE COLLECTED, SET IN LINE 50 
GOTO 180 
PRINT : PRINT "STARTING TIME = “ΤΟ, 0 
PRINT : PRINT "STOPFED AT T = “shy ts 
REM Ξ τ ΞρΕορορορρποο OUR a ee ee 
INVERSE 
PRINT + INPUT "DO YOU WANT A PRINTOUT Ὁ Y OR N 
"SAE 
IF Af = "N" THEN 510 
PRINT = INPUT "IS PRINTER ON 7 Y OR N 7 "SA 
NORMAL. 
1Ὲ As < > "Y" THEN 4so 
PRY L 
PRINT : PRINT "OUTPUT FROM DSC " 
IF Ks << > "Y" THEN 500 
PRINT "STORED IN FILE "sG¢ 
PRINT "STARTING TEMF/K = "skO: PRINT "STOFFING 


EMP/K = "s3KF: PRINT “HEATING RATE 7 K PER MIN 
Re PRINT 

PRINT " ποτα τα τε τα τα τα τὸ τα το τα i ee ee Μὴ δα μα Βα ΜῈ μὰ τα το γον μτϑα σον 
PRINT :X = 1: GOTO 130 

REM PUTTING X~1 SETS FLAG FOR FRINTOUT, SEE 
ES 150 AND 420 
MAX = O:B = V(O) 

FOR I = 0 TO W 

IF At = "N" THEN 550 


ΞΟ 3. 


sy 
ON)" 


es 


ome. 
yO) 


== 
Ly ans ἢ 
a 


‘ ¥ 
Sadia cs 


Ss ea 


a3: 


2 


as 
a ΓΣ 


i 

ΝΥ 

ΩΝ 
he 


sao 
sso 
Boo 
$70 
SAG 
B90 
“ΩΦ 
600 
619 
620 
630 


65 
ὁ. 1 


652 
654 
658 
659 
640 


bbl 
6709 
675 
680 
685 
690 
7oOO 
720 
730 
740 
74S 
750 
760 
770 
780 
790 
BOO 
810 
820 
B20 
B40 
845 


850 
8460 
B45 
870 


880 
B90 


PRINT ¢ FRINT 13" 
ΤῸ MAX THEN 570 
MAX = V(T) 

If Vl) » BR THEN 590 


"eT(T) RIT) VET) 


Roe V(T) 

NEXT I 

NORMAL. 

PRINT : PRINT "MAX = “sMAX,"MIN = "SR 

IF Af =~ "N" THEN 652 

PRINT + PRINT "DATA FROM CHANNEL "τὸ 

PRINT + PRINT "READ AT TIME INTERVALS OF “ety” 8 
FCS": PRINT 

PRY O 


PRINT : FLASH + PRINT “WHEN ERINTER HAS FINISHED 
PRESS ANY KEY TO CONTINUE ": NORMAL : GET VF 

IF K€ ¢ > "¥" THEN 659 

ROTO 1590 
REM -----τ- τ - - Ριιοτ----- -΄------------ 
1NVERSE 
PRINT τ INPUT "DO YOU WANT A PLOT OF YOUR READIN 
as 7 Y OR N 7 "sSAt 

NORMAL 

IF AF ¢ > "VY" THEN 1990 

If HR = © THEN D1 = Τί): GOTO 685 
Di = KF - KO - 
RANGE = MAX - B 
γ9 = MAX + 0.05 # RANGE: YO = B = 0.05 * RANGE 
v2 = Y9 -- YO 


ΙΕ D2 2 "5 Ὁ THEN 940 
REM ~-QUTLINE FRAME 
HGR 

HCOLOR= 3 


HPLOT 9,0 TO 279,90 TO 279,150 TO 9,150 TO 9,0 
REM ~TIC INTERVALS 
G4 = INT (270 / 10) 
Ω5 = INT (150 / 10) 

FORK = 1 7170 9 

HPLOT 9 + kK * 64,150 TO 9 + K * 64,147 

HPLOT 9,147 - K * GS TO 12,147 - Κακὰ G5 

NEXT K 

REM -PLOTTING OF THE POINTS 

For l= Oo TOW 

IF HR = © THEN X = INT (270 * ΤῸ / Di) + 9: GOTO 


840 
y= INT (270 * (R(T) -- KO) / pi) + 9 
Vom 1S0 = INT (50 Ἢ CV Clore yo) / D2) 


{Fox 4) OF OR % $1279 OR Y <0 OR Y > 150 THEN 900 


IF re OORV= OORU® © THEN HPLOT X,Y: GOTO 
890 

HPLOT V,U TO X,Y 
Ym XU) ¥ 


910 PRINT 

920 PRINT "X TIC = "spt / fog" ¥Y TIC = Me B22) 10 

930 GOTO 1000 

940 PRINT "Y MAX ~ Y MIN 15 ZERO FLOT ABANDONED" 

950 GOTO 2000 

1000 GOSUB 7000 

1010 ΣΕ ze < > "¥" THEN 2000 

1141 PR# 1 

1142 PRINT : PRINT "DATA FOR GRAFH": IF K* = DO ae 
PRINT "STORED AS FILE "368: PRINT 

1150 IF HR < 5 Ὁ THEN 1243 

1160 PRINT + PRINT "ORIGIN AT ZERO TIME AND Y = "3 
: PRINT 

4170 PRINT "TIME INTERVAL IS "sT(W)s" SECS ": FRIN 


1180 GOTO 1245 

1243 PRINT : PRINT “ORIGIN AT T = "skOs" AND Y = 
O: PRINT 

1244 PRINT "T RANGE IS “εκὸς 


"TO "SKF," koe FRI 


1245 PRINT "Y RANGE 1S “s¥Os" TO “sV95" ON 12 RIT 
ALE FROM © TO 4095" 

1246 IF HR = © THEN FRINT "TIC INTERVAL ON TIMF αὶ 
8 IS "sT(W) 7 105" SECS ": GOTO 1248 

1247 PRINT "TIC INTERVAL ON T AXIS IS ὙΨ091 7 103” 
᾿ς PRINT "TIC INTERVAL ON Y AXIS IS "402 / 10 

1248 FRINT 

1249 FRINT “HEATING RATE WAS "SHRs" K/MIN"? PRINT 

1250 FRINT : PRINT "ACTUAL MAX Y = “ἜΜΑΧ 

1300 PRINT "ACTUAL MIN Y = "sR: PRINT 

1400 PRO 

1580 GOTO 2000 

S65 REM τ ΞΜ WRITE: Τὸ DISKs=—--—">=~— 

1590 PRINT ¢ FLASH : FRINT "PLEASE WAIT - WRITING 
DISK ": NORMAL 

1600 PRINT D#; "OPEN": Gf 

1610 PRINT D#: WRITE": G¥ 

1620 FOR I= 0 TOW 

1622 IF HR = 0 THEN 1628 

1624 PRINT R(T) 

1626 GOTO 1430 

1628 PRINT ΤΩΙ) 

1630 PRINT Ν(1) 

1640 NEXT I 

14645 FRINT 9999: PRINT 9999 

1650 PRINT D#: "CLOSE" 

1660 GOTO 659 

1990 PRINT 

2000 FLASH : PRINT "RUN ENDED “= NORMAL 

2010 END 


pasos ae 


Se 2S 


A299 


ζῶ 


TOW 
TOD 
1 
FORO 
ΖΛ 


Γ' 


FO50 
760 
7FOI7O 
7ORO 
FORO 


REM - 


DUMF SUBROUT INE------~-> > > 


PRINT "IT TS POSSIBLE TO GET A PRINTOUT OF YOUR 


GRAPH " 

FOR T =~ 1 τὸ 30008 
PRINT : PRINT "DO 
Oh Vi ΒΝ ae 


NEXT T 


YOU WANT TO ATTEMPT A PRINTOU 


INPUT Z7#: IF Z# < > "Y" THEN 7210 

PRINT + INVERSE τ INPUT "ENTER PRINTER TYPE τ: Ε 
SOM = E , C-ITOH = C "3YF 

NORMA 

It Yt: τ “E" AND Y€ Φ > "C" THEN 7040 


IF Yt = "E" THEN 7190 


Df = CHR (4):E# = 


CHR# (27) 


INVERSE τ FRINT "REMEMRER APPLE FARALLEL TNTERF 
ACE IN SLOT 1 ": NORMAL 


aa is 


7100 FOR I = 1 TO 1000: NEXT I 

7110 PRINT Df: "PR#I": PRINT E¢; "FP" REM FROFORTIONA 
L. PRINT 

7120 PRINT D#;"FPRHO" 

7130 PRINT D#;"RLOAD CITOH.ORJ" 

7140 CALL 6400 

7150 PRINT D#;"PRA1" 

7160 PRINT E#;"N": REM RESET TO PICA 

7170 PRINT D#_"PRHO" 

7180 GOTO 7710 

7190 FR 1 

7200 PRINT CHR# (17) 

7210 RETURN 


B.3 Display of data 


Although the program DSC provides for display of the data captured and 
stored, it is often convenient to have a separate program for displaying the 
contents of the data files which soon accumulate. Decisions can then be taken 


on further 


procedures to be followed. 


Program FASTDISP: This is a general-purpose program for displaying the 
contents X, Y of an APPLE sequential file. 


Input: 
Output: 


Name of the file containing the data. 

Display on screen ofa plot of X from XMIN to XMAX against Y on 
a scale extended 10% above and below the ὙΜΙ͂Ν and YMAX 
values. The points are joined (see line 1080) and the display is 
automatically dumped to an Epson printer. 


Program READISP: This isa more versatile, but slower, version of FASTDISP 
and includes option for deletion and addition of points, 


1 REM EQUUS «VERE MST τ 


2 IF PEEK (104) = 96 THEN 10 
πὶ POKE 103,12 POKE 104,96: ΡΌΚΕ A096 * 4,0 
4 PRINT CHRE (4) "RUN FASTDISF" 


10 HOME 

12 PRINT GORE IIE III SII RE” 
1523 PRINT "* *" 
14 PRINT "* PROGRAM 'FASTDISF’ *" 
15 PRINT "#*FOR QUICK DISPLAY OF DATA FROM A FILE®" 
18 PRINT "* " 


19 PRINT Το EIEEEA III IER? 

πο PRINT - PRINT "VERSION AS AT 1/08/86": PRINT 

50 PRINT + PRINT "Ν.Β. PRINTER MUST BE AN EPSON AN 
MUST BE ON ": FRINT - PRINT "POINTS WILL RE aor 
D SEE LINE 1080" 

60 D# = CHR (4) 

70 DIM A(5O2) ,R(SO2) 

BO PRINT =: INPUT “ENTER FILENAME "3G 

90 PRINT D#; "OPEN" GF 

100 PRINT D#;"READ":G6F 

105 PRINT + INVERSE + PRINT "FLEASE WAIT - READING 
NORMAL : FRINT 

110 I = 1 

120 INPUT ACI) ,B(T) 

140 IF ACI) = 9999 GOTO 170 

i50 I= 1+ 1 

160 GOTO 120 

170 PRINT D#;"CLOSE” 

180 PRINT "END OF ELE 

190 FRINT : PRINT “OUTFUT TO FRINTER NEXT" 

640 XMAX = A(1)2MX = ACI) 

670 YMAX = B(i):MY = B(1) 

680 FOR J = 2 TO 1 

690 IF οὖ) = 9999 GOTO 746 

700 IF ACJ) > XMAX THEN XMAX = ACI) 

710 IF A(J) « MX THEN MX = ACJ) 

720 IF B(J) > YMAX THEN YMAX = FCI) 

730 IF B(J) < MY THEN MY = B(J) 

740 NEXT J 


746 PR 1 
749 PRINT : PRINT “FILENAME = "GF 
750 PRINT : PRINT "XYMAX = ",XMAX 


760 PRINT "XMIN = "Mx 

770 PRINT "YMAX = ",YMAX 

780 PRINT "ὙΜΙ͂Ν = ",MY 

785 FRAO 

1030 X9 = XMAX:XO = Mx 

1035 RANGE = YMAX - MY 

1040 ὙΦ = YMAX + 0.05 κα RANGE: YO = MY - 9.05 # RAT 
1050 W = I: REM NO. OF FOINTS 


182 


1960 
1o70 
1080 
1990 
1100 
1110 
1120 
1120 
1140 
1142 
1144 
1146 
1148 
1150 
1152 
1160 
11980 
1190 
1210 
1220 
1230 
1240 
8 
1250 


1260 
1270 
1280 
1200 
1540 
1550 
1560 
1570 
1580 
1590 
0 
1400 
1610 
1612 
1614 
14620 
1620 


ETUC TTTUNAS weeee 


Di = X? - XO 
p2 "νῷ = YO 
24 = "¥": REM POINTS JOINED 
NORMAL 
REM -OUTLINE FRAME 
HGR 
HCOLOR= 3 
HPLOT 9,0 TO 279,90 ΤΌ 279,150 TO 7,150 TO 9,0 
REM TIC INTERVALS 
G4 = INT (270 / 10) 
GS = INT (150 ἡ 10) 
FoR K = 1 T0 9 
HPLOT 9 + K * 64,150 TO 7 + kK * 64,147 
HPLOT 9,147 ~ K * GS To 12,147 -Κὶ * GS 
NEXT Κὶ 
FoR t= 1 TOW 
A= ACT) 
R= ACI) 
REM -PLOTTING OF THE POINTS 
xX = INT (270 # (A - XO} ἡ) δ᾽ + 9: 
yer 150 - INT (150 * (Β - YO) / D2) 
IF X <4 O OR X > 279 OR Y < O OR Υ » 150 THEN 12 
fa) 
IF t= 10RV=0 OoRU=0O OR J#< > PVG: ΤΗΕΝ 
HPLOT X,Y: GOTO 1270 
HPLOT V,U TO X,Y 
v= X2U = Y 
NEXT I 
PRINT 
PRW 1 
PRINT CHR (17) 
PRY O 
PR 1 
PRINT ¢ PRINT "DATA FOR GRAPH STORED AS "3G6¢ 
FRINT = PRINT “ORIGIN AT X = "χοὸς" AND Y = "3Y 


PRINT + PRINT "X RANGE IS ";XO;" TO "sx? 
PRINT : PRINT "Y RANGE 1S. "YOu" Ta) bE AE 
PRINT + PRINT "X TIC INTERVAL = "401 / 10 
PRINT : PRINT "Y TIC INTERVAL = ":02 / 10 
PRE O 

END 


1 REM CHECK THAT PROGRAM IS LOADED ABOVE SCREEN# 

2 IF PEEK (104) = 96 THEN 10 

3 ΡΌΚΕ 103,1: ΡΌΚΕ 104,96: FOKE 4096 * 6,0 

A PRINT CHR# (4) "RUN READISF" 

10 HOME 

12 PRINT [GROSSE IDRIS IIIS EET EEO 

13 PRINT "* FROGRAM 'READISP’ ἘΝ" 

14 PRINT "* TO READ X,Y DATA FROM A TEXTFILE *" 

15 PRINT "* Pay 

16 FRINT FOE IO EI EI IIE IER IRAE IRIE! 

19 PRINT : PRINT 

20 FRINT "VERSION AS AT 1/08/86": PRINT 

50 PRINT : FRINT "POINTS MAY RE DELETED OR ADDED " 

60 D& = CHRE# (4) 

70 DIM A(352) ,B(352) 

80 PRINT : INPUT "ENTER FILENAME "36 

90 PRINT D#; "OPEN"; G* 

100 PRINT D#; "READ": GF 

110 Tel 

120 INPUT ACI),BCI) 

130 PRINT Δ(1),Β(1) 

140 IF ACI) = 9999 GOTO 170 

150 I= 1+ 1 

160 GOTO 120 

170 PRINT D$;"CLOSE" 

180 PRINT "END OF FILE" 

161 PRINT τ INPUT "LIST CONTENTS OF FILE ON PRINTE 
2? YORN 7? "3Z# 

162 IF z#< > “Y" THEN 190 

163 PR# 1 

184 PRINT : PRINT : PRINT “CONTENTS OF FILE ":G#: 

185 PRINT "AC(1)","BC(I)": PRINT 

186 FORN=1 T0T 

187 PRINT A(N),B(N) 

188 NEXT N 

189 PR# O 

190 PRINT τ INPUT "DO YOU WANT TO DELETE ANY POINT 
? "SF 

200 IF F# < >» "Y" GOTO 360 

202 PRINT : INPUT "BY POINT NUMBER OR VALUE 7? P OF 
2? "SFS 

204 IF FS < > "FR" THEN 210 

206 PRINT : INFUT "ENTER POINT NO. (REMEMBER THAT E 
H DELETION IS FOLLOWED BY REORDERING) > "sk 

207 IF K > rORK< = 90 THEN 205 

208 GOTO 280 

210 PRINT : PRINT "ENTER THE POINT TO BE DELETED : 
y" 
᾽ 

220 INPUT X,Y 


230° FOR J = 1 TOT 


IF δῷ) = X AND Β() = Y THEN 270 

NEXT J 

PRINT : PRINT "POINT NOT FOUND" 

GOTO 190 
κε 

FoR J =k +1 TOT 
αὖ - 1) = AGI) 
B(J - 1) = Bd) 

NEXT J 

jie a ce | 

PRINT + INPUT "DO YOU WANT TO DELETE ANY MORE PO 
INTS 7 "SHE 

IF Ht = "N" GOTO 360 

BOTO 202 

PRINT + INPUT "DO YOU WANT TO ADD ANY FOINTS Ὁ" 
εἰς 

IF L$ = "N" GOTO 560 

PRINT : PRINT “ENTER THE POINT TO BE ADDED : X,Y 
INPUT X,Y 

FOR J = 1 TO! 

IF ACI) > X GOTO 430 

NEXT J 
K = 


= K - 1 GOTO 500 
1) = ACJ) 

1) = BJ) 

= i! 

GoTo 450 


~ 
π 
ate une 


PRINT + INPUT "DO YOU WANT TO ADD ANY MORE FOINT 
gs 7 "sms 

IF ΜῈ = "N" GOTO 560 

BOTO τοῦ 

Sins Gate COW] | ee a 1.14) 4.4.» 

ἜΤ ΤΩΣ 

PRINT D&, "WRITE" SF 

FoR J = 1 TOT 

PRINT AGI): PRINT Β(5) 

NEXT J 

PRINT D#; "CLOSE"; G# 

PRINT + FRINT "CORRECTED DATA STORED IN FILE "4G 
+ 

PRINT + INPUT "DO YOU WANT A PRINTOUT OF MAX AND 
MIN VALUES ΤΥ OR N 7 "sC# 

Ir Ce = "NY BATO 790 

YMAN eT EMX π AGL) 
YMAX = B(1):MY = BCL) 

FOR a= 2 TOT 


690 IF Αἰ) = 9999 BOTN 745 

700 IF AJ) >» XMAX THEN XMAX = AGI) 

710 IF ACJ) “ MX THEN MX = ACI) 

720 IF B(J) >» YMAX THEN YMAX = BGI) 

730 IF Bid) < MY THEN MV = Rid) 

740 NEXT J 

745 #IF 7# = "N" THEN 750 

746 PR# 1 

750 PRINT : FRINT "XMAX = " )xXMAX 

760 PRINT "XMIN = ",MX 

770 PRINT "YMAX = ",YMAX 

780 PRINT "YMIN = ",MY 

705 PRO 

790 REM END OF ORIGINAL "READF ILE" 

B00 REM START OF ORIGINAL"FILEDISPLAY" 

1000 PRINT : INPUT "DO you WANT TO DISPLAY FLOT Ό 
SCREEN 7 Y OR N Ὁ "“s7# 

FOI. We ZS  ὝΣ THEN 1340 

1020 ΜῈ = G#: REM FILENAME 

1030 ΧΦ = XMAX:XO = Mx 

1035 RANGE = YMAX - MY 

1040 Y9 = YMAX + 0.0% * RANGE: YO = MY - 0.05 * RANE 

1050 W = I: REM NO. OF FOINTS 

1060 Di = X9 - XO 

1070 D2 = Y? - YO 

1080 PRINT =: INFUT "pO YOU WANT FOINTS JOINED 7? Y 
Ν Ὁ "sds 

1090 NORMAL 

1100 REM —OUTLINE FRAME 

1110 HGR 

1120 HCOLOR= 3 

1130 HPLOT 9,0 TO 279,0 TO 279,150 TO 9,150 TO 9! 

1140 D# = CHR (4) 

1150 PRINT D#; "OPEN"; MF 

11460 FOR 1 Ὁ 1. TOW 

1170 PRINT D#3 "READ": Mt 

1180 INPUT A 

1190 INPUT B 

1200 IF A = 9999 THEN 1320 

1210 REM -PLOTTING OF THE POINTS 

1220 X = INT (270 κ (A — XO) ἡ 51) +9 

1230 Y = 150 - INT (150 * (B - YO) / D2) 

ike hal x < 0 OR X > 279 OR Y 40 OR Y > 150 THEN 

1250 IF I= 10RV=#=0 ORU=O OR J#K > "Ὑ" TH 
HPLOT X,Y: GOTO 1270 

1260 HPLOT V,U TO X,Y 

1270 V=s X:U = Y 

1280 NEXT IT 

1290 PRINT DS: “CLOSE τς 

1300 PRINT 

1310 PRINT : GOTO 1330 


΄ 


1720 «PRINT + PRINT “END OF FILE REACHED ": GOTO 1290 


12320 GOSUB 1350 


1740 END 
1750 PRINT "IT TS POSSIBLE TO GET A FRINTOUT OF YOUR 
BRAPH " 


17460 FOR I = t TO 3000: NEXT I 

1270 PRINT + PRINT "DO YOU WANT TO ATTEMPT A FPRINTOU 
᾽ς τὴν ΘΕ Alin 

1380 INFUT Z#: IF ΖΦ { > "Ὑ" THEN 1570 

1290 PRINT : INVERSE ¢ INPUT "ENTER PRINTER TYPE : Ε 
PSON “Εἰς C-ITOH = C ":Y# 

1490 NORMAL. 

1440 IF Ye ὁ > "E" AND YS « > "C" THEN 1390 

1420 IF Ye = "E" THEN 1540 

1430 Dt = CHRF (4):E¢ = CHR# (27) 

1440 INVERSE : PRINT "REMEMBER APPLE PARALLEL INTERF 
ACE IN SLOT 1 ": NORMAL 

1450 FOR I = 1 TO 1000: NEXT I 

14460 PRINT DF: "PRL"? PRINT E#s "FP": REM FROPORT IONA 
ι. PRINT 

1470 PRINT D#_"FR#O” 

1480 PRINT D#;"BLOAD CITOH.ORJ" 

1490 CALL 6400 

1500 PRINT D#; "PRI" 

1510 PRINT E#:"N": REM RESET TO FICA 

1520 PRINT Df; "FRHO" 

1530 GOTO 1570. 

154 Ὁ PRH 1 

1550 PRINT CHR# (17) 

1560 PRH O 

1570 PR# 1 

1580 PRINT FRINT "DATA FOR GRAPH STORED AS ":Mt 

1590 PRINT + PRINT "ORIGIN AT X = ΧΟ AND VY = "SY 
ΠῚ 

1400 PRINT + FRINT "X RANGE 15. "gxXO3" TO "PX? 

14610 PRINT + PRINT "Y RANGE IS "εὐ" TO "3Y9 

14?0 rRe oo 

1470) RETURN 


B.4 Modification of the data 
(a) Baseline correction 


Program OPERFILE; When the baseline for a run is not flat, it is advisable to 
record the baseline separately, using program DSC. The baseline can then be 
subtracted from the sample data using program OPERFILE. The program is 
more general in that if allows addition, subtraction, multiplication or division 
of the contents of two files. 


ἘΠΕΒΣΌΜΕΡ. FEE ee er ee 
output. 

Output: Results of the operation stored in the file specified. Summ 
operations on the printer. 


10 HOME ᾿ 

“πὸ FRINT FO ORG IORI EIT δὲ 
30 PRINT "* PROGRAM ‘OPFRFILE' ΤΌ ADD, SURTRACT, 
40 FRINT "* MULTIPLY OR DIVIDE THE Y VALUES FROM 
SO PRINT "* TWO TEXTFILES 

55 PRINT errr rrr rrr rrr r rere tet... 2 2k eee 
546 PRINT : PRINT "VERSION AS AT 28/06/85 MER": Ρ 


57 FRINT "FILES MAY BE EDITED USING ‘READFILE '": 


" SMOOTHED USING ‘SMOOTHFILE’": FR 
ue OR DIFFERENTIATED USING ‘NUMDIFILE ° 
: PRINT " BEFOREHAND οἱ 


PRINT + PRINT 

60 Dt = CHR (4) 

65 PRINT : INFUT "IS PRINTER ON 7 Y OR N 7 "274 

66 IF 78 < > “Vv THEN 45 

70 DIM A(SO2) ,B(502) 

75 DIM C(500) ,D(500) ,AA(SO0) .Ε (500) 

80 PRINT : INFUT "ENTER FIRST FILENAME, NOTE ORDE 
F OPERATION IS FIRST+,-,*,/SECOND » "168 

90 PRINT D#; "OPEN"; Gt 

100 PRINT D#;"READ":GF 

110 I = 1 

120 INPUT ACI) ,B(T) 

120 PRINT ACI) ,B(T) 

140 IF ACI) = 9999 GOTO 165 

iso le I τ 

160 GOTO 120 

165 GM = 1-1 

170 PRINT D#;"CLOSE" 

180 PRINT "END OF FILE" 

190 PRINT + INPUT "ENTER SECOND FILENAME > "sk 

195 PRINT D#;"OFEN":H* 

200 PRINT D#;"READ"sHF 

210 I = 1 

220 INPUT C(I) ,D(1) 

230 PRINT C(I) ,D(1) 

240 IF C(I) = 9999 GOTA 245 

250 l= I+1 

260 GOTO 220 

265 HM = 1-1 

270 PRINT D#;"CLOSE" 

280 PRINT "END OF FILE" 


290 PRH 1 

z00 POKE 1657,80 

πιὸ PRINT : PRINT "FILE NO.1 = "2,05 

320 PRINT : PRINT "FIRST POINT IS "sA(1)3 SFC( 4)3B¢ 
1) 

330 PRINT = PRINT "LAST POINT IS "sA(GM)s SPC( 4)5B 
(GM) 

335 PRINT : PRINT "NO. OF FOINTS IN FILE,EXCLUDING 9 
999,0 WAS "3GM 

340 PRINT : FRINT ε PRINT 

250 PRINT : PRINT "FILE NO.2 = "3H 

πρὸ PRINT =: FRINT “FIRST POINT IS "sCC1)4 SFC( 4γε0 .( 
1) 

370) PRINT + PRINT "LAST POINT 15 "“:O(GM), SPCC 44D 
(GM) 

375 PRINT +: PRINT "NO. OF POINTS IN FILE,EXCLUDING 9 
999,00 WAS "3HM 

380 PRINT + PRINT - PRINT 

400 PR O 

410 PRINT + PRINT "CHOOSE THE OPERATION REQUIRED": PRINT 
"(1) ADDITION "s PRINT " (2) SUBTRACTION ": PRINT 
"(τ MULTIPLICATION "s PRINT "(4) DIVISION ": INPUT 
CH 

420 IF CH < 1 OR CH > 4 THEN 410 

450 PRINT : INPUT "ENTER NUMBER OF STARTING POINT IN 
FILE 1 > ΟΣ 

460 PRINT τ INPUT "ENTER NUMBER OF STARTING POINT IN 
FILE 2 > "sH1 

465 GO = GisHO = Hi 

470 PRINT =: INFUT "ENTER NO. OF FOINTS TO ΒΕ TREATED 

> ΝΡ 

490 FOR I = 1 TO NF 

500 ON CH GOTO 510,520,530 

505 GOTO 540 


S10 AACI) = BCBL) + δ(Η1): REM ADDITION 
515 GOTO S50 
520 AACI) = R(GL) - D(Hi): REM SUBTRACTION 
\, 575 8nTO 550 
he «πὸ AACT) = BIG1) κα D(H1i): REM MULTIPLICATION 
Ν 535 GOTO 550 
540 AACT) 5 Β(Θ1) ἡ D(Hi): REM DIVISION 


sso i = GO+ I= i 
S52 F(T) = ACK) 
Gi = Gi + 1:Ht = H1 τὶ 


Es 
u 
a 
> 


Se 
‘. 7 


TN 

Sh Sos PRINT F(t): SPCC 4);AACT) 

i {i 540 NEXT T 
As ἃ ceo PRINT + INPUT "DO YOU WANT TO STORE OUTPUT IN A 


Vy FILE? Y OR N 7"5C# 
{ 563 IF ΟΣ = “Ν" THEN 640 
565 PRINT + TNPUT "ENTER A FILENAME FOR STORAGE Ὁ " 
ΓΕ 
PRINT Df: "OREN" l¢ 


= 
ΡΝ 
τὸ 


ies 


= 
a 
Ἢ 


ταν 


SEU MLN Fx τὰν νος 8h 
590 FOR I = 1 TO NF 
600 PRINT F(1): PRINT AAT) 
610 NEXT I 
615 PRINT 9999: PRINT ὁ 
620 PRINT D#;"CLOSE";L* 
630 PRINT : PRINT "DATA STORED IN FILE "3l# 
640 PRINT : INFUT "DO YOU WANT A FRINTOUT OF MAX & 
MIN VALUES 7 Y OR N 7 "τὸς 
650 IF C# = "N" GOTO 790 
660 XMAX = F(1)2MX = FCL) 
670 YMAX = AACL)sMY = AACL) 
680 FOR J = 2 TO NF 
690 IF F(a) " 9999 GOTO 75ὸ 
700 IF Ε(Ὁ) >» XMAX THEN XMAK © Fd) 
710 IF F(J) < MX THEN MX = Fad) 
720 IF AACI) > YMAX THEN YMAX = AACI) 
730 IF AACI) < MY THEN MY = AAG) 
740 NEXT J 
750 PRINT : PRINT "XMAX = " SXMAX 
760 PRINT "XMIN = ",MX 
770 PRINT "YMAX = ",YMAX 
780 PRINT "YMIN = ",MY 
790 PRINT + INPUT "DO YOU WANT A SLUIMMARY ON THF FI 
TER 7 N.B. SWITCH ON ! YORN ? "3C# 
ΘΟ IF C#< > "Y" THEN 1000 
g02 PRINT : INPUT " FULL=F, SHORT#S > "3Q* 
B10 PR# 1 
820 POKE 1657,80 
830 PRINT : FRINT : PRINT “GUMMARY OF OPERATION O 
ILES ";G#;" AND "SH 
840 PRINT 1 gg τα ἐσ πα πὸ πε τα τ τα a ee τν τὸ μα τπ ον τα απ τ πα τε τα δε με φετο πε πραγ εττααν 
mescosen': PRINT 
850 PRINT "THE FILES WERE "3 
860 ON CH GOTO 880,890,900 
870 GOTO 9710 
@80 PRINT "ADDED ": PRINT : GOTO 920 
890 PRINT "SUBTRACTED "+ PRINT : GOTO 920 
900 PRINT "MULTIPLIED "+ FRINT : GOTO 920 
910 PRINT "DIVIDED ": PRINT 
920 PRINT “ORDER WAS FILE "sGH;" +,-,% OR / "3" F 
"SH#: PRINT 
930 PRINT "STARTING FOINTS WERE POINT "3G0;" IN 
"3G," AND POINT "sHO;” IN FILE ";H#: PRINT 
NPs" POINTS WERE TREATED ": PRINT 
940 PRINT : PRINT 
945 IF Of = "S" THEN 990 
950 FOR I = 1 TO NP 
960 PRINT F(T)3 SFC 4) AAT) 
970 NEXT I 
990 PRINT : PRINT "XMAX = ",XMAX 
991 PRINT "XMIN = ",MX 


n 


ΘΕΦΕΘ 


Bs 


} 


-- 
ς 


9992 FRINT "YMAX = ",YMAX 

99% PRINT "YMIN = ",MY 

994 IF τὰ = “ " THEN 999 

995 FRINI + PRINT "DATA STORED IN FILE "sla PRINT 
997 FRAO 

«οὐ END 


(b) Smoothing of data 


There are many approaches to this problem (chapter 12 and Fig. 12.3). Three 
programs are given below for smoothing of noisy data, using different 
algorithms. The algorithms are given in the references quoted in the listings. 


Programs SMOOTH FILE, BINARY FILTER, ROBUST FILTER 


Input: Filenames and adjustable parameters as specified. 
Output: Smoothed data ina file, summary of operations on the printer. 


10 HOME 

20 PRINT ΠΤ IEE EEE EEE IIIA IEE 
πὸ ΡΗΙΝΤ "* DATA SMOOTHING FROGRAM ae 
AQ) PRINT "* BASED ON SAVITZKY GOLAY FILTER i 
50 PRINT "* REF: ANAL. CHEM. ,36 (1964) 1627 *" 


60 PRINT "* USING 5 POINT CONVOLUTION TABLE I ΝΣ 
7O PRINT GEIB HIEIHIESE ΤΥ ΣΝ δ με FEIT 
@O PRINT : PRINT 

90 De = CHRE (4) 

100 DIM Υ (250) ,Z (250) ,D(5) ,C (5) 

110 I = 0: REM COUNTER 

120 PRINT "DATA MUST BE AT EQUALLY-SPACED INTERVALS" 


130 PRINT 

140 PRINT + PRINT “DATA MUST END WITH 9999,0 " 

{so PRINT + INPUT "ENTER TIME INTERVAL IN MINUTES "; 
x 

160 PRINT : INPUT "ENTER NAME OF DATAFILE "1:65 

170 PRINT D#; "OPEN" 3G 

190 PRINT D#; "READ"; G* 

190 FOR I = 1 TO 350 

200 INPUT T 

210 INPUT Y(T) 

220 IF T = 9999 THEN 240 

230 NEXT I 

240 PRINT D#"CLOSE"s 6 

250 PRINT + PRINT " CALCULATING — PLEASE WAIT” 

260 M= I - 1: REM NO OF POINTS ENTERED 

270. ΠῈΜ - CONVOLUTING FACTORS FOR 5 POINT SMOOTHING 

280 c(i) = - :0(2) = 122003) = 17:π|4) = 12:C0(5) = 


270 N= οὺυὲ pS PEEL STE τον 

300 REM - READ 5 Y VALUES INTO TEMPORARY STORE 

310 K = O: REM - FOINT COUNTER 

320 S = O: REM - SUMMATION 

330 REM - MULTIPLY Y VALUES BY CONVOLUTING FACTORS, 
UM AND THEN NORMALISE TO GET NEW VALUE OF THE MI 
-FOINT 

240 FOR J =1 705 

350 DJ) = Y(K + J) 

360 D(J) = CJ) * Dd) 

370 S = § + Did) 

380 NEXT J 

3290 IF K > 5 THEN 430 

400 REM - STORE FIRST TWO POINTS 

410 Z(1) = YC1) 

420 7(2) = Y(2) 

430 7K +3) =S/N 

440 Kak + 1 

45. IF K< = (M—- 5) THEN 320 

460 Z(M - 1) = Y(M - 1 

470 Z«(M) = Y(M) 

480 REM - LAST TWO FOINTS STORED 

490 PRINT +: INPUT "ENTER FILENAME FOR SMOOTHED VAL 
S “sks 

500 PRINT Df; "OPEN" Ss FS 

510 PRINT D#;"WRITE"SFF 

520 FORL = 1 TOM 

530 FRINT L * yoy att ED) 

540 NEXT L 

550 FRINT "9999": PRINT "oO" 

5460 PRINT Dt; "CLOSE" ;F# 

570 PRINT : INFUT "pA YOU WANT MORE SMOOTHING ΤΥ 


N ? "3B 
seo IF ΒΦ = "Y" THEN 400 
590 END 


600 FOR L = 1 TO M 

610 Y(L) = Z(L) 

620 NEXT L 

630 REM - SMOOTHED VALUES NOW RECOME INFUT 
640 GOTO 310 


10 HOME 

20 PRINT ΝΗ ΝΗ ΝΝΝΑΥΝΚΝΚΚΝΝΝΝ ΚΝ IEEE TEEIIEE GREET 
30 PRINT "* DATA SMOOTHING PROGRAM * 
40 FRINT "* BASED ON RINOMIAL FILTER * 
πὸ PRINT "" REF? REV. SCT. INST. ,54(8) (1984) 1034 * 
60 PRINT "* MARCHAND AND MARMET * 


70 PRINT ΤΥ GIBEEIEEEEE AHI EGEE 


—_ 
ΞΟ 

( i 

eee - 
a 


τς 
of Y 


ΞΟ; 


x 


HOD! 
—s 


TOY OX 
τοῦς πως 


ss 


— 


pee 


A SS 
fs) 
Ἐς ὦ 


ΘῈ 


EB, 


ree OT 


75 PRINT : PRINT "VERSION 7/5/86" 

AO PRINT + PRINT 

90 Dt = CHRt (4) 

100 DIM Y(SOO),Z(500) ,X (500) 

410 1 = Of REM COUNTER 

115 J πῶς REM PASS COUNTER 

{20 PRINT "DATA MUST BE AT APPROXIMATELY EQUALLY-SPA 
CED INTERVALS" 

τὸ FRINT 

140. PRINT + PRINT "DATA MUST END WITH 9999,0 " 

160 PRINT + INPUT "ENTER NAME OF DATAFILE "0: 

170 PRINT D#;"OPEN"3G6¢ 

190 PRINT D#;"READ";G¢ 

190 FOR I = 1 TO 500 

200 INPUT ΧΕΙ) 

210. INPUT Y¢T) 

220 «IF ΧΙ) = 9999 THEN 240 

230 «NEXT T 

240. PRINT D#;"CLOSE";G$ 

mas PRINT : INPUT "SET NO. OF PASSES EG. 25 ":Ν 

250 PRINT : PRINT " CALCULATING - PLEASE WAIT" 

340 M = I - it REM’ NO OF POINTS ENTERED 

265 Je 9 +1 

270 FOR I = 1 TO (M - 1) 

BB0 Z(1) = (YC) + VCT + 19) 7 2 

290 NEXT I 

300 FOR I = 2 TO (M - 1) 

Z10 VOI) = (200 = 1) + 20D) 7 2 

πρὸ NEXT I 

320 IF J < = N THEN 265 

ΧΕ PRINT : INFUT "OUTPUT ON SCREEN (S) , PRINTER (F 
) OR ONLY TO DISK (D) 7 "3Z# 

πρὸ IF 7# = "D" THEN 490 

395 IF 7# = "S" THEN 400 

ποθ PRH 1 

aoo FORL = 1 TOM 

aio. PRINT X(L),Y(L) 

420 NEXT L. 

A25 PRH Ὁ 

ago PRINT + INPUT "ENTER FILENAME FOR SMOOTHED VALUE 
S "1 

500 PRINT D#: "OPEN"; FS 

S10 PRINT D#;"WRITE"S FS 

520 FORL = 1 TOM 

530 PRINT X(L): PRINT Y(L) 

540 NEXT L 

550 PRINT "9999": PRINT "ὁ" 

πρὸ PRINT Df: "CLOSE"sFS 

590 END 


aw 
20 


τὸ 


35 


110 
120 


130 
140 
150 
160 
170 
1860 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 


re 


FRINT TEI IEEE AEE με ΕῚ 


FRINT "* PROGRAM "ROBUST FILTER’ “" 
PRINT : FRINT SU GEE EM 
κπ" 
PRINT : PRINT "VERSION AS AT 3/07/85 MWE/MEB ed 
PRINT : PRINT "BASED ON B-M.RUSSTAN AND Ww. HARDL 
J. APPL. SPECTROSC. ,28(1984) 309-313 ": PRINT 
PRINT "SMOOTHING IS TUNED BY TWO PARAMETERS "=: 
"4. SIGNAL-TO-NOISE RATIO": FRINT "2. HALF-WIDTH 
F DATA PEAKS " 
PRINT "THESE PARAMETERS HAVE TO BE SET IN LINE 
(aS 5 AND P ) ACCORDING TO THE SCALE OF THF Ih 
T DATA AND RY TRIAL AND ERROR " 
S§ = 1:P = 50 
INVERSE 
PRINT "PRESENT SETTING + § = "ge" P= "SP 
NORMAL 
PRINT : FRINT "Δ΄ WINDOW SI7E ° IS NOW SET ad] 


REM — WINDOW DIMENSION 

INPUT "INFUT WINDOW HEIGHT (E.G. 50)" 

INFUT "INPUT WINDOW WIDTH (ODD NO. 3-9 }"1Ὲ 
DIM X (500) «Υ (501) ,W(23) ,2 (25) 

PRINT + INFUT "ENTER FILENAME FROM WHICH DATA 
TO BE READ > "3A 

PRINT : INPUT "ENTER FILENAME IN WHICH DATA 
TO BE STORED >» "BF 

PRINT : FLASH +: PRINT “READING DATA FROM DISK 
":; NORMAL 
D# = CHR# (4) 

PRINT D#; "OPEN": At 

PRINT Df: "READ"; At 

FOR I = 1 TO S00 

INPUT X(T): INFUT Y(T) 

te Xr) = 9999 THEN Zoo 

NEXT T 
Ges t= 1 

PRINT δες "CLOSE": At 

HOME 

GOSUB 450 

FOR N = INT (EF ἡ 3) TO 9 - INT (E ἡ 2) 
GOSUB SSO 
Y(N) = UNT (ΤΊ 

PRINT X(N) ,Y CN) 

NEXT N 

PRINT +: FLASH : FRINT "OPENING FILE "3h: NOF 
FRINT D#:; "OPEN": BE 

PRINT D#, "WRITE" S RF 

FOR Σ = 2 ΤΟΝ = 14 

PRINT X(T): PRINT Υ{Τ) 


= 
CY 


Ἑ 


OHO 
=e 


350 


340 
350 
τοῦ 
270 


πρὸ 
πὸ 
Aw 
aio 


420 


42 
427 
το 
440 
450 
460 
470 
480 
ag) 
500 
510 


s20 


NEXT T 

PRINT 9999: PRINT 9999 

FRINT Df; "CLOSE"; Bf 

PRINT + PRINT "DO YOU WANT A SUMMARY ON PRINTER 
> Y OR N 7 "es INPUT 7: 

Ir 7# ¢ > "Ὑ" THEN 410 

ree 1 

FOVE 1657,80 

PRINT +: PRINT "SUMMARY OF OFERATTONS": PRINT "== 
neccennenecorecmanacsmm=': PRINT 

PRINT I - t3" DATA FOINTS FROM FIIE "SAG: PRINT 
“GMOOTHED USING ROBUST FILTER"? PRINT "OUTPUT STO 
RED IN FILE "sRt 

PRINT + PRINT "WINDOW HEIGHT SET AT "sk 

PRINT "WINDOW WIDTH SET AT "3E 

FRINT : PRINT 

PRO 

END 

REM 

REM - WINDOW SUBROUT INE 
f=, INT 7 24 1) 

FOR r= 1TOM~ 1 

X=S* I/M 
We 0.75 # (ho X κ X) 
wim + T) = WeW(M — ID #W 


NEXT I 
» WIM) = 0.75 
RETURN 
REM 
REM - SMOOTHING SUBROUT INE 
L=N- INT (Ε / 2) 


FOR I = 1 TOE 

2. 1, tgs ACL et ea 4) 
NEXT T 
T= Ύ Νὴ πὸ 
Θ1 = 0:82 πὸ 

FOR I= 1 TOF 
Viena Zh a= 
wom wr) 
Y2 = Κιγξ "ὦ 

IF Y1 > K THEN 710 
y2= -κ 

IF Yi ¢ - THEN 710 
v2 = νιν = 1 
81 =~ Si + (ὦ κ Υ2) 
52 = 82 + ὦ κ YS) 
NEXT I 

IF 82 = 0 THEN 52 = 0.75 
H = Si 7 52 

IF ( ARS (H) < = FP) THEN RETURN 
TeT+H 

BOTO 420 


(c) Sealing of data 


It is sometimes useful, c.g. in kinetic analyses. ἴο be able to scale either the 
the Y data, or both. 


Program SCALEFILE 


Input: The filename, the X and Y scale factors and the Y base corrects 


O 


offset. A filename for storage of the scaled data. 


utput: Scaled data in a file. 

HOME 

PRINT SORES IG IIIS TAHITI IE 
PRINT "* PROGRAM "SGCALEFILE’ TO MULTIPLY * 
PRINT "* THE X AND Y VALUES FROM A TEXTFILE - 
ΡΕΙΝΤ "* BY CONSTANT FACTORS ” 


PRINT RHEE ROI OIEE EEE IIB RIIB EES # 
PRINT : PRINT "VERSION AS atTi9/07/985 MER": FRI 


D# = CHR# (4) 
DIM A(SOO) ,B(500) 
DIM τ (500) ,D(500) ,AA(SOO) ,F (500) 
PRINT : INPUT "ENTER FILENAME > "36% 
PRINT + INPUT "ENTER X SCALE FACTOR > "3XF 
PRINT + INPUT "ENTER Y BASE CORRECTION FACTOR 
"3 VB 
PRINT : INPUT “ENTER Y SCALE FACTOR > "3 YF 
PRINT : INVERSE - PRINT "PLEASE WAIT , READING 
NORMAL =: FRINT 
PRINT Df; "OPEN"; G* 
PRINT D#; "READ"; GF 
im 1 
INPUT A(I),B(T) 
IF ACI) = 9999 GOTO 145 


ACT) = ACI) * XF 
RB(I) = (BCI) -- YR) * YF 
τα τὰ 
GaTO 120 
GM = I- 1 


PRINT D#; "CLOSE" 

PRINT "END OF FILE" 

PRINT t INPUT "ENTER A FILENAME FOR STORAGE 
sls 

PRINT : INVERSE =: FRINT “PLEASE WAIT , WRITIN 
: NORMAL : PRINT 

PRINT D#; "OPEN"3Lf 

PRINT D#3 "WRITE" LF 

FOR I = 1 TO GM 

PRINT A(I): PRINT R(T) 

NEXT I 

PRINT 9999: PRINT ὦ 


196 THermat anaryany suremurs 


620 PRINT D¥s "CLOSE" ει 4 
420) ΓΕΙΝῚ - PRINT "DATA STORED IN FILE "sl* 
τοῦ PND 


BS Processing of the data 
(a) Numerical differentiation 


The DSC curve is already an example ofa differential measu rement and so is not 
usually differentiated further. (fan integral measurement such asa TG curve has 
been recorded, then a DTG curve can be obtained by numerical differentiation. 
The usual method is that of Savitsky and Golay (reference given in the listing of 
program NUMDIFILE). 


Program NU MDIFILE 


Input: Option of entering data manually or from a file. The X data must be 
equally spaced. Option for preliminary smoothing as in program 
SMOOTHFILE, and then differentiation may be based on either (1) 
a quadratic or (2) ἃ cubic polynomial. 

Output: Tables of values on printer and storage of values in a file. 


10 HOME 
20) PRINT wa eg Ha ERIE ATI RIOR II IE” 


δὸ FRINT "* NUMERTCAL DIFFERENTIATION PROGRAM mt 


AQ FRINT "* RASED ON SAVITZKY AND GOLAY *" 
πο ΓΠΙΝΤ "* PERL ANAL. CHEM. 536 (1964) 1627 ἣν" 
60 PRINT "» USING & POINT CONVOLUTION TABLE ia 


70 PRINT τ UL GEIS ΕΝ IE EEEEBEOESEEIE EI TEE” 
AO PRINT : FRINT 

90 DIM y (950) «7 (250) .Ὁ (5) ,O(5) 

100 PRINT “YOULL MAY TNPUT DATA MANUALLY (M) (J 

410 PRINT "OR FROM A DATA FILE(F) " 

120 INPUT "ENTER M OR F "sO 


130 IF Of ¢ . "EF" THEN 340 
140 PRINT = INPUT “ENTER FILENAME "sF# 
150 DF ~ CHR (a) 


160 PRINT D#; "OREN" P# 

170 PRINT DF; "READ" ΡῈ 

160 FOR I = 1, TOS 

190 INPUT X(T) 

poo INFUT Υ(Τ) 

210 NEXT I 

220 IF (x(2) - XC) = (ees) - ΧΟ) Σ 0.00001 THEN 
B10 

πὸ FOR I = 3 TO 250 

240 INPUT C 

250 INFLUT YT) 

260. ρῶν 9999 THEN 280 

270 «NEXT I 

2eac X = XC?) x (1) 


290 
300 
310 
320 


330 
340 


350 
360 
370 
780 


390 
400 
410 
420 
430 
440 
450 
460 


470 
480 
490 
500 
510 


520 
530 
540 
550 
560 


570 
580 
590 


600 


610 
620 
630 


640 


650 
660 
670 
680 


PRINT D#; "CLOSE": PF 

GOTO 450 

PRINT D#; "CLOSE"; F* 

PRINT : PRINT "X DATA MUST BE AT EQUALLY SFACED 
INTERVALS ! " 

GOTO 100 

PRINT "DATA MUST RE AT EQUALLY~SPACED INTERVALS’ 
PRINT 

INPUT "ENTER X INTERVAL "3X 

PRINT 

PRINT "Y VALUES (MAX NO=250) SHOULD NOW RE ENTFRI 


D IN SEQUENCE" 

PRINT “END WITH 9999" 

FOR I = 1 TO 250 

INPUT Y(I) 

IF YC(I) = 9999 THEN 450 

NEXT I 

PRINT "MAX MO OF Y VALUES REACHED " 
M= 1 —- 1: REM NO OF POINTS ENTERED 


PRINT τ PRINT "OPTIONS AVAILABLE +: ": PRINT i 
(1) STRAIGHT TO NUMER ICAL DIFFERENTIATION": FRIN 
᾿ς (2) SMOOTHING FIRST" 

INPUT "ENTER 1 OR 2 "30 

IF @ = 1 THEN 510 

IF @ = 2 THEN 550 

GOTO 460 

PRINT : PRINT "FURTHER OPTIONS τὴς PRINT " a 
QUADRATIC FOLYNOMTAL”: FRINT " (2) CUBIC POLYNC 
TAL" 

INPUT "ENTER 1 OR 2 ay 

IF V = 2 THEN 6350 

GOTO 590 

REM - CONVOLUTING FACTORS FOR 5 POINT SMOOTHING 
cq) = - 3:€(2) = L22E(S) «1710 (Ά}) = 1210 (5) " 
- 3 
N = 35: REM - NORMALISING FACTOR 

GOTO 4460 

REM - CONVOLUTING FACTORS FOR 5 PT NUMERICAL D 


FERENTIATION (QUADRATIC) TABLE ITI 


CG) = ~ 2b = = 1:C(3) = O20(4) 5 1:05) 
2 

N = 10: REM - NORMALISING FACTOR 

GOTO 660 

REM - CONVOLUTING FACTORS FOR 5 PT NUMERICAL D 


FERENTIATION (CUBIC) TABLE IV 

CY = 10 2). > @:c(3) = o:c(4) = 8:0 (5) = 
1 

N = 12: REM — NORMALISING FACTOR 

REM - READ 5 Y VALUES INTO TEMPORARY STORE 
FRINT : PRINT "CALCULATING - PLEASE WAIT " 


kK = O: REM - POINT COUNTER 


usu 


690 S πῶς: REM - SUMMATION 

700 REM - MULTIPLY Y. VALUES BY CONVOLUT ING FACTORS,S 
UM AND THEN NORMALISE TO GET NEW VALUE OF THE MID 
-POINT 

710 FOR J =1 705 

720 NII) = YCK + J) 

730 DCJ) = CUI) * DJ) 

7405285 + DtJ) 

750 NEXT J 

7640 IF K > S THEN 820 

770 REM — STORE FIRST TWO POINTS 

7e0 IF Q< > 2 THEN BOO 

790 X wit REM - VALUE OF X INTERVAL NOT REQUIRED IN 

SMOOTHING .SEE | INF. B90 

poo Z(1) = YL) 

p10 Z(2) = Y(2) 

e20 ΖΚ + 3) = SS (x κ N) 

aso k =k + 1 

δ IF <= (M - 5) THEN 690 

eso ΖΜ — 1) = V«(M = 1) 

@60 Z(M) = Υ(Μ) 

870 REM - LAST TWO POINTS STORED 

B80 REM - OUTPUT FOLLOWS 

e890 FRINT 2 INPUT “pO YOU WANT A PRINTOUT 7 "sd 

900 IF J# = "N" GOTO 1150 

910 PRINT : PRINT "IS PRINTER ON IY BRING fo 

920 GET AF 

920 IF ΔῈ <¢ > "Vv" THEN 910 

940 FRY 1 

950 FRINT CHRE (137) 3 "BON" 

960 IF @ = 2 THEN 1050 

970 IF V = 2 THEN 990 

980 7t = " (QUADRATIC APPROX) ": GOTO 1000 

φῳφὸ 7% = " (CUBIC APFROX) " 

1000 PRINT + PRINT uy FIRST DIFFERENTIAL S":7# 

1070 PRINT πα em et air TOT eS STURN E a=E ” 


joDo PRINT τ PRINT "X","¥","DY/DX" 


pou PRINT " τ ΙΑ, pelts artes eed crccers anger lar Neuse ee eat ee c= paris 
1040 PRINT : GOTO 1110 

1OUO PRINT ς: FRINT " SMNOIHED VAL WES? 
1oae PRINT x gesneneerersrcemscctreremen!! 

170 PRINT 


1000 FRINT " xX " "ORTG Y","SMOD THED Vv" 
19m PRONT " sem essa ee oe mane “Ὁ μενον eens tipo mm meew mene 


1100 PRINT 
4410 FOR tL = 1 TOM 
112% PRINT 1. * Χο Gees Ge) 


1130 NEXT lt. 

1140 PR ὦ 

1150 PRINT 

1140 INFUT "DO YOU WANT TO STORE THIS OUTPUT IN A I 
TAFILE 5. "30F 

τ 70 “LF OS κα >» "Y" THEN 1350 

1180 PRINT = INPUT "ENTER FILENAME ἡ: ΜῈ 

1190 Df = CHR (4) 

1200 PRINT DF: "OPEN EME 

1210 PRINT D#; "WRITE sit 

1220 IF @ < » 2 THEN 1250 

1230 Bi = 1:82 - Μ 

1740 GOTO 12460 


1270 A=L* X 


4 1000022 = INT (Z)2Z7(L) * Z / 10000 
DECIMAL FLACES 

PRINT Z(1L) 

NEXT ι. 

PRINT 9999 

PRINT ὁ 

PRINT Dt: "CLOSE": ΜῈ 

PRINT : PRINT "FINISHED 7 YOR N 7" 
GET BS 

IF Be = "Y" THEN 1440 

IF « » 2 THEN 440 

FOR L = 1 TOM 

1400 Y(L) = Z(L) 

1410 NEXT L 

1420 REM — SMOOTHED VALUES NOW BECOME INFUT 
1430 GOTO 460 

1440 END 


(b) Peak integration 


Peak area determination is usually done by numerical integration using either 
Simpson's rule or the trapezoidal rule. For a large number of closely-spaced 
points, the trapezoidal rule is adequate and simpler. Limits of integration have 
to be set after examination of the data using programs FASTDISP or 
READISP, Note that the X values may he unevenly spaced. 


Program TRAPSUM 


Input: Filename for data, starting and stopping values of X. 
Output: Results on printer, i.e. area at return, area correction, and corrected 
aren as defined in Fig. Β.}. 


at 


ener" 


x 


shaded region = area at return = area AEDA -- area OCBD 

triangle ABC = area correction = area ADBA + area DCBD 

corrected area = shaded region + correction 

area AEDBA = area AEDA ~ area DCBD + area ADBA + area DCBD 
- area AEDA + area ADBA 


Figure B.1 Peak area determination, 


1. HOME 

2 PRINT SEE IEIE IESE IEEE IRIE 

= PRINT "* PROGRAM TRAPSUM FOR AREA INTEGRATION ἊΨ 

4 PRINT “* USING THE TRAFEZOIDAL RULE x" 

5s PRINT HOGER FEEEHEEIEEEA HEIGHT! 

ὁ PRINT ¢: PRINT “YERSION AS AT 15/7/85 MER" 

7 PRINT : PRINT 

8 PRINT "THIS VERSION IS FOR ": FLASH : PRINT "UINEOU 
AL": NORMAL +: PRINT "X INTERVALS ": PRINT 

9 PRINT "THE PRINTER MUST BE ON ": PRINT " 


eeomernere':s PRINT 


10 INVERSE : PRINT "DO YOU WANT CUMULATIVE ARRA rr 
Tout 7 Y OR N 2": INPUT Bt: NORMAL 

“a, ΤΕ ΒΕ Ὁ ΣῪ" THEN 29 

12 PR# 1 

1% POKE 1657,80 

14 PRINT "PEAK INTEGRATION ΒΥ TRAPEZOIDAL RULE " 


16 PRINT Vpeeere deere στο το πα τα το τε πε sesame Ν᾽ 


18 PRINT 
20 PRINT “POINT X VALUES Y VALUES Y INC Cc 
ULATIVE AREA " 
oy) PRINT Wereee = Soe ἀπο εεινινεροθησασέ κε εν τυ λΡ πε τ ΡΥ ΘΕ Υ ταν 
"“" 
25 ΡΒ O 


Φ PRINT : PRINT 
πὸ INPUT "ENTER FILENAME CONTAINING DATA > “sAt 
90 PRINT 


100 INFUT “ENTER STARTING X > EF 
110 PRINT 
120 INPUT "ENTER STOPPING X > “ss 
130 PRINT 


270 F = O15 = 0 

240 Dt = CHR# (4) 

249 FRINT D#£" OPEN": AF 

250 PRINT D#; "READ"; At 

260 INPUT X1,Y1 

265 PRINT X1,Y1 

270 1 = 1 

280 INPUT X2,Y2 

290 PRINT X2,Y2 

soo Ir =I+ 1 

310 IF X2 = 9999 THEN FRINT "END OF FILE DETECTE 
: GOTO 2000 

315 IF F = 1 THEN 500 

S20) IF X2 > = PF THEN 340 

255 GOTO 280 

360 F = i:l = Y2:6 = I:M = X2 

370 GOTO S520 

ποῦ T= (Y2 - Ll) + (x2 = XP) 

soz S=S+T 

504 IF BF ¢ > "y" THEN 518 

505 PRH I 

510 PRINT Is" "1 X23 Weg 
Υ στ." "3S 

515 PRY O 

516 IF X2 > = @ THEN 400 

520 X1 = X2 

520 Yi = Y2 

545 GOTO 280 

600 8 = S§ - (((v2 - L) 5 «M2 — X1)) / 2) 

601 REM END CORRECTION 

610 FR = ¢¥2 — 0) ἢ (X2 - M) 1 2 

611 REM AREA CORRECTION 


ave 


615 PRH 1 
620 PRINT + FRINT 
425 PRINT “OUTPUT FROM TRAPSUMS": PRINT "<--eomeree= 


“παπποῦ PRINT 
πὸ PRINT "AREA AT RETURN = "3S 
640 PRINT "AREA CORRECTION = 
450 PRINT “CORRECTED AREA = "3S ~ R 
ADO PRINT * eer rns em nS VN es et |" 
291: PRINT 
ἡ 692 PRINT "DATA FILENAME - "08 
t 494 PRINT : PRINT "INTEGRATED RETWEEN X = "3F3" AND 
a y= "30 
ὶ 


cheat a 


rr, 
~ 
=-{ 


a3 


oe 


χοῦ PRA O 
2000 PRINT D+: "CLOSE" 
2100 END 


Ξ, + er 


πρμεπιννν. 


Explanation of the symbols used in tl 
text 


(SEE INDIVIDUAL SECTIONS FOR MOR 
DETAIL) 


=area; amplitude of vibration: Arrhenius pre-exponential factor: c« 
ponent of binary mixture 

=calibration factor (heat capacity): bulk modulus; sample thicknes 
component of binary mixture 

=heat capacity; calibration constant; Curie constant 

=thermal diffusivity 

= Young's modulus; emanating power, Arrhenius activation energ) 

=force: fraction melted 

=shear modulus 

=enthalpy; magnetic field strength 

=instrument constant; heat transfer coefficient 

=crysoscopic constant 

=length 

=magnetization; molar mass 

=period of oscillation 

=heat 

=gas constant. thermal resistance: recorder sensitivity 

=entropy; chart speed; area 

= temperature 

=volume 

=general ordinate 


> 


πο Ὁ ΞΟ ΞΟ ΜΠ ΘΟ ὦ 


=constant; partial area: activily 

=constant; temperature integral (section 13.5) 
=specific heat capacity 

= diameter 


urs 


ς exponential base 
1{() ~funetion 
“(0 —funetion 


᾿ 
᾿ (g) pas 
᾿ " h() =funetion 
ike yy h  —baseline displacement; height: relative amplitude of vibration 
’ ‘4 k — =rate coefficient 
Oi: (i) =liquid 
Ἔ ΝῊ "  =mass; constant exponent in rate equation 
ἢ "  =apparent order of reaction; number of moles 
i ἱ Ρ = pressure, constant exponent in rate equation 
a p(x) integral defined in section 13.5 (x= F/RT) 
4 q = =heat 
+, Κη ς  =partial area 
) ἢ (s) =solid 
14 1  =time 
hi χ =mole fraction; dimension 
{ x, =mole fraction of solvent 
x, =mole fraction of solute 
xt =make fraction of impurity 
y  =dimension 
z  —=dimension 


Δ =change or finite difference, damping 

a  =coeflicient of linear expansion: fractional reaction 
β  =see chapter 4 

ε =strain; dielectric constant; undetermined premelting 
n = viscosity coefficient 

0 =angle; Weiss constant 

@(7) =function of temperature 

κ =volume magnetic susceptibility 

λ =thermal conductivity 


pe =chemical potential 

» 5 Poisson's ratio 

\, = Fresonance frequency 

p  =density 

a τεϑίγοος 

@  =heating rate 

y  =mMass magnetic susceptibility 


Subscripts 
c  =calibrant 
f —furnace: fusion 


2 ae ase 


= glass transiuon 

= geometric mean 
=melling 

=initial 

=constant pressure 
= reference 

= sample 


Index 


Activation energy 20, 60, 101, 127, 
145, 149 

ADC see Analogue-to-digital converters 

Alloys 42, 60, 61 

Analogue computers 117 

Analogue-to-digital converters (ADC) 
119-21, 174 

Annealing 111 

Antiferromagnetism 103 

Area under the curve 23, 33, 53, 
157, 158, 171, 200 

Arrhenius equation 60, 127, 129, 136, 
137 

Atmosphere control 8, 12, 14, 19. 24, 
52, 91, 106, 150 

Atmosphere, self-generated 13 

Automation of instruments 120, 122-5 

Avrami-Erofe’ev equation 131, 135. 
145 


122. 


Baseline 30, 32, 37, 38, 121, 174, 186 

Binary coded decimal (BCD) output 
120 

Boersma DTA see Calorimetric DTA 

Borchardt and Daniels kinetic method 
136-8 

Bulk modulus 70 

Buoyancy effects 12, 122 


Calibration 27, 33, 34, 38. 57. 63. 66, 
85, 91, 113, 171 
Calorimetric DTA (or heat-flux DSC) 
25, 28, 29, 34, 38. 50 
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